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Abstract: We present a new perturbative formulation of non-equilibrium thermal field 
theory, based upon non-homogeneous free propagators and time-dependent vertices. Our 
approach to non-equilibrium dynamics yields time-dependent diagrammatic perturbation 
series that are free of pinch singularities, without the need to resort to quasi-particle ap- 
proximation or effective resummations of finite widths. In our formalism, the avoidance 
of pinch singularities is a consequence of the consistent inclusion of finite-time effects and 
the proper consideration of the time of observation. After introducing a physically mean- 
ingful definition of particle number densities, we derive master time evolution equations 
for statistical distribution functions, which are valid to all orders in perturbation theory. 
The resulting equations do not rely upon a gradient expansion of Wigner transforms or 
involve any separation of time scales. To illustrate the key features of our formalism, we 
study out-of-equilibrium decay dynamics of unstable particles in a simple scalar model. 
In particular, we show how finite-time effects remove the pinch singularities and lead to 
violation of energy conservation at early times, giving rise to otherwise kinematically for- 
bidden processes. The non-Markovian nature of the memory effects as predicted in our 
formalism is explicitly demonstrated. 
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1 Introduction 



As modern particle physics continues to advance at both the energy and intensity fron- 
tiers, we are increasingly concerned with transport phenomena in dense systems of ultra- 
relativistic particles, the so-called density frontier. One such system is the deconfined phase 
of Quantum Chromodynamics, known as the quark-gluon plasma [1], whose existence has 
been inferred from the observation of jet quenching in Pb-Pb collisions by the ATLAS [2], 
CMS [3] and ALICE [4] experiments at the CERN Large Hadron Collider (LHC). In ad- 
dition to laboratory experiments, an understanding of such ultra-relativistic many-body 
dynamics is of interest in theoretical astro-particle physics and cosmology. Predictions 
about the evolution of the early Universe rely upon our understanding of the dynamics of 
states at the end and shortly after the phase of cosmological inflation. 

The Wilkinson Microwave Anisotropy Probe (WMAP) [5, 6] measured a baryon-to- 
photon ratio at the present epoch of rj = ns/n-y = 6.1161q 249 x lO"^*^, where is the 
number density of photons and ns = ni, — is the difference in number densities of 
baryons and anti-baryons. This observed Baryon Asymmetry of the Universe (BAU) is also 
consistent with the predictions of Big Bang Nucleosynthesis (BBN) [7]. The generation of 
the BAU requires the presence of out-of-equilibrium processes and the violation of baryon 
number (B), charge (C) and charge-parity {CP), according to the Sakharov conditions [8]. 
One such set of processes is prescribed by the scenarios of baryogenesis via leptogenesis [9- 
11], in which an initial excess in lepton number (L), provided by the decay of heavy right- 
handed Majorana neutrinos, is converted to a baryon number excess through the B + L- 
violating sphaleron interactions. The description of such phenomena require a consistent 
approach to the non-equilibrium dynamics of particle number densities. Further examples 
to which non-equilibrium approaches are relevant include, for instance, the phenomena of 
reheating and preheating [12-15] and the generation of dark matter relic densities [16]. 

The classical evolution of particle number densities in the early Universe is described by 
Boltzmann transport equations; see for instance [17-24]. A semi-classical improvement to 
these equations may be achieved by substituting the classical Boltzmann distributions with 
quantum-statistical Bose-Einstein and Fermi-Dirac distribution functions for bosons and 
fermions, respectively. However, such improved approaches take into account finite-width 
and off-shell effects by means of effective field-theoretic methods. Hence, a complete and 
systematic field-theoretic description of quantum transport phenomena would be desirable. 

The first framework for calculating Ensemble Expectation Values (EEVs) of field oper- 
ators was provided by Matsubara [25]. This so-called Imaginary Time Formalism (ITF) of 
thermal field theory is derived by interpreting the canonical density operator as an evolution 
operator in negative imaginary time. Real-time Green's functions may then be obtained 
by subtle analytic continuation. Nevertheless, the applicability of the ITF remains limited 
to the description of processes occurring in thermodynamic equilibrium. 

The calculation of EEVs of operators in non-static systems is performed using the 
so-called Real Time Formalism (RTF); see for example [26, 27]. In particular, for non- 
equilibrium systems, one uses the Closed Time Path (CTP) formalism, or the in-in for- 
malism, due to Schwinger and Keldysh [28, 29]. The correspondence of these results with 
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those obtained by the ITF in the equihbrium hmit are discussed extensively in the hter- 
ature [30-36]. A non-perturbative loopwise expansion of the CTP generating functional 
is then provided by the Cornwall-Jackiw-Tomboulis (CJT) effective action [37, 38], which 
was subsequently applied to the CTP formalism by Calzetta and Hu [39, 40]. 

Recently, the computation of the out-of-equilibrium evolution of particle number den- 
sities has received much attention, where several authors put forward quantum-corrected or 
quantum Boltzmann equations [41-50]. These approaches rely upon the Wigner transfor- 
mation and gradient expansion [51] of a system of Kadanoff-Baym (KB) [52, 53] equations, 
originally applied in the non-relativistic regime [54, 55]. Often the truncation of the gradi- 
ent expansion is accompanied by quasi-particle ansaetze for the forms of the propagators. 
Similar approaches have recently been applied to the quark-gluon plasma [56]. Dynamical 
equations may also be derived using functional renormalization techniques [57]. 

In this article, we present a new approach to non-equilibrium thermal field theory. 
Our approach is based upon a diagrammatic perturbation series, constructed from non- 
homogeneous free propagators and time-dependent vertices, which explicitly break space- 
time translational invariance. In particular, we show how the systematic inclusion of finite- 
time effects and the proper consideration of the time of observation render our perturbative 
expansion free of the pinch singularities, thereby enabling a consistent treatment of non- 
equilibrium dynamics. Unlike other methods, our approach does not require the use of 
quasi-particle approximation or other effective resummations of finite- width effects. 

A key element of our formalism has been to define physically meaningful particle num- 
ber densities in terms of off-shell Green's functions. This definition is unambiguous, as it 
can be closely linked to Noether's charge, associated with a partially conserved current. 
Subsequently, we derive master time evolution equations for the statistical distribution 
functions related to particle number densities. These time evolution equations do not rely 
on the truncation of a gradient expansion of the so-called Wigner transforms, neither do 
they involve separation of various time scales. Instead, they are built of non-homogeneous 
free propagators, with modified time-dependent Feynman rules, which enable us to analyze 
the pertinent kinematics fully. Our analysis shows that the systematic inclusion of finite- 
time effects leads to the microscopic violation of energy conservation at early times. Aside 
from preventing the appearance of pinch singularities, the effect of a finite time interval of 
evolution leads to contributions from processes that would otherwise be kinematically disal- 
lowed on grounds of energy conservation. Applying this formalism to a simple scalar model 
with unstable particles, we show that these evanescent processes contribute significantly 
to the rapidly-oscillating transient evolution of these systems, inducing late-time memory 
effects. We find that the spectral evolution of two-point correlation functions exhibits an 
oscillating pattern with time-varying frequencies. Such an oscillating pattern signifies a 
non-Markovian evolution of memory effects, which is a distinctive feature governing truly 
out-of-thermal-equilibrium dynamical systems. 

The outline of the paper is as follows. In Section 2, we review the canonical quan- 
tization of a scalar field theory, placing particular emphasis on the inclusion of a finite 
time of coincidence for the three equivalent pictures of quantum mechanics, namely the 
Schrodinger, Heisenberg and Dirac (interaction) pictures. In Section 3, we introduce the 
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Table 1: Glossary for clarifying polysemous notation. 

OTP formalism, limiting ourselves initially to consider its application to quantum field 
theory at zero temperature and density. This is followed by a discussion of the constraints 
upon the form of the OTP propagator. With these prerequisites reviewed, we proceed in 
Section 4 to discuss the generalization of the OTP formalism to finite temperature and den- 
sity in the presence of both spatially and temporally inhomogeneous backgrounds. In the 
same section, we derive the most general form of the non-homogeneous free propagators for 
a scalar field theory. The thermodynamic equilibrium limit is outlined in Section 5. In Sec- 
tion 6, we demonstrate the absence of pinch singularities in the perturbation series arising in 
our approach. Subsequently, in Section 7, we define the concept of particle number density 
and relate this to a perturbative loopwise expansion of the resummed OTP propagators. In 
Section 8, we derive new master time evolution equations for statistical distribution func- 
tions, which go over to classical Boltzmann equations in the appropriate limits. Section 9 
studies the thermalization of unstable particles in a simple scalar model, where particu- 
lar emphasis is laid on the early-time behaviour and the impact of the finite-time effects. 
Finally, our conclusions and possible future directions are presented in Section 10. 

For clarity, a glossary that might be useful to the reader to clarify polysemous no- 
tation is given in Table 1. Appendix A provides a summary of all important propagator 
definitions, along with their basic relations and properties. In Appendix B, we describe 
the correspondence between the RTF and ITF in the equilibrium limit at the one-loop 
level for a real scalar field theory with a cubic self-interaction. The generalization of our 
approach to complex scalar fields is outlined in Appendix C. Appendix D contains a series 
expansion of the most general non-homogeneous Gaussian-like density operator. In Ap- 
pendix E, we summarize the derivation of the so-called Kadanoff-Baym equations. Lastly, 
in Appendix F, we describe key technical details involved in the calculation of loop integrals 
with non-homogeneous free propagators. 
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2 Canonical Quantization 



In this section, we review the basic results obtained within the canonical quantization 
formalism for a massive scalar field theory. This discussion will serve as a precursor for our 
generalization to finite temperature and density, which follows in subsequent sections. 

As a starting point, we consider a simple self-interacting theory of a real scalar field <I>(x) 
with mass M, which is described by the Lagrangian 

C{x) = ldf,^{x)d^'<^{x) - iM^'^^x) - - ^A$^(x) , (2.1) 

where g and A are dimensionful and dimensionless couplings, respectively. Throughout this 
article, we use the short-hand notation: x = x^ = (x'^,x), for the four-dimensional space- 
time arguments of field operators, and adopt the signature (+,—,—,—) for the Minkowski 
space-time metric r/^j^. 

It proves convenient to start our canonical quantization approach in the Schrodinger 
picture, where the state vectors are time- dependent, whilst basis vectors and operators 
are, in the absence of external sources, time-independent. Hence, the time-independent 
Schrodinger-picture field operator, denoted by a subscript S, may be written in the familiar 
plane-wave decomposition 

$s(x;ti) = jdli^ (as(p;ti)e*P-" + 4(p; t;)e-^P-" ) , (2.2) 

where we have introduced the short-hand notation: 

for the Lorentz-invariant phase-space (LIPS). In (2.3), -E'(p) = \/p^ + is the energy 
of the single-particle mode with three- momentum p and 0{pq) is the generalized unit step 
function, defined by the Fourier representation 



6{po) = i 



2tt i + ie 



1 , > ; 

\, Po = 0; (2.4) 
, PO < , 



where e = 0^. It is essential to stress here that we define the Schrodinger-, Heisenberg- and 
Interaction (Dirac)-pictures to be coincident at the finite microscopic boundary time tj, 
such that 

$s(x;ti) = $H(ti,x;ti) = $i(ii,x;ti), (2.5) 

where implicit dependence upon the boundary time ti is marked by separation from explicit 
arguments by a semi-colon. 

The time-independent Schrodinger-picture operators ag(p;tj) and as(p;tj) are the 
usual creation and annihilation operators, which act on the stationary vacuum |0), re- 
spectively creating and destroying time-independent single-particle momentum eigenstates. 
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Their defining properties are: 

4(p;ti)|0) = \p;ii) , (2.6a) 
as{p;ii)\p';U) = {27rf2E{p)6^^\p - p')\0) , (2.6b) 
«s(p;ii)|0) = 0. (2.6c) 

Note that the momentum eigenstates |p;tj) respect the orthonormahty condition 

{p';ii\p;ii) = {27rf2Eip)S^'\p - p') . (2.7) 

We then define the time- dependent interaction-picture field operator ^i(x;ti) via 

^l{x;ii) = e^^s(^o-t.)$g(x;t.)e-*^s(^o-*«) , (2.8) 

where i?g is the free-part of the Hamiltonian in the Schrodinger picture. Using the com- 
mutators 

[Hi asip;ii)] = -Eip)asip;ii) , (2.9a) 
[Hi alip;ii)] = +E{p)alip;ii) , (2.9b) 

the interaction-picture field operator may be written 

$i(x;t~) = y"dnp(as(p;t;)e-*^(P)("o-*>)e^P-" + 4(p;ti)e^^(P)("«-*')e-^P-'') , (2.10) 

or equivalently, in terms of interaction-picture operators only, 

<^iix;ti) = ^dHp (ai(p,0;ti)e-^^(P)^Oe^P-'' + aJ(p,0;t;)e^^(P)^°e~^P-") . (2.11) 

Notice that in (2.11) the time- dependent interaction-picture creation and annihilation op- 
erators, a\{p,i;ii) and ai{p,i;ii), are evaluated at the microscopic time i = 0, after 
employing a relation analogous to (2.8). We may write the four-dimensional Fourier trans- 
form of the interaction-picture field operator as 

Mp;ii) = Jd''xe'P-'^<^i{x;U) = 27r5(/ - M^) ( 0(po)ai(p, 0; t^) + 0(-po)4(-p, 0; t^) ) . 

(2.12) 

In the limit where the interactions are switched off adiabatically as ti — )• — oo, one 
may define the asymptotic in creation and annihilation operators via 

aJP(p) = lim a\^\p,0;ii) = Z-^/\lim a^^\p;ii)e+^-^'^^P^'\ (2.13) 

ii — > — oo ii — > — oo 

where Z = 1 + 0{h) is the wavefunction renormalization. Evidently, keeping track 
of the finite boundary time ti plays an important role in ensuring that our forthcoming 
generalization to perturbative thermal field theory remains consistent with asymptotic field 
theory in the limit ti — t- — oo. Hereafter, we will omit the subscript I on interaction-picture 
operators and suppress the implicit dependence on the boundary time ti, except where it 
is necessary to do otherwise for clarity. 
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We start our quantization procedure by defining the commutator of interaction-picture 
fields 

mx), <^{y)] = iA\x,y;M^) , (2.14) 

where iA^{x,y; M'^) is the free Pauh-Jordan propagator. Herein, we denote free propa- 
gators by a superscript 0. The condition of micro-causahty requires that the interaction- 
picture fields commute for space-like separations (x — y)^ < 0. This restricts iA^{x, y; M^) 
to be invariant under spatial translations, having the Poincare-invariant form 

iA\x,y;M^) = J dUp g-^^^P^^^" " ^''^e^P'^'' " - (x ^ y)) . (2.15) 

Observe that iA^{x, y; M^) represents the difference of two counter-propagating packets of 
plane waves and vanishes for (x — y)^ < 0. 

It proves useful for our forthcoming analysis to introduce the double Fourier transform 

iA^{p,p';M^) = JJd^xdS e'P-^'e-'P'-y iA\x,y;M^) (2.16a) 
= 2Tre{po)dip^ - M^){2tt)^6^^'>{p - p') , (2.16b) 

where £{po) = 0{pq) — 0(— po) is the generalized signum function. Note that we have 
defined the Fourier transforms such that four-momentum p flows away from the point x 
and four-momentum p' flows towards the point y. 

From (2.14) and (2.15), we may derive the equal-time commutation relations 

[$(t,x), Hi,y)] = 0, (2.17a) 
[^(t,x), cD(t,y)] = -i5(3)(x - y) , (2.17b) 
[7r(t,x), 7r(t,y)] = 0, (2.17c) 

where 7r(x) = dxQ^{x) is the conjugate-momentum operator. In order to satisfy the canon- 
ical quantization relations (2.17), the creation and annihilation operators must respect the 
commutation relation 

[a(p,t), a\p',i')] = (27r)32ii;(p)<5(3)(p _ p')e-i^(p)(t - , (2.I8) 

with all other commutators vanishing. Here, we emphasize the appearance of an overall 
phase e-^^(p)(*-*') in (2.18) for i / i' [cf. Section 4.1]. 

The vacuum expectation value of the commutator of Heisenberg-picture field operators 
may be expressed as a superposition of interaction-picture field commutators by means of 
the Kallen-Lehmann spectral representation [58, 59]: 

(0|[$h(x), $H(y)]|0) = iA{x,y) = / dsa(s)iAO(x,y;s) , (2.19) 

Jo 

where iA^{x,y] s) is the free Pauli- Jordan propagator in (2.15) with replaced by s 
and iA{x,y) is the dressed or reswmmed propagator. The positive spectral density (j{s) 



iAO(x,y;M2)|^o^^o^- = 
a,.„iA0(x,y;M2)|^„^^„^- = 

5xo9yo«A°(x,y;M2)|^o= o = t = 
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contains all information about the spectrum of single- and multi-particle states produced 
by the Heisenberg-picture field operators $h- Note that for a homogeneous and stationary 
vacuum |0), cr(s) is independent of the space-time coordinates and the resummed Pauli- 
Jordan propagator maintains its translational invariance. If a{s) is normalized, such that 

o 

ds a{s) = 1 , (2.20) 

the equal-time commutation relations of Heisenberg-picture operators maintain exactly 
the form in (2.17). In this case, the spectral function cannot depend upon any fluctu- 
ations in the background. Clearly, for non-trivial 'vacua', or thermal backgrounds, the 
spectral density becomes a function also of the coordinates. The spectral representation 
of the resummed propagators may then depend non-trivially on space-time coordinates, 
i.e. a = a{s;x,y); see for instance [41]. In this case, the convenient factorization of the 
Kallen-Lehmann representation breaks down and the equal-time commutation relations of 
Heisenberg-picture field operators do not maintain their canonical form. 
The retarded and advanced causal propagators are defined as 

iAn{x,y) = e{xo - yo)iA{x,y) , iAA{x,y) = -9{yo - xo)iA{x,y) . (2.21) 

Using the Fourier representation of the unit step function in (2.4), we introduce a convenient 
representation of these causal propagators in terms of the convolution 

'dko iA{po -_ko,p'Q - ko;p,p') 
to 

The absolutely-ordered Wightman propagators are defined as 



^Ak(A) {P,P) = ^| — T— TZy- • (2-22) 



iA>(a;,y) = ($(x)$(y)) , iA<{x,y) = {^{y)^{x)) . (2.23) 

We note that the two-point correlation functions A(x, y), A>_<(x, y) and A'^^x{x, y) satisfy 
the causality relation: 

A{x,y) = A>{x,y) - A<(x,y) = AR(x,y) - AA{x,y) . (2.24) 

Our next step is to define the non-causal Hadamard propagator, which is the vacuum 
expectation value of the field anti-commutator 

iAi{x,y) = ({^(x), Hy)}). (2.25) 

Correspondingly, the time-ordered Feynman and anti-time-ordered Dyson propagators are 
given by 

iAF(x,y) = (T[$(x)$(2/)]) , iAD(x,y) = (t[$(x)$(y)]) , (2.26) 

where T and T are the time- and anti-time-ordering operators, respectively. Explicitly, 
Af(x, y) and Ad(x, y) may be written in terms of the absolutely-ordered Wightman prop- 
agators A>(x,y) and A^{x,y) as 

AF{x,y) = e{xo - yo)A>(x,y) + 9{yo - xo)A<(x,y) , (2.27a) 
Auix,y) = 9{xo - yo)A<(x,y) + 9{yo - xo)A>(x,y) . (2.27b) 
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The propagators Ai{x,y), A>^<(a;,y) and Apj:){x,y) obey the unitarity relations: 

Ai(x,y) = AF(x,y) + A^{x,y) = A>(x,y) + A<(x,y) = 2ilm Af(x, y) . (2.28) 
Finally, for completeness, we define the principal-part propagator 

Ap(x,y) = ^(AR(x,y) + AA(x,y)) = ReAF(x,y). (2.29) 
Here, we should bear in mind that 

Re(lm)AF(x,2/) / e^f-e^^'-^ Re(lm) Af(p,p') , (2-30) 

implying that 

Ap(p,p') / ReAF(p,p') , (2.31) 

unless Af(p,p') = Af(— — p')) which is not generally true in non-equilibrium thermal 
field theory [cf. (A.4f)]. 

The definitions and the relations discussed above are valid for both free and resummed 
propagators. In Appendix A, we list the properties of these propagators in coordinate, 
momentum and Wigner (see Section 4.2) representations, as well as a number of useful 
identities, which we will refer to throughout this article. More detailed discussion related 
to these propagators and their contour- integral representations may be found in [60]. In 
Appendix C, these considerations and the analysis of the following sections are generalized 
to the complex scalar field. 

3 The CTP Formalism 

In this section, we review the Closed Time Path (CTP) formalism, or the so-called in-in 
formalism, due to Schwinger and Keldysh [28, 29]. As an illuminating exercise, we consider 
the CTP formalism in the context of zero-temperature quantum field theory and derive 
the associated 2x2 matrix propagator, obeying basic field-theoretic constraints, such as 
CPT invariance, Hermiticity, causality and unitarity. Finally, we discuss the properties of 
the resummed propagator in the CTP formalism. 

In the calculation of scattering-matrix elements, we are concerned with the transition 
amplitude between in and out asymptotic states, where single-particle states are defined in 
the infinitely-distant past and future. On the other hand, in quantum statistical mechanics, 
we are interested in the calculation of Ensemble Expectation Values (EEVs) of operators 
at a fixed given time t. Specifically, the evaluation of EEVs of operators requires a field- 
theoretic approach that allows us to determine the transition amplitude between states 
evolved to the same time. This approach is the Schwinger-Keldysh CTP formalism, which 
we describe in detail below. 

For illustration, let us consider the following observable O in the Schrodinger picture: 

0(s)(t/;ti) = /[d$] g($(z),^7;^,|$s(x;^0^s(y;^^)|^(z),^/;^^)s , (3.i) 
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where [d$] represents the functional integral over all field configurations '&(z). In (3.1), 
|$(z), fj^; ij)g is a time- dependent eigenstate of the time-independent Schrodinger-picture 
field operator $s(x;ti) with eigenvalue $(x), where the implicit dependence upon the 
boundary time ti has been restored. 

We should remark here that there are seven independent space-time coordinates in- 
volved in the observable (3.1). These are the six spatial coordinates, x and y, plus the 
microscopic time tf. In addition, there is one implicit coordinate: the boundary time ti. 
As we will see, exactly seven independent coordinates are required to construct physical 
observables that are compatible with Heisenberg's uncertainty principle. We choose the 
seven independent coordinates to be 

t = if - U, X = ^ (x + y) , p , (3.2) 

where t and X are the macroscopic time and central space coordinates and p is the Fourier- 
conjugate variable to the relative spatial coordinate R = x — y. 

In the interaction picture, the same observable O in (3.1) is given by 

Oii){if;ii) = y[d$] j($(z),t7;ti|$i(t7,x;ti)^>i(t7,y;t~)|^>(z),t/;ti)j ; (3.3) 

and, in the Heisenberg picture, by 

0(H){if;U) = y"[d$] jj($(z);ti|$H(t/,x;ti)$H(t/,y;ti)|$(z);ti)H • (3.4) 

Notice that the prediction for the observable O does not depend on which picture we are 
using, i.e. 0(g)(tj; tj) = 0^i'^{tf\ii) = 0(Y{-^{if;ii). This picture-independence of O is only 
possible because the time-dependent vectors and operators in O are evaluated individually 
at equal times. Otherwise, any potential observable, built out of time-dependent vectors and 
operators that are evaluated at different times, would be picture- dependent and therefore 
unphysical. Moreover, the prediction of the observable O should be invariant under time 
translations, i.e. 

0{if;ii) = O{t;0) = 0{t), (3.5) 

where t = tf — ti is the macroscopic time, as we will see below. Herein and throughout the 
remainder of this article, the time agruments of time translationally-invariant quantities 
are written in terms of the macroscopic time t only. 

3.1 The CTP Contour 

In order to evaluate equal-time observables of the form in (3.4), we first introduce the in 
vacuum state |Oin;ti), which is at time ti a time-independent eigenstate of the Heisenberg 
field operator ^u{x; ti) with zero eigenvalue; see [39, 40]. We then need a means of driving 
the amplitude {Om] ii\Om, U) , which can be achieved by the appropriate introduction of 
external sources. 

The procedure may be outlined with the aid of Figure 1 as follows. We imagine 
evolving our in state at time ti forwards in time under the influence of a source J+{x) to 
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Imt 



z(0) = u 



Ret 



z(l) = U -ie 



z(l/2) = if ~ te/2 



macroscopic time t = Re z(u) — ti 



initial conditions: 
macroscopic time t = 



observation: 
macroscopic time i — tf — ti 



Figure 1: The closed-time path, C = C+ U C_, running first along C+ from ti to tf — ie/2 
and then returning along C_ from — to ti — ie. The relationship between microscopic 
and macroscopic times is indicated by a dashed black arrow. 

some out state at time t/ in the future, which will be a superposition over all possible future 
states. We then evolve this superposition of states backwards again under the influence of 
a second source J-{x), returning to the same initial time ti and the original in state. The 
sources J±{x) are assumed to vanish adiabatically at the boundaries of the interval [ii, if]. 
We may interpret the path of this evolution as defining a closed contour C = C+ U C_ in 
the complex-time plane (t-plane, i G C), which is the union of two anti-parallel branches: 
C+, running from ii to if — ie/2; and C_, running from if — ie/2 back to ii — ie. We 
refer to C+ and C_ as the positive time-ordered and negative anti-time-ordered branches, 
respectively. As depicted in Figure 1, the small imaginary part e = 0^ has been added 
to allow us to distinguish the two, essentially coincident, branches. We parametrize the 
distance along the contour, starting from ti, by the real variable u € [0, 1], which increases 
monotonically along C. We may then define the contour by a path z{u) = t{u) — if{u) G t, 
where t(0) = = ii and i(l/2) = if. Thus, the complex CTP contour z{u) may be 
written down explicitly as 

z{u) = e{\-u) ii + 2u[if-ii) +e{u-\) ii + 2 [l - u) {if - ii) - ieu , (3.6) 

with 6^(0) = 1/2, as given in (2.4). 

In the presence of external sources J±{x), the eigenstate |<I>(x), i; j of the Heisen- 
berg field operator $H(i)X; ii) becomes explicitly time-dependent, satisfying the eigenvalue 
equation 

$H(t,x;ti)|$(x),t;t~)H,j = $(x)|^>(x),t;ti)H,j • (3-7) 
Also, |<I>(x), i; J forms a complete orthonormal basis, respecting the orthonormality 
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condition 

J[d^] = loo- (3.8) 

We may then write the CTP generating functional Z [J± , t] as 

Z[J±,t] = j[d<^] ^_(Oin,ti;ti|$(x),t/;ti)y X ^ j^{<^{x),if;ii\Oin,ii;ii)j^ 
= /[d$] {Oin,ti;ti\T(e ^ ^ M|«>(x),t/;ti) 

J -7— V / H, J_ 

X (a>(x),t,;t,|T(e^4''^'^*«(^)^-(^))|0i.,t,;t,) , (3.9) 

H, \ / J+ 

where the integrations run from tj to and the 'latest' time (with u = 1) appears 
furthest to the left. 

In order to preserve the correspondence with the ordinary S-matrix theory in the 
asymptotic limit tj — t- —do, we take tf = —ti. With this identification, the CTP gen- 
erating functional Z\J±,t] becomes manifestly CPT invariant. This enables one to easily 
verify that microscopic CPT invariance continues to hold, even when time translational 
invariance is broken by thermal backgrounds, as we will see in Section 4. Given that ti is 
the microscopic time at which the three pictures of quantum mechanics coincide and the 
interactions are switched on, it is also the point at which the boundary conditions may be 
specified fully and instantaneously in terms of on-shell free particle states. The microscopic 
time ti is therefore the natural origin for a macroscopic time t = tf — ti = 2tf. This 
macroscopic time is also the total interval of microscopic time over which the system has 
evolved. This fact is illustrated graphically in Figure 1. 

We denote by $±(x) = ^{x^ G C±,x) fields with the microscopic time variable 
confined to the positive and negative branches of the contour, respectively. Following [39, 
40], we define the doublets 

$'^(x) = (cl>+(x), cl>_(x)) , <^a{x) = Vab^Hx) = (<I>+(x), -<I>_(x)) , (3.10) 

where the CTP indices a, 6 = 1, 2 and r]ab = diag (1, —1) is an SO (1, 1) 'metric' Inserting 
into (3.9) complete sets of eigenstates of the Heisenberg field operator at intermediate times, 
we may derive a path-integral representation of the CTP generating functional: 



Z[Ja,t] = AA/[d$"] exp [z + / d^x Ja(x)$"(x) 



(3.11) 



where M is some normalization and fit ~ [—t/2, t/2] X is the Minkowski space-time 
volume bounded by the hypersurfaces x^ = ±t/2. In (3.11), the action is 

= [ d^x \^r,abd^^''{x)d^^\x) - \{M^r]ab - ielab)^^{x)^\x) 

- ^gvabc'^''{x)'^''{x)^'{x) - j,Xl^abcd^''{x)^\x)<^%x)^''ix)] , (3.12) 

where rjabc- = +1 for a = 6 = c = • • • = 1, r]abc- = — 1 for a = 6 = c = • • • = 2, and 
Vabc— = otherwise. In (3.12), the e = 0"^ gives the usual Feynman prescription, ensuring 
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convergence of the CTP path integral. We note that the damping term is proportional to 
the identity matrix 1^^ and not to the 'metric' ijab- This prescription requires the addition 
of a contour-dependent damping term, proportional to — u), which has the same sign 
on both the positive and negative branches of the contour, C+ and C_, respectively. 

In order to define properly a path-ordering operator Tq, we introduce the contour- 
dependent step function 



where 
function 



and 



y 



Ocix'^-yn = e{ux-uy) , (3.13) 
z{uy). By analogy, we introduce a contour-dependent delta 

6{Ux - Uy) _ 5{Ux - Uy) _ 6{Ux - Uy) 



I — I 

i du i 



2\tf - ti\ 2t 
As a consequence, a path-ordered propagator A.c{x,y) may be defined as follows: 



iAc(x,y) = (Tc [$(x)$(y); 



(3.14) 



(3.15) 



For and on the positive branch C+, the path-ordering Tc is equivalent to the 
standard time-ordering T and we obtain the time-ordered Feynman propagator iAF(a;,y). 
On the other hand, for x^ and y^ on the negative branch C_, the path-ordering Tc is 
equivalent to anti-time-ordering T and we obtain the anti-time-ordered Dyson propaga- 
tor iA-Q^x, y). For x^ on C+ and y^ on C_, x^ is always 'earlier' than and we obtain the 
absolutely-ordered negative-frequency Wightman propagator iA<(x,y). Conversely, for 
on C+ and x^ on C_, we obtain the positive-frequency Wightman propagator iA>(x, y). In 
the SO (1, 1) notation, we write the CTP propagator as the 2x2 matrix 



iA'^''(x,y) ^ {Tc[^''{x)^\y)]) = i 



AF(x,y) A<(x,7/) 
^>{x,y) ^Y){x,y) 



(3.16) 



In this notation, the CTP indices a, h are raised and lowered by contraction with the 
'metric' r/ab, e.g. 



iAab{x,y) = r]aciA'"''ix,y)r]db 



AF{x,y) -A<(x,y) 
-A>(x,y) AD(x,y) 



(3.17) 



Notice the difference in sign of the off-diagonal elements in (3.17) compared with (3.16). An 
alternative definition iA"'^{x, y) of the CTP propagator uses the so-called Keldysh basis [35] 
and is obtained by means of an orthogonal transformation: 



A'^^(x,y) = 0\0\, 



. cd 



{x,y) 



AA(x,y) 
^■R.{x,y) Ai(x,y) 



O 



ah 



1 

V2 



1 1 
1 -1 



. (3.18) 



In Section 4, we will generalize these results to macroscopic ensembles by incorpo- 
rating background effects in terms of physical sources. In this case, the surface integral 
§dnt '^'^A' ^°'{x)d^^^^{x) may not in general vanish on the boundary hypersurface dO^t of the 
volume However, by requiring the '-|-'- and '—'-type fields to satisfy 
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we can ensure that surface terms cancel between the positive and negative branches, 
and C- , respectively. In this case, the free part of the action may be rewritten as 



50[<i>^^] = // d^xd^y if»(x)A°r'(x,y)$''(y) 



where 



A' 



0,-1 
ab 



{x,y) = 5^^\x - y)\-{Ul + M^)r^ab + ielab 



(3.20) 



(3.21) 



is the free inverse CTP propagator and = is the d'Alembertian operator. 

Note that the variational principle remains well-defined irrespectively of (3.19), since we 
are always free to choose the variation of the field 6^°'{x) to vanish for on dQ.t. 

We may complete the square in the exponent of the CTP generating functional Z[Ja-, t\ 
in (3.11) by making the following shift in the field: 



$"'(x) - f d 



d^y A^'''\x,y)My) , 



(3.22) 



where zA'^'"''(x, y) is the free CTP propagator. We may then rewrite Z[Ja,t] in the form 
Z[Ja,t] = ^0[0,t]exp / d^xC'^^'i-^- 

d^xd^yJaix)A'^''\x,y)My) 



X exp 



(3.23) 



where >C'°* is the interaction part of the Lagrangian and ^"^[0,^] is the free part of the 
generating functional 



ZO[0,t] = Mj [d$"] e^^ , (3.24) 

with the free action S^[^"',t\ given by (3.20). We may then express the resummed CTP 
propagator iA"-^{x,y,if;ii) as follows: 



zA^\x,y,tj-U) = ^^ir^^jYT^^jYr-^^iJa,t] 



(3.25) 



Ja=0 



where 2^[0, t] is the generating functional with the external sources Ja set to zero. The 
functional derivatives satisfy 

6 



5Ja{x) 



d^yj\y) = ^-'d^'^Xx - y), 



(3.26) 



with x^, y^ G ilf We will see in Section 4.3 that the resummed CTP propagator 
A°'''{x,y,if;ii) is not in general time translationally invariant. 

In the absence of interactions, eigenstates of the free Hamiltonian will propagate unin- 
terrupted from times infinitely distant in the past to times infinitely far in the future. As 
such, we may extend the limits of integration in the free part of the action to positive and 
negative infinity, since 







(3.27) 
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i.e. the sources Ja{x) vanish for ^ [—t/2, t/2]. The free CTP propagator A^'"'''{x,y) is 
then obtained by inverting (3.21) subject to the inverse relation 



dSA°-i(x,z)AO'^-(z,y) = rj,-6^^\x - y) 



(3.28) 



where the domain of integration over is extended to infinity. We expect to recover 
the famihar propagators of the in- out formahsm of asymptotic field theory, which occur 
in S-matrix elements and in the reduction formalism due to Lehmann, Symanzik and 
Zimmermann [61]. The propagators will also satisfy unitarity cutting rules [62-65], thereby 
maintaining perturbative unitarity of the theory. Specifically, the free Feynman (Dyson) 
propagators iAp^j^-^{x,y) satisfy the inhomogeneous Klein-Gordon equation 



- {Ol + M2)iA0(D)(x,y) = {-)i5^'\x - y) ; 
and the free Wightman propagators iA^(x,2/) satisfy the homogeneous equation 

-{Ul + M2)iA0(<)(x,y) = 0. 



(3.29) 



(3.30) 



In the double momentum representation, the free part of the action (3.20) may be 
written as 

d^p dV l^a/ 



-cI>-(p)A°-^(p,pO<f''(-y) , 



27r 27r 



where 



m + ie 




{p^ — rm?' — ie) 



(2^)^5(4) (p - p') 



(3.31) 



(3.32) 



is the double momentum representation of the free inverse CTP propagator, satisfying the 
inverse relation 



d^g 



A'l;-\p,q)A^^''{q,p') = r/,^(27r)V^)(p - p') 



(3.33) 



(27r)4 "f- 

Given that the free inverse CTP propagator is proportional to a four- dimensional delta 
function of the two momenta, it may be written more conveniently in the single Fourier 
representation as 



p^ — m? + ie 







—{p^ — m? — ie) 

Hence, for translationally invariant backgrounds, we may recast (3.32) in the form 



A 



0,-1, 
ab 



where 



,0,ab 



(3.34) 

(3.35) 
(3.36) 



is the single-momentum representation of the free CTP propagator. 
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3.2 The Free CTP Propagator 

We proceed now to make the following ansatz for the most general translationally invariant 
form of the free CTP propagator, without evaluating the correlation functions directly: 

_ - M2 + ie)~^ + Ci(p)(5(p2 - M2) C3(p)5(p2 - M2) 

C2(p)(5(p2 - - M2 - ie)~^ + C4(p)(5(p2 - M^) 

(3.37) 

The Ci{p) = 0{po)ci{p) + ^(— Po)c'j(p) are as yet undetermined analytic functions of 
the four-momentum p, which may in general be complex. The diagonal elements are the 
Fourier transforms of the most general translationally invariant solutions to the inhomoge- 
neous Klein-Gordon equation (3.29), whilst the off-diagonal elements are the most general 
translationally invariant solutions to the homogeneous Klein-Gordon equation (3.30). 

The remaining freedom in the matrix elements of A^'°'''{p) is determined by the fol- 
lowing field-theoretic requirements: 

(i) CPT Invar iance. Since the action should be invariant under CPT, the real scalar 
field $ should be even under CPT. From the definitions of the propagators in (A.l), we 
obtain the CPT relations in (A. 3). Consequently, the momentum representation of these 
relations in (A. 4) imply that 

ci(4)(p) = ci(4)(-p) , C2{p) = h{-p) ■ (3.38) 

(ii) Hermiticity. The Hermiticity properties of the correlation functions defined in (A.l) 
give rise to the Hermiticity relations outlined in (A. 4). These imply that 

C4(p) = -c\{p) , Up) = -cU-p) ■ (3.39) 

In conjunction with (3.38), we conclude that C2{p) and C3{p) must be purely imaginary- 
valued functions of the four-momentum p. 

(iii) Causality. The free Pauli- Jordan propagator A^{x,y) is proportional only to the real 
part of the free Feynman propag ator ReAO(x,y) [cf. (A.6b)]. The addition of an even- 
parity on-shell dispersive part to the Fourier transform of the free Feynman propagator 
Ap(p) will contribute to the free Pauli- Jordan propagator terms that are non- vanishing 
for space- like separations (x — y)^ < 0, thus violating the micro-causality condition 
outlined in Section 2. It follows then that ci{p) and C4(p) are also purely imaginary- valued 
functions. We shall therefore replace the Ci{p) by the real- valued functions fi{p) through 
Ci{p) = —2TTifi{p), where the minus sign and factor of 27r have been included for later 
convenience. The explicit form of the free Pauli- Jordan propagator in (2.16b), along with 
the causality relation (2.24), give rise to the constraint 



f2ip) - hip) = e{po) 



(3.40) 



(iv) Unitarity. Finally, the unitarity relations in (2.28) require that 



Mp) + Mp) = 1 + flip) + Mp) ■ 



(3.41) 



Solving the system of the above four constraints for /i,2,3,4(p)) we arrive at the following 
expression for the most general translationally invariant free CTP propagator: 



(p2 _ m2 + 



"2 M2 - 



[P 



M2) 
-1 



27ri/(p)5(p2 



1 1 
1 1 



All elements of A ' (p) contain terms dependent upon the same function 
Kp) = flip) = 0{po)fip) + Oi-po)fi-p) . 



(3.42) 



(3.43) 



These terms correspond to the vacuum expectation of the normal-ordered product of fields 
(:$(a;)$(y):), which vanishes for the trivial vacuum |0). Therefore, we must conclude that 
f{p) also vanishes in this case. We then obtain the free vacuum CTP propagator A^''^^{p), 
which contains the set of propagators familiar from the unitarity cutting rules of absorptive 
part theory [62, 65]: 



A°'"''(p) 



+ ie) 



-2Trie{po)S{p'^ 



-27ri9{-po)S{p'^ 



[P 



M2 



(3.44) 



We may similarly arrive at (3.42) by considering the free CTP propagator in the 
Keldysh representation A^'°'^{p) from (3.18). The constraints outlined above allow us 
to add to the free Hadamard propagator (p) any purely- imaginary even function of p 
proportional to 5{p'^ — M^), that is 



A^^''\p) 





2 „2 



[iPo + ^ej" - P 
- 27ri2f{p)6{p^ 



M2 

2 



[ipO - i6)2 - p2 - M2 

-2TTi6{p^ - M^) 




1 



(3.45) 



We note that there is no such freedom to add terms to the free retarded and advanced 
propagators, Ap^(p) and A\{p), which is a consequence of the micro-causality constraints 
on the form of the free Pauli- Jordan propagator A^{p). Employing the fact that 



O' 






o = 


1 


1 1 






1 




2 


1 1 



(3.46) 



we recover (3.42), which serves as a self-consistency check for the correctness of our ansatz 
for the free CTP propagator. 
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3.3 The Resummed CTP Propagator 

In order to obtain the resummed CTP propagator, we must invert the inverse resummed 
CTP propagator on the restricted domain [—t/2, t/2] subject to the inverse relation 

d'zA-,\x,z,if-k)A'>%z,y,ij-k) = ii:5^'^\x - y) , (3.47) 

t 

for x^, G i}f We shall see in Section 4.1 that this restriction of the time domain implies 
that a closed analytic form for the resummed CTP propagator A"'^{x, y, if;ii) is in general 
not possible, for systems out of thermal equilibrium (see also our discussion in Section 5). 
The double momentum representation of the inverse relation (3.47) takes on the form 

A-,i(p,g,t/;tO(2vr)%^'^(<?-g')A'V,P',t/;iO = " P) , 



{2.r {2.y 

(3.48) 

where we have defined 
Si'\p-p') = Stipo-p'o)d^^Hp-p') = ^|^d4xd4ye^f-e-*f'-^j(^)(x-y). (3.49) 

The restriction of the domain of time integration has led to the introduction of the analytic 
weight function 



t 



27r 



5t{po - p'o) = T^sinc [^-^-^]t , (3.50) 



2 



which has replaced the ordinary energy-conserving delta function. As expected, we have 

lim 6t{po - p'o) = S{po - p'o) , (3.51) 

r — > oo 

SO that the standard description of asymptotic quantum field theory is recovered in the 
limit ti = —t/2 — )• — oo. Moreover, the weight function 5t satisfies the convolution 

/ + 00 
d^o StiPo - 9o)'Jt(go - p'o) = h{Po - p'o) ■ (3.52) 
-oo 

The emergence of the function 5t is a consequence of our requirement that the time evo- 
lution and the mapping between quantum-mechanical pictures (see Section 2) are governed 
by the standard interaction-picture evolution operator 

U{if,ii) = Texp(^iJ^^^diH"'\i;ii)^ . (3.53) 

This evolution is defined for times greater than the boundary time ti, at which point the 
three pictures are coincident. We stress here that 5t is neither a prescription, nor is it an 
a priori regularization of the Dirac delta function. 

As we will see later, the oscillatory behaviour of the sine function in 6t is fundamentally 
important to the dynamical behaviour of the system. Let us therefore convince ourselves 
that these oscillations persist, even if we smear the switching on of the interactions, or 
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Figure 2: Comparison of 5t{po — Pq) (black dotted) and 5a{po — p'o) (blue dashed). The 
arbitrary mass M is included so that axes are dimensionless. 

equivalently, if we impose an adiabatic switching off of the interaction Hamiltonian for mi- 
croscopic times outside the interval [—t/2, t/2]. To this end, we introduce to the interaction 
Hamiltonian i?™* in (3.53) the Gaussian function 

Mi) = exp(^-^) , (3.54) 
such that the evolution operator takes the form 

lilt / 



Uitf,ti) = Texpl^iJ dtAt{t)H'''\t;ti)j . (3.55) 

Clearly, for t ^ t, At{t) — )• 0, whereas for t <^ t, At{t) — t- 1. To account for the effect of 
At{t) in the action, the following replacement needs to be made: 

St{po-Po) ^ ^a{po - p'o) = ^ r'\ie-'^P^-P'^y'At{i) 

= * ^-Upo-p'.?t'\^,^( ^ - 2i(po -P())t \ I ^^J ^ + nPo-Po)t \ 
2^/2^ L V 2^/2 / V 2^/2 J_ 

Due to the error functions of complex arguments in (3.56), the oscillatory behaviour re- 
mains. The analytic behaviour of both the functions 6t{po — Pq) and 5a (po — Po) is shown 
in Figure 2 in which we see that the smooth smearing by the Gaussian function At{t) has 
little effect on the central region of the sine function, as one would expect. 
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4 Non-Homogeneous Backgrounds 



Until now, we have considered the vacuum to be an 'empty' state with all quantum numbers 
zero. In this section, we replace that 'empty' vacuum state with some macroscopic back- 
ground, which may in general be inhomogeneous and incoherent. This non-trivial 'vacuum' 
is described by the density operator p. Following a derivation of the CTP Schwinger-Dyson 
equation, we show that it is not possible to obtain a closed analytic form for the resummed 
CTP propagators in the presence of time-dependent backgrounds. Finally, we general- 
ize the discussions in Section 3.2 to obtain non- homogeneous free propagators in which 
space-time translational invariance is explicitly broken. 

The density operator p is necessarily Hermitian and, for an isolated system, evolves in 
the interaction picture according to the von Neumann or quantum Liouville equation 

= p{i-U)] , (4.1) 

a(t-ti) 

where W^^{t] ti) is the interaction part of the Hamiltonian in the interaction picture, which 
is time-dependent. Notice that the time derivative appearing on the LHS of (4.1) is taken 
with respect to the time translationally invariant quantity t — ti. Developing the usual 
Neumann series, we find that 

p{i;ii) = U-^{i,ti)p{hu)U{i,U) , (4.2) 

where U is the evolution operator in (3.53). Hence, in the absence of external sources and 
given the unitarity of the evolution operator, the partition function Z = Tr p is time- 
independent. On the other hand, the partition function of an open or closed subsystem is 
in general time-dependent due to the presence of external sources. 

We are interested in evaluating time-dependent Ensemble Expectation Values (EEVs) 
of field operators (•)^ at the macroscopic time t, which corresponds to the microscopic 
time If = t/2, where the bra-ket now denotes the weighted expectation 

(•>. ^ ^-^Wj^ ■ f^-^' 

Ticp{tf;ti) 

In this case, EEVs of two-point products of field operators begin with a total of nine 
independent coordinates: the microscopic time of the density operator and the two four- 
dimensional space-time coordinates of the field operators. As discussed in the beginning 
of Section 3 [cf. (3.2)], this number is reduced to the required seven coordinates, i.e. one 
temporal and six spatial, after setting all microscopic times equal to t/ = t/2. Hence, 
physical observables in the interaction picture are, for instance, of the form 

m, x; «.)*(«/, y; y>. ^ Tr(,ii,i.m,^m0f.ym , 

In the presence of a non-trivial background, the out state of Section 3 is replaced by 
the density operator p at the time tf = t/2. Consequently, the starting point of the 
corresponding CTP generating functional of EEVs is 

Z[p,J±,t] = Tr [(Te-*^"t''*"^-("^*"("))pH(t/;iO (Te'^"'''*"'^^^"^*''^"^)] • (4.5) 
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Within the generating functional Z in (4.5), the Heisenberg-picture density operator pn 
has exphcit time dependence, as it is built out of state vectors that depend on time due 
to the presence of the external sources J±. In the absence of such sources, however, the 
state vectors do not evolve in time, so /jh and the partition function Z\p^ -/i = 0, t] = TV p 
in (4.5) become time-independent quantities. 

In contrast to existing approaches, the explicit microscopic time = t/2 of the 
density operator p}i{if;ii) appearing in the CTP generating functional (4.5) is the time 
of observation and not the initial time ti = —t/2, at which the boundary conditions are 
specified. As we shall see in Section 6, this subtle difference will contribute to the absence 
of pinch singularities in the resulting perturbation series. 

4.1 The Schwinger Dyson Equation in the CTP Formalism 

In order to generate a perturbation series of correlation functions in the presence of non- 
homogeneous backgrounds, we must derive the Schwinger-Dyson equation in the CTP 
formalism. Of particular interest is the explicit form of the Feynman-Dyson series for the 
expansion of the resummed CTP propagator. We will show that, in the time-dependent 
case, this series does not collapse to the resummation known from zero-temperature field 
theory. In particular, we find that a closed analytic form for the resummed CTP propagator 
is not attainable in general. 

We proceed by inserting into the generating functional Z in (4.5) complete sets of 
eigenstates of the Heisenberg field operator $h at intermediate times. In this way, we 
obtain the path-integral representation 

2[p,Ja,t] = jid^"] H($_(x),t/;ti|pH(t/;i»)|^+(x),t/;ti)H 



X exp 



^f5[<I>^^] + / d^x Ja(x)cl>'^(x 



(4.6) 



As before, we may extend the limits of integration to infinity in the free part of the action 
and the J-dependent term, due to the fact that the external sources vanish outside the 
time interval [—t/2, t/2]. It is only in the interaction part of the action that the finite 
domain of integration must remain. 

Following [40], we write the kernel ^{^-{x),tf;ti\pii{tf;ti) |$+(x), tj; as an infi- 
nite series of poly-local sources: 



^{^.{^),if,U\pH{if;ii)\^+{^)Jr,ii)ii = exp(iir[$«,t]) , (4.7) 



where 



= K+[ d^x Ka{x,if;ii)^%x) + l[ [ d^xd^x' Kabix,x', if- U)^''{x)^\x')+--- 
J fit ^JJnt 

(4.8) 

is a time translationally invariant quantity. The poly-local sources Kab- encode the state 
of the system at the microscopic time of observation tf = t/2, i.e. the time at which the 
EEV is evaluated, according to Figure 1. It follows that these sources must contribute only 
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for xq = x'q = ■ ■ ■ = tf and therefore be proportional to delta functions of the form 
5{xo — if)5{xQ — if) ■ ■ ■ . For instance, the bi-local source Kab in the double momentum 
representation must have the form 

Kab{x,x',if;k) = If^-^e-'P-^e'P'-^' e"^P^-P'o^'fKai>{p,p',t) , (4.9) 



J J (27r)4 (27r)4 

so that the po and pg integrations yield the required delta functions. Here, it is understood 
that the bi-local Kab sources occuring on the LHS and RHS of (4.9) are distinguished by the 
form of their arguments. We emphasise that Kab{x,x' ,tf;ti) is not a time translationally 
invariant quantity due to the explicit dependence upon if on the RHS of (4.9). In contrast, 
Kabip,p' ,t) is time translationally invariant. 

Notice that we could extend the limits of integration to infinity for the time integrals in 
the expansion of the kernel given in (4.8) also. Nevertheless, for the following derivation, all 
space-time integrals are taken to run over the hypervolume for consistency. We should 
reiterate here that the limits of time integration can be extended to an infinite domain iloo 
in all but the interaction part of the action. We will also suppress the time-dependencies 
of the action S and sources Kab - for notational convenience. 

We now absorb the constant K in (4.8) into the overall normalization of the CTP 
generating functional Z and the local source Ka into a redefinition of the external source Jq. 
Then, Z may be written down as 

Z[Ja,Kab,---] = y"[d$1 exp i(sm + j d^X Ja{x)^''{x) 

+ \jj<^^x Kab{x, x'W{x)^\x') (4.10) 

+ Ijjj d^xdVdV Kabcix,x',x")<^''{x)<^\x')<^^{x") + ■■■^ . 

The Cornwall-Jackiw-Tomboulis (CJT) effective action [37] is given by the following 
Legendre transform: 

p[$a^gab^gafec^...] ^ W[Ja, Kab, Kabc] - J d^X Ja{x)$'' (x) 

d^xdV Kab{x,x') (^^%x)^\x') + ihg''\x,x')'^ 
d^xdVd^x" Kabc{x,x\x") (^'^''{x)^\x')^%x") 

+ mg^''\x,x')'%''\x") - h^g''^^ix,x',x")^ + ■■■ , (4.ii) 

where VV[Ja, Kab, Kabc, ■ ■ ■ ] = —ifi\nZ[Ja, Kab, Kabc, • • • ] is the generating functional of 
connected ensemble Green's functions. We obtain an infinite system of equations: 

^%^) = jjj^^ = m^)) , (4.12a) 

ihg''\x,x') = 2—^^^— - $"(x)$^(x') (4.12b) 
dKab{x,x') 
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6 



5W 



6Kabc{x,x',x") 

6^W 



-h 



6Ja{x)6Jb{x')SJc{x") ' 



3ihg''''''{x,x')^''\x") - 

(4.12c) 



and 



6T 



6^<^{x) 



6r 



-Jaix) - J d^x' Kabix,x')<i>\x') 

d^x' d^x" Kabc{x,x' ,X") (^^\x')^''{x") + ihg^%x\x") 
d\" Kabc{x,x\x"W{x") 



^^Kbix x') - - 

2 J^ab{X,x ) ^ 



"6 



(4.13a) 

' * * * ; 

(4.13b) 
(4.13c) 



where the parentheses (abc) on the RHS of (4.12c) denote cychc permutation with respect 
to the indices o, b, c. 

The above infinite system of equations (4.12) and (4.13) may be simphfied by assuming 
that the density operator p is Gaussian, as we will do later in Section 4.3. In this case, 
the tri- local and higher kernels {Kabc, Kabcd-, ■ ■ ■ ) can be set to zero in (4.13). We may 
then eliminate the three- and higher-point connected Green's functions {Q"'^'^, gabcd^ . . . -j 
as dynamical variables by performing a second Legendre transform 

Pjja^gabj = , (4.i4) 

where the ^'s are functionals of and Q"^^ given by 

6T 



5g 



abc--- 



r^a^gab^gabc^ 



. 



(4.15) 



The effective action is evaluated by expanding around the constant background field 
<I>q(x) = ^°-{x) — h^/'^(f>"-{x), defined at the saddle-point 



<I> = <I>n 



+ Ja{x) + jd*x' Kab{x,x')^l{x') = 0. 



(4.16) 



The result of this expansion is well known [37, 38] and, truncating to order /i^, we obtain 
the two-particle-irreducible (2PI) CJT effective action 



T[^\ ga.] ^ ^[^a] ^ _ Tr, ^Ln, Det.^ Qab + [ G'J - K^b )*Q^'' - Va \ + ^'r2[$", , 

(4.17) 

where a subscript x and the *'s indicate that the trace, logarithm and products should be 
understood as functional operations. The operator is defined by 



G-J{<^^;x,x') 



6^''{x)5^''{x') 



0,-1 



{x,x') + 



+ Kab{x,x') 

5$"(a;)5$^(x') 



+ Kabix,X 



(4.18) 
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where A|^^ ^{x,x') is the free inverse CTP propagator in (3.21) and S*™' is the inter- 
action part of the action. Obviously, aU Green's functions depend upon the state of the 
system at the macroscopic time t through the bi-local source K^b- For the Lagrangian in 
(2.1), we have 



ab 



1 



gr^abc^^ix) - ^\r]abcd^''{x)¥{x)\^ + Kab{x,x') , 



1, Vabc- 



(4.19) 
1 for aU indices 



where rjabc--- = +1 for all indices a = b = ■ ■ ■ 
a = b = ■ ■ ■ = 2 and rjabc— = otherwise. 

The overall normalization (??„") of (4.17) has been chosen so that when Kab{x, x') = 0, 
we may recover the conventional effective action [66] by making a further Legendre trans- 
form to eliminate Q'^^ as a dynamical variable: 



P[^a] ^ = Sm + -Tr,Ln,Det,fcG-,i + 0{h^) 



(4.20) 



Here, Q^^ has been replaced by G^^ and Q""^ is a functional of 

In the CJT effective action (4.17), r2[<I>", ^"''] is the sum of all 2PI vacuum graphs: 



E 



a, 6= 1,2 




'Jab 



12 




(4.21) 



where combinatorial factors have been written explicitly and we associate with each n-point 
vertex a factor of ^ 

and each line a factor of z^"^($"; x, y). The three- and four-point vertices are 

iS'^^)($";x) = -ig7]aaa - iXTjaaaa^" {x) , iS'^f^ {^"^ ] x) = -iXfJaaaa ■ 



(4.22) 



(4.23) 



Upon functional differentiation of the CJT effective action (4.17) with respect to 



Q {x,y), we obtain by virtue of (4.13b) the Schwinger-Dyson equation 



= G-,^($'^;x,y) + ^,b(^^g'^^x,y), 



where 



(4.24) 



-2ih- 



-ih 



5g^\x,y) 



Sab + 



(4.25) 



is the one-loop truncated CTP self-energy. A combinatorial factor of ^ has been absorbed 
into the diagrammatics. 
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Suppressing the <^>" and Q"'^ arguments for notational convenience, the CTP self-energy 
Uabix, y) may be written in matrix form as 



Ii{x,y) -n<(x,y) 
-n>(x,y) -n*(x,y) 



(4.26) 



where n(x, y) and — n*(x, y) are the time- and anti-time-ordered self-energies; and n>(x, y) 
and n<(x,?/) are the positive- and negative-frequency absolutely-ordered self-energies, re- 
spectively. In analogy to the propagator definitions discussed in Section 2 and Appendix A, 
we also define the self-energy functions 

ni(x,y) = n>(x,y) + n<(x,y) = n(x,y) - n*(x,y) = 2nmn(x,y), (4.27a) 

np(x,y) = ^(nR(x,y) + nA(x,y)) = Ren(x,y), (4.27b) 

2iMr(x,y) = n>(x,y) - n<(x,y) = nR,(x,y) - nA(x,2/) , (4.27c) 

which satisfy relations analogous to those described in Appendix A. r(x, y) in (4.27c) is 
related to the usual Breit-Wigner width in the equilibrium and zero-temperature limits. 
The Keldysh representation [see (3.18)] Tlab{x.,y) of the CTP self-energy reads: 



^ab{x,y) 



ni(x,y) nR(x,y) 
nA(x,y) 



(4.28) 



In the limit 



\b{x,y,if;ii) 



0, the Schwinger-Dyson equation (4.24) reduces to 

K:b~\^^y) + Kab{x,y,if;ii) + nab{x,y,if;U) , (4.29) 

in which A~i^{x, y, if;ii) = Gabi^"" = 0; x, y, if, ii) and A|^|j ^(x, y) is the free inverse CTP 
propagator defined in (3.21). We have re-introduced the dependence upon tf and U for 
clarity. Notice that due to the explicit ty-dependence of the bi-local source Kab{x, y,tf; tj) 
in (4.9), the inverse resummed CTP propagator A^^^ {x , y , i f, ii) and the CTP self-energy 
^ab{x , y , t f, ti) are not time translationally invariant quantities. 

In order to develop a self-consistent inversion of the Schwinger-Dyson equation in 
(4.29), the bi-local source Kabix,y,if;ii) is absorbed into an inverse non-homogeneous 
CTP propagator 



^ab ^(2;,yJ/;ii) 



^i^^y) + Kab{x,y,if;ii), 



whose inverse, to leading order in Kab, is the free CTP propagator A'^'"'' 

A^^''\x,y,if-ii) + 0{K^) , 



D^^-\x,y,irrU) 



(4.30) 
x,y,if;ii), i.e. 

(4.31) 



as we will illustrate in Sections 4.3 and 5. The contribution of the bi-local source Kab 
is now absorbed into the free CTP propagator A^''^^{x,y,if;ii), whose time translational 
invariance is broken as a result. The Schwinger-Dyson equation (4.29) may then be written 
in the double momentum representation as 



A-,\p,p'Jf,ii) 



^ab ^(P^P) + '^ab{P,P,if;ii 



(4.32) 
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Since the stationary vacuum |0) has been replaced by the density operator p at the 
microscopic time tf = t/2, we must use the fohowing field-particle duality relation in the 
Wick contraction of interaction-picture fields: 

(0|$(x;t,)at(k,t/;?,)|0) = e'^^e^^e^)*/ . (4.33) 

Here, the extra phase e^^^^^^f arises from the fact that the creation and annihilation oper- 
ators of the interaction-picture field are evaluated at the microscopic time i = 
[cf. (2.12) and (2.18)], whereas the operator a"f(k, tj; tj), resulting from the expansion of 
the density operator p (see Section 4.3), is evaluated at the microscopic time tf = t/2. 
Analytically continuing this extra phase to off-shell energies, we must therefore associate 
with each external vertex of the self-energy Ilab{p,p' ,tf', U), appearing in (4.32), a phase: 

e'Po^f , (4.34) 

where po is the energy flowing into the vertex. This amounts to the absorption of an overall 
phase 

eiiPo-p'o)if (4.35) 

into the definition of the self-energy Ilabip,p' ,tf;ti). We recall from (4.9) that the bi-local 
source Kab{p,p' ,if;ii) already contains this free phase factor. 

Convoluting from left and right on both sides of (4.32) first with the weight func- 
tion {2tt)^61'^\p - p') from (3.50) and then with A'^'''^{p,p' Jfji) and A"^(p,p', t/; t^), 
respectively, we obtain the Feynman-Dyson series 

f f d^n Mn' r\'^n" A^n'" 

A-M,i,:Ui = A-'(P.p'.<,;fO - /-/(l^ jl^ ^ lip ("6) 
X A»'"(j>,,,f,;f.)(2l)''«!-"(5 - q')IlM,<l",ir,i,K2ir)'4'>{q" - «"') A*!?'",?', (/; «,) , 

where Sl-^\p — p') is defined in (3.49). Because of the form of 5t{po — p'q) in (3.50), we 
see that this series does not collapse to an algebraic equation of resummation, as known 
from zero-temperature field theory. As we will see in Section 4.2, one cannot write down a 
closed analytic form for the resummed CTP propag ator A"-''{p,p',if;ii). 

Given that 5t satisfies the convolution in (3.52), the weight functions may be absorbed 
into the external vertices of the self-energy Ilabip,p' ,if;ii); see Section 9. The Feynman- 
Dyson series may then be written in the more concise form 

A"^(p,p',V;t,) = A'''-'{p,p',if;U) (4.37) 

~ llwfWf^'^ "'(P, Q, h mcdiq, q, if, U)A''\q',p', if-U) . 

Note that for finite t, 5t{po — p'q) is analytic in the limit po — )• p'q. As we shall see 
in Section 6, the systematic incorporation of these finite-time effects ensures that the 
perturbation expansion is free of pinch singularities. 
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4.2 Applicability of the Gradient Expansion 

Here, we will look more closely at the inverse relation (3.47) that determines the resummed 
CTP propagator. We will show that a full matrix inversion may only be performed in the 
thermodynamic equilibrium limit. Hence, the application of truncated gradient expansions 
and the use of partially resummed quasi-particle propagators, particularly for early times, 
become questionable in out-of-equilibrium systems. 
We define the relative and central coordinates 

+ 

such that 



R% = - y^ X^y = ; ' , (4.38) 



— Xt^y + ^Rxy 1 — 2^^y ' (4.39) 

We then introduce the Wigner transform (see [51]), namely the Fourier transform with 

tion F{R,X) is 



respect to the relative coordinate Rxy only. Explicitly, the Wigner transform of a func- 



F{p,X) = Jd^Re'P-^F{R,X) . (4.40) 

The resummed CTP propagator A°'^{x, y) respects the inverse relation in (3.47). Here, 
we suppress the tf- and tj-dependence of the propagators for notational convenience. In- 
serting into (3.47), the Wigner transforms of the resummed and inverse resummed CTP 
propagators, the inverse relation takes the form 



11 (2^)' (2vr)4 ' 



^ip,.y ^^c(2^)4^(4) _ ^ 



(27r)4 (2^)4 

where the domain of integration is restricted to be in the range [— t/2\. 

In the case where deviations from homogeneity are small, i.e. when the characteristic 
scale of macroscopic variations in the background is large in comparison to that of the 
microscopic single-particle excitations, we may perform a gradient expansion of the inverse 
relation in terms of the soft derivative = d/dXxy,^- Writing = Xxy + Rzy/'^ 

and Xzy = Xxy — Rxz/"^ and after integrating by parts, we obtain 



// 



(27r)4 (27r)4 

d^P dV .-.p...y.v.c,o.^4.(4) 



X A-i(p,X)exp[-^(5-5j - S-^)]a'"^(p',X) 



(27r)4 (2^) 



^ e-^P-- -y r^^{2TTf5r (p - p') , (4.42) 



where X = Xxy 



We now define the central and relative momenta 



^ + P'" ^ = - p'M ^ (4.43) 
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which are the Fourier conjugates to the relative and central coordinates, and X^, 
respectively. It follows that 

= + , p''^ = g'^ - ■ (4.44) 

We may then Fourier transform (4.42) with respect to Rxy to obtain 

{0;,x + 20+,,)]{A-,Hq + Q/2,X)}{A''{q-Q/2,X)} 

Va'(27: f6l'\Q) , (4.45) 



X exp 



-iQ-X 



{2^Y 

where, following [42, 43], we have introduced the diamond operator 

<>tx{A]{B] = \{A, 3)1^ (4.46) 
and {A, denote the symmetric and anti-symmetric Poisson brackets 

^^J- ^ W^dx.^dx^W,- ^'-''^ 

For t > 0, we may perform the integral on the RHS of (4.45), yielding 

J4 



d Q -iQ-X /n_N4 r(4) 



(2vr) 



X exp 



4 



'^■^ {2T:f5r{Q) (4.48) 
(0,> + 20^,x)]{A;,n'? + Q/2,X)}{A^'^(g-Q/2,X)} = r^^e{t-2\X,\), 



In the above expressions, of particular concern is the (}q ^ operator, where the relative 
momentum Q and the central coordinate X are conjugate to one another. Thus, if the 
derivatives with respect to X are assumed to be small, then the derivatives with respect 
to Q must be large. In this case, all orders of the gradient expansion may be significant, 
so it is inappropriate to truncate to a given order in the soft derivative d'^ . 

As t — )■ oo, we have the transition 5\^\q) — )■ 5'^'^\Q) and (4.48) reduces to 



q,X 



{A-J{q,X)}{A"^{q,X)] = r^^. (4.49) 



Even for these late times, we can perform the matrix inversion exactly only if we truncate 
the gradient expansion in (4.49) to zeroth order. However, such a truncation appears 
valid only for time-independent and spatially homogeneous systems. Employing a suitable 
quasi-particle approximation to the Wigner representation of the propagators, it can be 
shown [55] that this inversion may be performed at first order in the gradient expansion. 
However, off-shell contributions are not fully accounted for in such an approximation. 

In conclusion, a closed analytic form for the resummed CTP propagator may only be 
obtained in the time- independent thermodynamic equilibrium limit. The truncation of the 
gradient expansion may be justifiable only to the late-time evolution of systems very close 
to equilibrium, even for spatially homogeneous thermal backgrounds. A similar conclusion 
is drawn from different arguments in [67]. 
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4.3 Non-Homogeneous Free Propagators 

Unlike the resummed CTP propagator, the free CTP propagator can be derived analyti- 
cally, even in the presence of a time-dependent and spatially inhomogeneous background. 
The non-homogeneous free propagator will account explicitly for the violation of space-time 
translational invariance. Our derivation relies on the algebra of the canonical quantization 
commutators of creation and annihilation operators described in Section 2. Subsequently, 
we make connection of our results with the path-integral representation of the CTP gen- 
erating functional in (4.10). Finally, we introduce a diagrammatic representation for the 
non-homogeneous free CTP propagator. 

We note that the derived propagators are 'free' in the sense that their spectral structure 
is that of single-particle states, corresponding to the free part of the action (see Section 3). 
Their statistical structure, on the other hand, will turn out to contain a summation over 
contributions from all possible multi-particle states. The time-dependent statistical distri- 
bution function appearing in these propagators is therefore a dressed object. This subtle 
point is significant for the consistent definition of the number density in Section 7 and the 
derivation of the master time evolution equations in Section 8. 

The starting point of our canonical derivation is the explicit form of the density oper- 
ator p. We relax any assumptions about the form of the density operator and take it to be 
in general non-diagonal but Hermitian within the general Fock space. We may write the 
most general interaction-picture density operator at the microscopic time tf = t/2 as 

p{if,ii) = Cexp - j dHk, W^io(ki : f))a\]^i,tf) - j dHk; VFoi(0 : K)a{]^^,tf) 

- 11 dHk, dHk. W^n(ki : k[)a\ki,if)a{k[,if) - ■■■ (4.50) 

1^1 / n^nkjindnk;) Wnm{{k,} : m})lla\k,,if)l[am,if) 

\i=l ^ ^ 7 = 1 ^ i=l 1 = 1 



1 1 

n 



where the constant C can be set to unity without loss of generality. The complex-valued 
weights VFnm({k}n; {k'}m) depend on the state of the system at time tf = t/2 and satisfy 
the Hermiticity constraint: 

Wnm{{k}n : {k'}„) = t^4„({k'}^ : {k}„) . (4.51) 

The density operator p may be written in the basis of momentum eigenstates by multiplying 
the exponential form in (4.50) by the completeness relation of the basis of Fock states at 
time tf = t/2: 

I = |0)(0| + f; [dU^?j\{p}i,if){{pUif\ , (4.52) 

where |{p}£,t/) is the multi-mode Fock state |pi,t/) ® |p2,t/) ® ••• ® |p£,t/). This 
usually gives an intractable infinite series of n-to-m-particle correlations. Taking all weights 
W„m({k}ri, : {k}m) to be zero if n + m > 2, i.e. taking a Gaussian-like density operator, 
it is still possible to generate all possible n-to-m-particle correlations. In Appendix D, 
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we give the expansion of the general Gaussian-hke density operator, where only sufficient 
terms are included to help us visualize its analytic form. 

Exploiting the remaining freedom in the commutation relations in (2.18), we may 
account for our ignorance of the series expansion of the density operator by writing the 
EEVs of two-point products of interaction-picture creation and annihilation operators as 

{a\p',if;ii)a{p,if;ii))t = 2 (p, p')/(p, P , i) > (4.53a) 
{a{p',if;ii)a{P,if^ti))t = 2 ^{p,p')g{p,p' ,t) . (4.53b) 

The energy factor 2(f(p,p'), having dimensions E^, arises from the fact that the 'number 
operator' {p,i;ii)a{p,i;ii) of quantum field theory has dimensions E~^, i.e. it does not 
have the dimensions of a number. Bearing in mind that the density operator is constructed 
from on-shell Fock states, a natural ansatz for this energy factor is 

^(p,p') = VE{p)E{p') . (4.54) 

The complex-valued distributions / and g have dimensions E~'^ and satisfy the identities: 

/(P,p',t) = r(p',P,t) , (4.55a) 
9{p,P,t) = g{p,p,t). (4.55b) 

We refer to / and g as statistical distribution functions. In particular, we interpret the 
Wigner transform 

n(q,X,t) = y|^e^Q-X/(q+Q,q_Q,i), (4.56) 

having dimensions E^, as the number density of particles at macroscopic time t in the 
phase-space hypervolume between q and q + d^q and X and X + d^X. Notice that 
n(q, X, f) is real thanks to the Hermiticity constraint (4.55a). 

The EEV of the two-point product {a{p,if;ii)a^{p',if;ii))^ follows from the definition 
(4.53a) and the canonical commutation relation in (2.18), giving 

{a{p,if;ii)aHp',if;ii))^ = {27rf2E{p)5^^Hp - p') + 2 ^(p, p')/(p, p', t) . (4.57) 

Hermitian conjugation of (4.53b) yields 

{a\p,if;ii)a\p',ij;ii))^ = 2 ^{p,p')g* {p,p' ,t) . (4.58) 

Note that (4.53), (4.57) and (4.58) are consistent with the canonical quantization rules in 
(2.14) and (2.17). 

When the linear terms in the exponent of the density operator p in (4.50) are non-zero, 
we may consider the EEVs of single creation or annihilation operators 

{aip,if;ii))^ = 2,/Eip)w{p,t) . (4.59) 
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In this case, we may define the connected distribution functions 

/con(P,p',i) = /(P,p',i) - W{p,t)w{p,t) , 
5con(P,p',i) = g(.P,p',t) - W{p,t)w{p',t) , 



(4.60a) 
(4.60b) 



which obey the same symmetry properties given in (4.55). 

We are now in a position to derive the most general form of the double momentum 
representation of the non- homogeneous free CTP propagator, satisfying the inverse rela- 
tion (3.33). Proceeding as in Section 3.2, we make the following ansatz for the most general 
solution of the Klein-Gordon equation in the double momentum representation: 



, 0, ab 



{p,p,if;ii) 



[P 



M2 + te)-^ 



-i2-Ke{-pQ)6{p^ 
_(p2 _ - 



M2) 
-1 



i27r\2po\^/^5{p^ - M^)f{p,p',t)e'^P'^-P'oy'f27T\2pX^^S{p' 



12 



M 



{2T:f6^^\p - p') 
(4.61) 



1 1 
1 1 



which we confirm by evaluating the EEVs directly, using the algebra of (4.53). 

In (4.61) the phase factor e*^^° ~Po)*/ arises from the fact that the creation and anni- 
hilation operators appearing in the Fourier transform of the field operator given in (2.12) 
are evaluated at the time t = 0. The density operator, on the other hand, is evaluated at 
the time tf. As a consequence, in the evaluation of the EEV, we have, for instance, 

(at(p',0;tXp,0;t.)), = 2 ^(p, p')/(p, p', t)e*[^(P) " , 



(4.62) 



which directly results from (4.53a). 

The form of the function f{p,p',t) is 

f{p,p',t) = e{po)e{p'o)fip,p',t) + 9{-po)ei-p',)r{-p,-p',t) 

+ e{po)e{-p'o)g{p,-p',t) + e{-po)e{p'o)g*{-p,p',t) 



(4.63) 



The function / satisfies the relations: f{p,p', t) = f{—p', —p, t) = f*{—p, —p', t), consistent 
with the properties in (A. 4). It also contains all information about the state of the ensemble 
at the macroscopic time t. For this reason, we refer to / as the ensemble function. 

In the double momentum representation, the retarded and advanced propagators are 

1 



^R(A) 



{p,p) 



(PO + (- 



le 



M2 



(2vr)^5(^)( 



p — p 



(4.64) 



The Pauli-Jordan A^{p,p'), Hadamard Ai{p,p' ,if;ii) and principal-part Ap{p,p') propa- 
gators become 

A\p,p') -- 
A\{p,p\if-U) -- 



i27re(po)<5(/ - M'^){2^f5^^\p - p') , (4.65a) 
z27r(5(p2 - M'^){2TTf5^^\p - p') (4.65b) 
z27r|2po|^/^5(/ - M2)2/(p,p',t)e^(Po-Po)*-f2^|2p'o|i/25(p'2 - M^) , 



a^(p,pO =V 



p^ 



M2 



{2T:f6^^\p - p') 



(4.65c) 
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Propagator 


Double Momentum Representation 


(Dyson) 


iA%^Jp,p',if;ii) = „ , (27r)V^)(p p') 

+ 2tt\2po\^/^5{p'^ - M2)/(p,p',t)e*(P«-Po)*/27r|2p'o|V25(p'2 _ ^2) 


+ (— )ve-freq. 

\A/ 1 crn'j'TTi fi n 
vv igiitiiicxii 


iA0(<)(p,y,t7;t,) = 2^0(+(-)po)'^(p' - M2)(2^)4<5(4)(p - p') 
+ 27r 2po ' '5(p - M^)f{p,p',t)e^^P° Poni2T:\2pQyi^5{p^ - M^) 


Retarded 
(Advanced) 


'^"l*)*"-"'' = (p„+(-)..P-p^-M2(2'''*'''(''-f') 


Pauli-Jordan 


iAO(p,p') = 27re(po)<5(p2 - M'^){2tiY5^*\p - p') 


Hadamard 


iA\{p,p',tf]ti) = 27r5(p2 _ M2)(27r)4M4)(p - p') 

+ 27r 12^01^/^-5(^2 - M2)2/(p,p',f)e*(Po-Po)*/27r|2p'o|i/25(p'2 - M^) 


Principal-part 


^A" (p,p') = P-^^^(2vr)^<5(^)(p - p') 



Table 2: The full complement of non-homogeneous free propagators, where f{p,p',t) is 
the ensemble function defined in (4.63) and V denotes the Cauchy principal value. 



Thus, at the tree level, only the Hadamard correlation function A^{p,p' , if;ii) depends ex- 
plicitly on the background and macroscopic time t, through the ensemble function f{p,p', t) 
in (4.63). This is a consequence of the causality of the theory, as we would expect from 
the spectral decomposition (2.22) of the retarded and advanced propagators A^^^^(x,y) 
in terms of the canonical commutation relation (2.14). Notice that the complex phase 
factor e'^^^" ~ 'P'o^^f has only appeared in the Hadamard propagator (4.65b) and so it does 
not spoil causality. Beyond the tree level, the background contributions are expected to 
modify the structure of the Pauli- Jordan and causal propagators, according to our discus- 
sion of the Kallen-Lehmann spectral representation in Section 2. The full complement of 
non- homogeneous free propagators is listed in Table 2. 

It would be interesting to establish a connection between these canonically derived 
non-homogeneous free propagators and those derived by the path-integral representation 
of the CTP generating functional Z. This will be achieved through the bi-local source 
Kab via the tree-level Schwinger-Dyson equation in (4.30). The role of the bi-local source 
Kab will be illustrated further, when discussing the thermodynamic equilibrium limit in 
Section 5. 

We proceed by replacing the exponent of the CTP generating functional Z in (4.10) by 
its double momentum representation. Subsequently, we may complete the square in this 
exponent by making the following shift in the field: 

$-(p) = <I>'«(p) - A°''^,(p)J^(-p) , (4.66) 
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». , 



• h 



Figure 3: The Feynman-diagrammatic interpretation of the non- homogeneous free CTP 
propagator iA^'"'^{p,p' ,if;ii) for the real scalar <I>, where the double line represents 
momentum-violating coupling to the thermal background through the bi-local source 
Kabip,p',if;ii)- 



where A^'"'^{p) is the free vacuum CTP propagator in (3.44) in which the ensemble func- 
tion / of (4.61) is set to zero. Notice that the normal-ordered contribution does not 
appear in A^'°-^{p), as it is sourced from the bi-local term Kab- Upon substitution of (4.66) 
into (4.10), the CTP generating functional Z takes on the form 



Z[J,K,t] = Z°[0,K,t]exp\iS 



lint 



l_6_ 

i 5Ja 



[d$ 



/ai \ i_ [ f d^p d'^p' 



(4.67) 



r{p)A^','{p)K,b{p,p',m'\-p') + ^"'{p)KUp,P,t)A'>^MpV''{-p) 
- r{p) (A0,(p)(2vr)^5W(p - p') - A'',^{p)K,,{p,p\t)A'^\{p'))j\-p') 



For the Lagrangian in (2.1), the cubic self-interaction part {—g^^) of the action may be 
written down explicitly as 



-w 

3! 



d^pi d^p2 d^p3 



(27r)4 (27r)4 (2^)^ 



llabc{2Tl)H[^\pi 



+ P2 + P3) 



1 



6 



C iS 



int 



1 6 

i 5Ja 



(4.68) 



i SJaipi) i 5Jb{p2) i SJcips) 

Hence, in the three-point vertex, the usual energy-conserving delta function has been re- 
placed by 5t, defined in (3.50), as a result of the systematic inclusion of finite-time effects. 
This time-dependent modification of the Feynman rules is fundamental to our perturbative 
approach to non-equilibrium thermal field theory and will be discussed further in Section 9 
in the context of a simple scalar model. 

In (4.67), the remaining terms linear in the external source J yield contributions to the 
free propagator proportional to upon double functional differentiation with respect to J. 
As we shall see in Section 5, these contributions may be neglected. Employing (3.25), we 
find that the non-homogeneous free CTP propag ator A°'"^(p,p',t/;ti) may be expressed in 
terms of the free vacuum CTP propagator A'''"*(p) and the bi-local source Kab{p,p', if, ii) 
as follows: 



iA''''\p,p'rtr,ii) 

where 



iA°''^''(p)(2^)45(4)(p - p') + iA''^-'^ip)iK,a{p,P,ir,U)iA''^^\p') 



Kab{p,p',if,k) = e^(Po-Po)*/i^„,(p,p',t) 



(4.69) 
(4.70) 
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The form of the free CTP propagator A^'"'^{p,p' ,if;ii) in (4.69) is consistent with a per- 
turbative inversion of (4.30) to leading order in the bi-local source Kab- It is also consistent 
with the canonically derived form of the non-homogeneous free propagators in (4.61). 

The result in (4.69) may be interpreted diagrammatically, where the non-homogeneous 
free CTP propagator iA^'°'^{p,p' , if, ti) for the real scalar <1> is associated with the Feynman 
diagram displayed in Figure 3. The bi-local source Kab{p,p' ,tf;ti) plays the role of a three- 
momentum- violating vertex that gives rise to the violation of translational invariance, thus 
encoding the spatial inhomogeneity of the background. 

5 The Thermodynamic Equilibrium Limit 

In this section, we derive the analytical forms of the free and resummed CTP propagators 
in the limit of thermal equilibrium. The results of this section are of particular importance 
for the discussion of pinch singularities in Section 6. We also show the connection between 
the equilibrium Bose-Einstein distribution function and the bi-local source Kab introduced 
in Section 4.3. 

In the limit of thermal equilibrium, the density operator p is diagonal in particle 
number, so all amplitudes except Wu vanish in (4.50). In this limit, the general density 
operator p, given explicitly in (D.l), reduces to the series 

p = |0)(0| + IJdU^.dUy^ ((2^)32i?(ki)5(3)(ki - k'l) - VFn(ki ik'i) (5.1) 
+ IJdU^, iyn(ki :qi)I^ii(qi ik'i) + ••• )|ki)(k;| 

+ IJ ■■■ ydHk, dHk, dnk. dHk. ((27r)32i^(ki)<5(3)(ki - k[) - T^n(ki : k'^) 

+ dHq, I^n(ki : qi)I^ii(qi : k'l) + • • • ) ( {27Tf2E{k2)6^^\k2 - k'^) 

- Wn{k2 : k'2) + dHq, W^ii(k2 : q2)I^ii(q2 : k'2) + • • • ) |ki, k2)(k;, + • • • , 

where the time arguments in the multi-mode Fock states have been omitted for notational 
convenience. In the equilibrium limit, the statistical distribution function g{p,p',t) is 
trivially zero. Instead, /(p, p',t) is calculated from (4.53) and takes the form of the series 

2£{p,p')fip,p',t) = (27r)32S(p)<5(3)(p - p') - H^n(p : p') (5.2) 

+ dHq Wii{p : q)II^ii(q : p) + ■ ■ ■ 
+ y"dnq ( {27Tf2E{p)6^'\p - q) - iyii(p : q) + . . . ) 

X ((27r)32£;(q)<5(3)(q - p') - Wuiq-.p') + •••)+•••, 

where disconnected parts have been cancelled order-by-order in the expansion by the nor- 
malization Trp in (4.3). The factor £{p,p') is defined in (4.54). 
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The equilibrium density operator peq must also be diagonal in the momenta and is 
thus obtained by making the replacement: 



M/ii(k:k') ^ (3E{k){2iTf6^^\k - k') (5.3) 
in (5.1), where /3 = 1/T is the inverse thermodynamic temperature. In detail, we find 



Peq = |0)(0| + ^ - / ... m(dnk,/^(k,))(g)|k,)(g)(k, 

n=l 1=1 i=l i=l 



(5.4) 



where the amplitudes are the Boltzmann distributions fpik) = e ^^C'). This last expres- 
sion of peq can be shown to be fully equivalent to the Gaussian form 



Peq = exp(^-/3ydnkS(k)at(k)a(k)^ , 



(5.5) 



which corresponds to the standard Boltzmann density operator 

Peq = e-^^° , (5.6) 

where is the free part of the interaction-picture Hamiltonian. Note that our convention 
for the normalization of the density operator p, including p^q, is chosen so that the canonical 
partition function Z{f3) = Tr e~^^° appears explicitly in the definition of the EEV in (4.3). 

We may now substitute the limit (5.3) into the series expansion of the statistical 
distribution function / in (5.2). Using the identities of summation 

oo 



n=0 n=l 

we find the following correspondence in the equilibrium limit: 

/(p,p',t) /eq(p,p') = i27Tf5^'\p - p')/B(i^(p)) , (5.8a) 

eq. ^ ' 

g{p,p',t) <7eq(p,p') = 0, (5.8b) 

eq. 

where 

/b(po) = ^^3^ (5-9) 

is the Bose-Einstein distribution function. The equilibrium statistical distribution func- 
tions in (5.8) depend only on the magnitude of the three-momentum p via the on-shell 
energy E{p). This is a consequence of the homogeneity and isotropy implied by thermo- 
dynamic equilibrium. Moreover, the multiplying factor ^(^^(p — p') on the RHS of (5.8) 
necessarily restores translational and rotational invariance. 

It is well-known that the pinch singularities present in perturbative expansions cancel 
in the equilibrium limit (see Section 6) and we can safely take the limit t — t- oo throughout 
the OTP generating functional, as we should expect for a system with static macroscopic 
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properties. Working then in the single momentum representation, we obtain from (3.42) 
and (3.43) the free equihbrium CTP propagators 



HP 



+ ie 



+ 2^/b(|po|)<J(/ - 



(5.10a) 



iAlip) = 27r{e{po) + fB{\po\)]S{p' - M') = 27re(po)(l + /b(po))5(p 



(5.10b) 

iAlip) = 27r(e{-po) + fB{\po\))s{p^ - M^) = 27re{po)fB{po)S{p^ - M^) . (5.10c) 

The form of the Wightman propagators written in terms of the signum function e(po) prove 
very useful in the calculation of loop diagrams, as detailed in Appendix B. 

Returning to the free CTP generating functional Z in (4.67), it follows from the results 
above that in equilibrium the bi-local source Kab{p,p' , tf, ti) must be proportional to a four- 
dimensional delta function of the momenta, i.e. 



Kab{p,p',if;ii 



eq. 



(5.11) 



In addition, K^'^{p) must satisfy 



1 1 
1 1 



(5.12) 



Solving the resulting system of equations, keeping terms to leading order in e, and noting 
that K'^{p) should be written in terms of the three- momentum only, we find 



K2{p) = 2ieh{E{p)) 



1 1 
1 1 



By virtue of the limit representation of the delta function 

1 e 



5{x) 



hm 



(5.13) 



(5.14) 



we can verify that we do indeed recover (5.13) and the correct free CTP propagator by 
means of (4.69). We also confirm that the terms linear in J remaining in (4.67) may safely 
be ignored, since they yield contributions to the free propagator proportional to ~ 
upon double functional differentiation with respect to J. 

Alternatively, interpreting the Boltzmann density operator in (5.6) as an evolution 
operator in negative imaginary time and using the cyclicity of the trace in the EEV, we 
derive the Kubo-Martin-Schwinger (KMS) relation (see for instance [68]) 



y^x-y) = A<(x"-yVi/3,x-y) . 
In the momentum representation, the KMS relation reads: 

A>(p) = e^P^A^ip) , 



(5.15) 



(5.16) 
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which offers the final constraint on f{p) in (3.43): 

Hp) = 0ipo)fB{po) + 9{-Po)fB{-Po) = /b(IpoI). (5.17) 
Furthermore, the KMS relation also leads to the fluctuation-dissipation theorem 

Ai(p) = (l + 2/b(po)) A(p) , (5.18) 



relating the causality and unitarity relations in (2.24) and (2.28). Subsequently, by means of 
the KMS relation, we may write all propagators in terms of the retarded propagator Ar(p): 



Re Af(p) 
Im Af(p) 
A>(p) 
A<(p) 



ReAR(p) , 

e(po)(l + 2/B(|po|))lmAR,(p) , 
2ieipo)(9{po) + /B(|po|))lmAR(p) , 
2ie{po)(9{-po) + /B(|po|))lmAR(p) . 



(5.19a) 
(5.19b) 

(5.19c) 

(5.19d) 



In the homogeneous equilibrium limit of the Schwinger-Dyson equation in (4.29), the 
inverse resummed CTP propagator is given by 



^atiP) 



Al'~\p) + UM ■ 



(5.20) 



In the absence of self-energy effects, the free equilibrium CTP propagator is obtained by 
inverting the equilbrium limit of (4.30): 



K,-\p) + K2{P) ■ 



(5.21) 



Knowing that ~ e from (5.13), the inversion of (5.21) can be done perturbatively 

to leading order in K^{p), in which case the expression (4.69) gets reproduced. Beyond 
the tree- level, however, the contribution from the bi- local source K'^{p) may be neglected 
next to the self-energy term Iiab{p) and the inverse resummed CTP propagator is explicitly 
given by 



(5.22) 



p2 - + n(p) -n<(p) 
-n>(p) + m2 - w{p) 

In this equilibrium limit, (5.22) may be inverted exactly, yielding the equilibrium resummed 
CTP propagator 



m2 + RenR(p) + ImnR(p) 



p2 - m2 + n*(p) 
-n>(p) 



-n<(p) 

+ M2 - U{p) 



(5.23) 



In this single momentum representation, the self-energies satisfy the unitarity and causality 
relations 

ni(p) = n>(p) + n<(p) = U{p) - U*{p) = 2ilmU{p) , (5.24a) 
2iMr(p) = U^ip) - n<(p) = nR(p) - UAip) = 2iImnR(p) , (5.24b) 
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where T{p) is the Breit-Wigner width, relating the absorptive part of the retarded self- 
energy IIr{p) to physical reaction rates [69, 70]. Notice that the KMS relation (5.15) leads 
also to the detailed balance condition 

n>(p) = e^P^U^ip) . (5.25) 

Given the relations in (5.24), we find, in compliance with (5.19), an analogous set of 
relations for the elements of the OTP self-energy: 

ReU{p) = RenR(p) , (5.26a) 

Imn(p) = e(po)(l + 2/B(|po|))lmnR(p) , (5.26b) 

n>(p) = 2ie{po){e{po) + /B(|po|))lmnR(p) , (5.26c) 

n<(p) = 2ieipo){ei-po) + /B(|po|))lmnR(p) . (5.26d) 

Ignoring the dispersive parts of the self-energy, we expect to recover the free OTP 
propagators given in (5.10) in the limit lmll(p) — )• e = 0^. This limit is equivalent to 

ImnR(p) ^ eR = e{po)e. (5.27) 

Expressing the equilibrium resummed CTP propagator in (5.23) in terms of the retarded 
absorptive self-energy ImnR(p), we can convince ourselves that we do indeed reproduce 
the free equilibrium CTP propagators (5.10) in the limit (5.27). 

In Appendix B, we discuss the correspondence of the results of this section with the 
Imaginary Time Formalsim (ITF), clarifying the analytic continuation of the imaginary- 
time propagator and self-energy. 



6 Perturbative One-Loop Expansion without Pinch Singularities 

Pinch singularities normally spoil the perturbation expansion of non-equilibrium Green's 
functions [71]. These pathologies arise from products of delta functions with identical argu- 
ments. In this section, we explicitly demonstrate, in contrast to [72], how pinch singularities 
do not occur in our perturbative approach. 

The avoidance of pinch singularities is achieved by the systematic inclusion of finite- 
time effects, which results in finite upper and lower bounds on time domains, as shown in 
Figure 1. Thus, our approach is different from the semi-infinite time domains employed 
in the existing literature [73, 74]. In particular, in our approach the finite-time effects are 
not included a priori and the removal of the pinch singularities is achieved without any 
ad hoc regularization or obscure resummation. Most importantly, proper consideration 
is given to the dependence upon the time of observation t. As noted in Section 4.3, the 
statistical distribution functions appearing in the non-homogeneous free propagators are 
dressed objects, evaluated at the macroscopic time t. Hence, the role of the Feynman- 
Dyson series is to dress the spectral structure of the propagators only. The evolution of 
the statistical distribution functions is determined by master time evolution equations to 
be derived in Section 8. 
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For early times, the appearance of pinch singularities is prevented by microscopic 
violation of energy conservation. At late times, the thermalization of the system ensures 
that pinch singularities cancel by virtue of the KMS relation (5.15). In this section, we 
demonstrate this cancellation explicitly at the one-loop level and describe the absence of 
pinch singularities in our approach for all times. 

In order to proceed perturbatively, we truncate the Feynman-Dyson series in (4.37) 
to leading order in the couplings and set the bi-local source Kab to zero. This corre- 
sponds to keeping free CTP propagators A^'"'''{p,p' ,if;ii) and one- loop CTP self-energies 
U^^jj {p,p' ,if;ii) on the RHS of (4.37). We then have the one-loop-inserted CTP propagator 



AW'"''(p,p',t^;t~) = A''^\p,p',if,U) 



It will prove convenient to work in a mixed CTP-Keldysh basis by inserting the transfor- 
mation outlined in (3.18) between the external legs and self-energies on the RHS of (6.1). 
In this mixed CTP-Keldysh basis, we may rewrite (6.1) as 

A('^'''\p,p',if,U) = A^''^\p,p',if-U) 

~\jjwf$^' A°;^(p,g,^>;^;)^«(g,g^^>;^OA°:if (g^p^^/;^ , (6.2) 



where, making use of the relations in (A. 6) 
A'i',,{p,q,tf,ii) -- 



AF{p,q,tf;ti) A<{p,q,tf;ti) 
^>{p,q,if;ii) -A^{p,q,if;ii) 



1 1 

-1 1 



^RiP,q,if;ti) ATx{p,q,if;ii) + 2A<(p, tj) 
^R{p,q,if;ii) -A^{p,q,if]ti) + 2A>(p, g, t/; t^) 



3.3) 



and 



A 



1 -1 




1 1 





AF{p,q,if;ii) A<{p, q,if;ii) 
^>{p,q,if;ii) -A^{p,q,if;ii) 



Ax{p,q,if;ii) 
AA{p,q,if;ii) + 2Ay{p,q,if;ii 



(6.4) 



AA{p,q,if;ii) 
-AA{p,q,if;ii) + 2A<(p, g, t/; 



In the same mixed CTP-Keldysh basis, the one-loop-inserted Feynman A^\p,p' , if, ii) 
and positive-frequency Wightman A^\p,p' ,if;ii) propagators may be written in the fol- 
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lowing forms: 

+ A^^{p,q)U^^\q,q',if;ii)Aliq',p',if;U) + A%ip,q,if,U)U^l\q,q' ,if;U)Al{q' ,p')^ , 

(6.5a) 

A^^\p,p',if,k) = AUp,p',if,U) - 11^^ (^Ai{p,q)U^^\q,q'Jj;k)Al{q',p') 

+ A'^{p,q)U^^\q,q',if;ii)Aliq',p',if;ii) + Alip,q,if,U)U^l\q,q' ,if;U)Al{q' ,p')^ , 

(6.5b) 

where we have used the identities 

ni(p,p',t/;ti) + IlYi{p,p,if;ii) + IIa{p,P ,if;ii) = 2Il{p,p ,if;ii) , (6.6a) 
Ui{p,p',if;ii) + nn{p,p',if;ii) - UA{p,p',if;U) = 2n>{p,p' ,if;U) , (6.6b) 

which can be derived from relations between the self-energies analogous to those in (A. 6). 

In the equilibrium limit, where translational invariance is restored, the expansions in 
(6.5a) and (6.5b) can be rewritten in the single momentum representation as follows: 

4'\p) = AUp) (6.7a) 

- (AO(p)n(i)(p)AO(p) + Al{p)U^l\p)Al{p) + A%{p)U^l\p)Al{p)) , 

A^^\p) = Alip) (6.7b) 

Ai{p)U^^\p)Al{p) + AO(p)njJ)(p)AO(p) + Al{p)U^l\p)Al{p)) , 



where the single momentum representations of the equilibrium propagators and self-energies 
are those discussed in Section 5. Using the identities in (5.19) and (5.26), derived by virtue 
of the KMS relation (5.15) and detailed balance condition (5.25), we can rewrite (6.7a) and 
(6.7b) entirely in terms of the retarded functions Ap^(p) and Il^\p). In this way, it can 
be shown that undesirable terms proportional to <5^(p^ — M^), which would signify pinch 
singularities, cancel exactly. We are then left with terms depending upon the following 
product: 

ReAO,(p)ImAO(p), (6.8) 
where ^ 

with V denoting the principal value integral, and 

ImA^(p) = -Tre{po)6{p'^ - M^) . (6.10) 

In order to show that (6.8) is also free of pinch singularities, we consider the limit 
representation of the Cauchy principle value: 

V- = lim . ^ ^ (6.11) 
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and the limit representation of the delta function in (5.14). We may then show that 

ReAO(p)ImA^(p) = |e(po)5'(p' - M^) , (6.12) 

where 5'{x) is the derivative of the delta function, satisfying 

+00 

dx 5'{x)y{x) = -y'{0) , (6.13) 

-00 

provided the function y{x) is analytically well-behaved at x = 0. Consequently, the ex- 
pression (6.8) is free of pinch singularities. 

Next, we use (6.12) and the free CTP propagators from (5.10) for the remaining terms 
in (6.7a) and (6.7b). This allows us to find a closed form for the one-loop-inserted equilib- 
rium Feynman and Wightman propagators 

A«(P) = : + + 2^^M\P0\)(S{P' - M^) + n*(i)(p)5'(/ - M^) 

(p^ — + ie) ^ 

(6.14a) 

(1) ^>\p) 



[p2 _ M2 ief 



27ri(e(±Po) + /b(|po|))(<5(/ - m2) + n*«(p)J'(/ - M^)) , 



where 



n(p) = RenR(p) + «e(po)ImnR(p) (6.15) 

is a self-energy-like function, bearing by itself no direct physical meaning. In (6.14b), we 
have included the negative- frequency Wightman propagator A^'*(p) for completeness. The 
results in (6.14) are consistent with the properties and relations in Appendix A and reduce 
to the expected form in the zero-temperature limit. In conclusion, we have shown that the 
one-loop perturbative expansion of the equilibrium CTP propagator is well-defined and 
that potential pinch singularities cancel exactly in the thermodynamic equilibrium limit as 
a result of the KMS relation. 

In non-equilibrium situations, the KMS relation does not hold. It would appear there- 
fore that the application of perturbation theory in the regime of out of thermodynamic 
equilibrium would be plagued by pinch singularities. However, for finite i, the usual 
energy-conserving delta functions have been replaced by the time-dependent energy-non- 
conserving weight functions 5t, defined in (3.50). Thus, for finite t, pinch singularities do 
not occur, even for general statistical distribution functions, since products of 5ts with 
identical arguments remain analytic. 

In the limit t — )• c«, one would naively expect pinch singularities to reappear, since 
6t{x) —7- 5{x) as given in (3.51). However, in the same limit, our boundary conditions are 
specified in the infinitely distant past. Hence, our system must have reached thermody- 
namic equilibrium. Thus, in the limit i — t- cxd, the time-dependent statistical distribution 
functions appearing in the non-homogeneous free propagators are replaced with equilib- 
rium Bose-Einstein distributions by the correspondence in (5.8). With this replacement, 
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the pinch singularities cancel by virtue of the KMS relation, as described at the beginning 
of this section. 

With the absence of pinch singularities established for early and infinitely-late times, 
we must now illustrate that the perturbation series remains free of these pathologies for all 
intermediate times. For this to be the case, we must demonstrate that the thermalization 
of the system takes place more rapidly than the growth of the pinch singularities discussed 
above. 

Potential pinch singularities arise from terms like 

gH'iix) , (6.16) 

where the argument x and the coupling constant g have mass dimensions. For large t, this 
term is dominated by the contribution from x — )• 0. By I'Hopitals rule and (3.51), we may 
show that these pathologies grow like 

g^^lix) -^^^ ^96{x). (6.17) 

t>i/Af 47r 

As a consequence, any order of pinch singularity will grow like {gt)^ , where n G N is 
some finite number. Therefore, in order for the perturbation series to be free of pinch 
singularities for all times, the deviation from thermodynamic equilibrium must decay more 
rapidly than a power law in the dimensionless parameter gt. 

In the spatially homogeneous limit, we may approximate the late-time behaviour of 
the statistical distribution function / by the following master equation: 

^ffl^ = -r>(p,t)/(|p|,t) + r<(p,t)(i + /(|p|,t)) , (6.18) 

where the four-momentum p'^ = (ii^(p), p) is on-shell. The master equation (6.18) is exactly 
the classical Boltzmann transport equation. The terms r>(p, t) = —iYly{p,t)/2M and 
r<(p, t) = — in<(p, t)/2M relate to the on-shell decay and production rates, respectively. If 
the system is close to equilibrium, which is reasonable for late but nonetheless finite times, 
we may assume that r>(<)(p, t) = r>(<)(p) are approximately constant. Expanding around 
the equilibrium Bose-Einstein distribution /b, we write 

/(|p|,t) = h{E{^)) + 5/(|p|,t). (6.19) 

Substituting this expansion into (6.18), the terms involving /b(-E(p)) cancel by virtue of 
the detailed balance condition (5.25) and the properties of the Bose-Einstein distribution. 
We then arrive at a first-order differential equation for the deviation (5i/(|p|,t): 

'^^^MA = _ (r>(p) - r<(p)) 5/(|p|,t) . (6.20) 

Notice that r>(p) — r<(p) is the total Breit-Wigner width r(p) ~ g"^, defined in (5.24b). 
Thus, 

<5/(|p|,t) = expf-r(p)tV/(|p|,0) , (6.21) 
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which decays faster than any power law in gt, as required. 

In summary, for early times, the analytic i-dependent vertices lead to microscopic 
violation of energy conservation. This energy non-conservation regularizes potential pinch 
singularities. For intermediate times, the t-dependent statistical distribution functions 
approach equilibrium faster than energy conservation is restored. In the limit t — >• c«, the 
time-dependent statistical distribution functions appearing in the non-homogeneous free 
propagators are replaced by their equilibrium forms via the correspondence in (5.8). We 
then obtain the well-known equilibrium thermal field theory in which energy conservation 
is fully restored and pinch singularities cancel exactly by virtue of the KMS relation. In 
conlusion, we have explicitly demonstrated how pinch singularities do not arise in our 
perturbative approach to non-equilibrium thermal field theory. 



7 The Particle Number Density 

It is important to establish a direct connection between off-shell Green's functions and 
physical observables. Such observables include the particle number density for which vari- 
ous interpretations have been reported in the literature [41-50]. In this section, we derive 
a physically meaningful definition of the particle number density in terms of the resummed 
CTP propagators. 

In order to count off-shell contributions systematically, we suggest to 'measure' the 
particle number density in terms of charges, rather than by quanta of energy. The latter 
approach would necessitate the use of a quasi-particle approximation to identify 'single- 
particle' energies, which we do not follow here. Instead, we begin by introducing the 
Noether charge Q{xo;ii) of the Heisenberg-picture scalar field ^}i{x;ii): 

Q{xo;ii) = ijd^x (^^l^{x;ii)7rl^ix;ii) - TTuix;ii)^Rix;ii)^ . (7.1) 

Here, we include all time-dependencies explicitly for clarity. We recall that 7r}i{x;ti) = 
dxo^li{x; ti) is the conjugate momentum operator. Up to the usual infinite T = vacuum 
contribution, the EEV of the charge operator Q{xq; ti) for a real scalar field <^ is zero, since 
particle and anti-particle contributions cancel. Therefore, we need to devise a method by 
which to separate the particle from the anti-particle degrees of freedom. 

For the general case of a spatially and temporally inhomogeneous background, we need 
therefore to generalize the Noether charge operator Q{xQ]ti) by writing it in terms of a 
charge density operator Q(p,X,Xo;fi) as 

Q(Xo;ti) = j d^X Q(p,X,Xo;t;) . (7.2) 

In the above, the three-momentum p is conjugate to the spatial part of the relative space- 
time coordinate = x^ — and X^^ = {X^, X) = {x^ + y^)/2 is the central space-time 
coordinate [cf. (3.2)]. To this end, we proceed by inserting into (7.1) unity in the following 
form: 

1 = y"d4y<5W(x -y) = jdS |-^e-P-(--y)5(xo - yo) ■ (7.3) 
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Observe that in what follows, = = thanks to 5{xq — yo). Subsequently, 
symmetrizing the integrand in x and y, we may write the charge operator as 

Q{Xo;ii) = '-jd^^ |d4y|-^e-*P-(--y)j(xo - yo) 

X ($^(y;fi)vrjj(x;ti) - TTu{y;ii)'^nix;ii) + (x ^ y)) . (7.4) 



In terms of the central and relative coordinates, X^ and R^, the charge density operator 
Q(p,X,Xo;ii) may be appropriately identified from (7.4) as 



Q{p,X,Xo;ii) = ^jd^Re-'^-^5{Ro) (7.5) 
X (^>|j(X-f;ti)vrli(X + f;ti) - ^n{X - f,ii)^Yi{X + f,ti) + (R ^ -R)^ 

Substituting for the definitions of the conjugate- momentum operators tth and vrjj, we may 
rewrite Q(p, X, Xq; ti) in the following form: 

Q(p,X,Xo;ti) = ijd^Re-'P'^6{R^) 

X dj^Ul{X - §;U)<^n{X + - $t^(X + f ; t,)^H(^ - f ; • (7.6) 



The EEV of Q(p,X, Xo;ti) at the macroscopic time t is then obtained by taking the 
trace with the density operator p as given in (4.3) in the equal-time limit Xq = tf. We 
have seen in Section 3 that the equal-time limit is necessary to ensure that the observable 
charge density is picture-independent and that the number of independent coordinates is 
reduced to seven as required. Thus, we have 



(Q(p,X,i7;ti))i = lim i /d^i? e'^P'^ 5(i?o) 



Xo^tf 

X dR,(^iA<{R,X,if;ii) - iA<{-R,X,if,k)'^ , (7.7) 

where we use the notation 

iA^{R,X,if;ti) = {^\^{X - f,ii)^H{X + (7.8) 

for the resummed CTP Wightman propagator. 

Let us comment on the two terms zA<(i2, X, tf; ti) and iA<(— i?, X, tj; ti) that occur on 
the RHS of (7.7). The first term iA^{R,X,tf;ti) comprises ensemble positive-frequency 
particle modes and ensemble plus vacuum negative-frequency anti-particle modes. The 
second one iA^{—R,X,tf;ti) comprises ensemble plus vacuum positive-frequency anti- 
particle modes and ensemble negative-frequency particle modes. Hence, we may extract 
the number density of particles by taking the sum of the positive-frequency contribution 
from iA^{R,X,tf;ti) and the negative-frequency contribution from iA^{—R,X,tf;ti). 

We separate the positive- and negative-frequency parts of (7.7) by decomposing the 
equal-time delta function 5{Ro) via the limit representation 
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with e = 0^. Thus, a physicahy meaningful definition of the number density of particles 
at the macroscopic time t is given by 



nip,X,if-ii) = - lim /d^Re^'P-^/ 

Xo ^if J J 



Xo ^if J J 2tt 

X ( J. \ . dR,i/\<{R,X,tf-ti) + ^ . dR,i^<{-R,X,tf-ti)] . (7.10) 

Using time translational invariance of the CTP contour, this observable may be recast in 
terms of the Wigner transform of the Wightman propagators as 

n(p,X,t) = n(p,X,t;0) 

= hm [^po(9{po)iA^{p,X,t;0) - e{-po)iA^{-p,X,t;0)) . (7.11) 

Xo t J ZTT \ J 

Note that the number density of anti-particles n'-'(p, X, f) is obtained by C-conjugating 
the two negative-frequency Wightman propagators in (7.11). Useful relations between 
correlation functions and their C-conjugated counterparts are given in Appendices A and C. 

We reiterate from (4.56) that the n(p, X, t) d^p d'^X is interpreted as the number of 
particles at macroscopic time t in the volume of phase-space between p and p + d^p and 
X and X + d'^X. The particle number per unit volume is obtained by integrating over all 
momentum modes, i.e. 

n(X,t) = y^n(p,X,t) (7.12) 
and finally the total particle number, by integrating over all space, i.e. 

X{i) = j d3Xn(X,t) . (7.13) 

By inserting the inverse Wigner transform 

zA<(p,X,t;0) = |^e-^^-^iA<(p+f,p-f,t;0), (7.14) 

into (7.11), the particle number per unit volume n(X, i) may be expressed in terms of the 
double momentum representation of the Wightman propagators via (7.12). After making 
the coordinate transformation p — )■ —p in the negative-frequency contribution, we then 
obtain 

n(X,t) = hm 2//^-^e-^^-^%o)wA<(p+f,p-f,i;0), (7.15) 



Xo^t J J (27r)4 {2'kY 

The particle number per unit volume ra(X, t) is also related to the statistical distribution 
function /, through 
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using (4.56). Equation (7.16) relies on our earlier observation that the statistical distri- 
bution function / appearing in the non-homogeneous free propagators of Section 4.3 is a 
dressed object, comprising all possible multi-particle background contributions. 

It is instructive to check that our definitions for the number density lead to the expected 
results for the free and quasi-particle equilibrium cases. Substituting the free equilibrium 
Wightman propagator (5.10c) of the real scalar field into (7.11), we obtain 

neq(p,X,i) = fB{E{p)) , (7.17) 

exactly as we would expect for the number density of free particles in thermodynamic 
equilibrium. Inserting instead the resummed equilibrium Wightman propagator given in 
(5.23) in the narrow width limit, we get 

nqp(p,X,t) = /b(^(p)) , (7.18) 

where £{p) is the solution to the gap equation, 

^^'(P) = p2 + m2 - RenR(f(p) +ie,p) , (7.19) 

and nqp(p,X,t) then represents the number density of quasi-particles. 



8 Master Time Evolution Equations for Particle Number Densities 

Having established a direct relationship between the non-homogeneous CTP propagators 
and the particle number density in Section 7, we are now in a position to derive in this 
section master time evolution equations for the particle number density n(X, t) and the 
statistical distribution function /. This is achieved in analogy to the derivation of the 
well-known Kadanoff-Baym equations [52, 53] by partially inverting the CTP Schwinger- 
Dyson equation obtained in Section 4.1. Our approach, however, differs significantly from 
other methods. As discussed in Section 4.2, we do not rely on a truncation of a gradient 
expansion of the resulting expressions. More details of the gradient expansion are given 
in Appendix E. In the next section, we will employ a loopwise truncation of the time 
evolution equations in terms of non-homogeneous free CTP propagators. As we will see, 
these dynamical equations are nonetheless resummed to all orders in a gradient expansion. 

We begin our derivation of the time evolution equations with the double momentum 
representation of the Schwinger-Dyson equation in (4.32). We convolute (4.32) consecu- 
tively from the right with the weight function (27r)'^(5j^^^(gi — q2) defined in (3.49) and 
then with the resummed CTP propagator A°'^{q2,p2, By making use of the inverse 

relation (3.48), we obtain the following expression: 

A0r'(pi,'?i)(2vr)%(')(gi - g2)A^,(g2,P2, t/; t.) = r^ai>{2^)Ht\pi - P2) 
dV ^..^^^o^^4.{4) 



(27r)4 {2tiY 



Uacipi,qi,tf;ti){27r) 61 ' {qi - q2)A%{q2,P2,tf;ti 
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where the contribution of the bi-local source Kab is neglected next to the self-energies. It is 
essential to remark that the LHS of (8.1) has the following coordinate-space representation: 



/ d^zAl'-\x,z)A''^{z,y,ij;U). (8.2) 



Substituting for the free inverse CTP propagator A[|^ ^{x,y) given in (3.21), we may then 
confirm via (8.2) that evaluating the qi and q2 integrals on the LHS of (8.1) yields 



// 



A0-l(pi,gi)(2vr)%(^)(<Zl - q2)A%{q2,P2,if,U) 

' p\ - M^)Aab{pi,P2,if;ii) ■ (8.3) 



Recalling that the self-energy Ilabip,p', if',ii) contains 5\^\p — p') functions in the vertices, 
we may perform the qi integral on the RHS of (8.1) by making use of (3.52). Consequently, 
(8.1) may be written down in the following concise form: 

(p2 _ M^'^Aab{pi,P2,if;ii) = Vab{2TT)^Si''\pi - P2) 

Uac{pi,q, if, ii)A%{q,p2, if, ii) . (8.4) 

With the definition of the particle number per unit volume n(X, t) from (7.12) in 
mind, we equate the element (a, h) = (1, 2) of each side of (8.4) to extract the interacting 
Klein-Gordon equation of the negative-frequency resummed Wightman propagator: 



pi - )A<{pi,p2,tfti) 

= (y^{Puq,if,ii)A<{q,p2,if,ii) - Il<{pi,q,if,ii)Aj){q,p2,if,ii)^ . (8.5) 

Using the decomposition of Aj:){p,p' ,if,ii) from (A.6e) and an analogous identity for 
Ii{p,p' ,if,ii) analagous, we may rewrite (8.5) as 

\ _ , /■ d^q 

pI - M'^jA<{pi,p2,if,ii) + jj^^Ti'p{pi,q,if,ii)A<{q,p2,if,ii) 

1 /■ d^g r ~ ~ ~ ~ 

= ~2Jj2tjY [^>*^^i'^'*^'**)^<*^^'^2>i/;*i) 

- U^{pi,q,if,ii) ( Ay{q,p2,if,ii) - 2A-p{q,p2,if,ii)) , (8.6) 



where the subscript V denotes principal-part evaluation of the functions given in (2.29) 
and (4.27b). 

Introducing the central and relative momenta, p = {pi + P2)/2 and P = Pi — P2, 
respectively, we write (8.6) in the following form: 

(po + f )' - E\p + f)]A<(p+f,p-f,t>;t,) 

+ ^{p+^,p- ^,ifii) = '^{p+^,p-^,if,ii) , (8.7) 



-47- 



where we have defined the following: 

d^g 



(27r)4 ^^'P^P^ i,Q,if;ii) iA<(g,p- ^,if]ii), 



in>(p+ ^,q,if;ii) iA^{q,p- ^,if;ii) 



iU<{p+ ^,q,tf;ti) [iAy{q,p- ^,tf;ti) - 2iA-p{q,p - ^ ,tf;ti 



(8.8b) 



With the aim of finding a master time evolution equation for the particle number 
density n(X, t) in (7.12), we integrate both sides of (8.7) with the measure 

d^p d^P ^_,p.^ 



Explicitly, this gives 



d^p d^P 



-iP-X 



(27r)4 (27r)4 



e{po) {PO + f)' - i^'(p+f) A<(p+f,p-f,t7;t,) 



0{po) . 



.9) 



(27r)4 (27r)4 



+ 



dV d^P 



(27r)4 (27r)4 

d^p d'P _,p.^ 



e{po)^{p+^,p-^,if-ii) 



(27r)4 (2vr)4 ^ ^(Po) ^(p + |,P - |, t/; ti) • (8.10) 

Adding to (8.10) the complex conjugate of the same expression with P — t- —P and using 
the identities in (A. 4), we may extract the terms proportional to p - P on the LHS of (8.10). 
In this way, we obtain 

d^p d'P 



+ 



(2^)4 (27r)4 

d^P dip ,p.^ 
(27r)4 (27r)4 

dV d^P „-.p.X 



0(po)p--PA<(p+ f ,p- f ,t/;ti) 



^.11 



^(po)( ^(p+f,p-f,t/;ti) + =^*(p-f,p+f,t/;t~) 



^(Po)('^(p+f,p-f,i/;ti) + ^*(p-f,p+f,t/;ti) 



(2^)4 (27r)4 

The first term on the LHS of (8.11) may be rewritten as 
d^p d^P 



e{pQ)[ ipodxo - P • P )e-'^-^ A<(p+ f ,p 



p 

2 



2 ) tf'^ ti) . 



i.l2) 



(27r)4 (27r)4 

Using time translational invariance of the CTP contour and taking the limit Xq — s- t 
in (8.12), we recognise that the derivative term with respect to t is precisely the time- 
derivative of the particle number density n(X, t) in (7.12). Hence, from (8.11) and (8.12), 
we arrive at our master time evolution equation for n(X, t): 

d^p d^P ^_,p.^ 



+ 



(27r)4 (27r)4 

d^p d^P ^_,p.x 



(27r)4 (27r)4 

d^d^ _iP. 
(27r)4 {2TTy " 



p-P%o)A<(p+f,p-f,t;0) 



0{P^){ ^(p+|,p-|,t;0) + ^*(p-|,p+|,i;0) 



^.13) 



■^e(po)('^(p+f,p-f,t;0) + <^*(p-f,p+f,t;0)) , 
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with Xq = t, where ^ and ^ are defined in (8.8a) and (8.8b). Comparing with the 
full non-truncated form of the Kadanoff-Baym kinetic equation in (E.4b), the series of 
nested Poisson brackets in (E.4b) has been replaced by a single convolution integral over 
the central momentum P in (8.13). 

In addition to the master time evolution equation for the particle number density 
n(X,t) in (8.13), we may respectively find a time evolution equation for the statistical 
distribution function /(p + — ^,t). Specifically, given the relation (7.16) and the 
differential equation (8.13), the following time evolution equation may be derived for /: 

9i/(p+ f ,P - f - //^ ^ e-^o*p- P0(po)A<(p + f f ,t;0) (8.14) 
+ |y'^^e-^o*0(po)(^(p+f,p-f,t;O) + ^*(p-f,p+f,t;0)) 

= //^^e-^"*0(po)(<^(p+|,p-f,t;O) + <r*(p-f,p+f,t;0)). 

It is important to stress here that (8.14) provides a self-consistent time evolution equation 
for / valid to all orders in perturbation theory and to all orders in gradient expansion. 
Nevertheless, it is a formidable task to provide an evaluation of (8.14) to all orders in 
perturbation theory. Therefore, in the next section we develop a systematic perturbative 
expansion to non-equilibrium dynamics based on the result (8.14) to all orders in gradients. 

We note that a physical intepretation may be attributed to the different terms con- 
tributing to (8.14). Specifically, all terms on the LHS of (8.14) may be associated with the 
total derivative in the phase-space (X, p), appearing in the classical Boltzmann transport 
equation [75]: 

Dt = 9t + vVx + F-Vp , (8.15) 

where v is the average non-relativistic velocity of the particle distribution and F is the 
force acting on this distribution. In particular, the ^ terms on the LHS of (8.14) are the 
force terms, generated by the potential due to the dispersive part of the self-energy. On 
the RHS of (8.14), the ^ terms represent the collision terms. 

It would be interesting to discuss the spatially homogeneous limit of (8.14) at late 
times. In this case, energy conservation holds to a good approximation, so the time- 
dependent weight functions 6t may be replaced by the standard Dirac delta functions, even 
in the vertices of the self-energies contained in the force ^ and collision terms of (8.14). 
Moreover, as a consequence of the assumed spatial homogeneity, all propagators and self- 
energies, which in general depend on two momenta pi and p2, will now be proportional to 
the four-dimensional delta function {2tt)'^5^'^\pi — P2)- Likewise, the statistical distribution 
function / takes on the form 

/(p+f,p-f,t) = (27r)V3)(P)/(|p|,t). (8.16) 

Because of the above simplifications, one can then work in the single momentum repre- 
sentation, by integrating over the three-momentum P. Thus, we find the following time 
evolution equation for f{\p\,t): 

dtfi\p\,t) = l'^9{po)[iUy{p,t)tA^ip,t) - iU^{p,t)iA^{p,t)) , (8.17) 
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where the purely-imaginary force term ^ and off-shell effects from Ap have vanished for 
= (0,0). This result corresponds to the semi-classical Boltzmann transport equa- 
tion in (6.18) or, equivalently, to the zeroth-order truncation of the gradient-expanded 
Kadanoff-Baym kinetic equation in (E.5b) with Xq = t. 



9 Thermalization in a Scalar Model 

We now apply the formalism developed in the proceeding sections to a simple scalar model. 
In particular, we introduce the modified Feynman rules that result from the systematic 
inclusion of finite-time effects and the violation of both energy conservation and space- 
time translational invariance. As a playground for studying the kinematics in the early- 
time energy-non-conserving regime, we calculate the time-dependent thermal width of a 
heavy scalar We show that processes, which would normally be kinematically disallowed, 
contribute significantly to the prompt shock regime of the initial evolution. We also show 
that the subsequent dynamics exhibit non-Mar kovian behaviour, acquiring oscillations with 
time-dependent frequencies. This evolution signifies the occurence of memory effects, as 
are expected in truly out-of-equilibrium systems. Finally, we look in more detail at the time 
evolution equations of this simple model and demonstrate the importance of the violation 
of energy conservation to the statistical dynamics. 

We consider a simple scalar theory, which comprises one heavy real scalar field ^ and 
one light pair of complex scalar fields {x^ , x), described by the Lagrangian 

C[x) = ldi,^{x)d''^{x) - lM^<^^{x) + d^,xHx)d''x{x) - m^x\x)x{x) 

- gHx)x\x)x{x) - lX[xHx)xix)]'^ , (9.1) 

where M ^ m. Appendix C describes the generalization of our approach to the complex 
scalar field x- 

We formulate a perturbative approach based upon the following modified Feynman 
rules: 

• sum over all topologically distinct diagrams at a given order in perturbation theory. 

• assign to each ^-propagator line a factor of 

p ;; p' = iAl''\p,p',ir,u). 

am ► -*- — ► - mb 

• assign to each x-propagator line a factor of 

P \\ p' = iA'^-'{p,p',if;U). 

a« ». — It — mb 
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• assign to each three-point $ vertex a factor of 



b 




X = -^gr]abc i2n)\ ( ^ Li Po,) S'^^') ( ELi pO , 



c 



where 5t is defined in (3.50) and the pre-factors riabc... are given after (3.12). 

• assign to each four-point x vertex a factor of 
a b 

X = -iXVabcd ( Ei=l PO,) S^'^ ( Ei=l Pi) ■ 

d c 

• associate with each external vertex a phase 

where po is the energy flowing into the vertex. 

• contract all internal CTP indices. 

• integrate with the measure 

J (27r)4 

over the four- momentum associated with each contracted pair of CTP indices. 

• consider the combinatorial symmetry factors, where appropriate. 

Notice that there are a number of modifications with respect to the standard Feynman 
rules. In particular, the familiar energy-conserving delta function has been replaced by 
5t in the vertices. This is indicated diagrammatically by the dotted line terminated in a 
cross, representing the violation of energy conservation. This loss of energy conservation 
is a consequence of Heisenberg's uncertainty principle, due to the finite macroscopic time 
of observation t, over which the interactions have been switched on. The time-dependent 
vertices vanish in the limit t — )■ 0, as we should expect. The loss of space-time translational 
invariance leads to a doubling of the number of integrations with respect to the zero- 
temperature and equilibrium cases. 
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Figure 4: One-loop local x self-energy: il^^^'^^ {q,q\if;ti) oc rjah- 

At the one-loop level, we have three diagrams. The local x self-energy shown in 
Figure 4: 



X {2T:f5[''\q - q' - k + k')7]abcdiA^^'''{k, k' ,i f,U) , (9.2) 
and the two non-local diagrams for the ^ and x self-energies shown in Figure 5: 

ill^^,,{q,q,tf,t,) - 2, [ZTTfi) e J j (27,)^ {2t,Y {2t,Y {2T,f 

X {2Tif 5\^\q - ki - k2){2T:f5\^\q' - k[ - k'^) 
X 7]acdi^^^''^{ki,k[,if;U)i^^'°''^^{k2,k'2,if-ii)7]efb , (9.3a) 

^iix,a6W'9'«/'*0 - 2! e y y (27r)'^ (27r)4 (27r)M2vr)4 

X (27r)^4^)(g - fei - A:2)(27r)^4^)(g' - k[ - k'^) 
X riacdi^'^^^'iki, k[,if; ii)i^Y^ {k2, k'^, if;ii)r]efb ■ (9.3b) 

A detailed description of the techniques required to perform these loop integrals is provided 
in Appendix F. 

9.1 Time-Dependent Width 

In this section, we study the time-dependent width of the heavy scalar <1>. In particular, we 
investigate the spectral evolution that results from the restoration of energy conservation, 
without solving explicitly the system of evolution equations for the statistical dynamics. 

We consider the following situation. We prepare two isolated but coincident subsys- 
tems and both separately in thermodynamic equilibrium and at the same temper- 
ature r, with all interactions turned off. The subsystem contains only the real scalar 
field <I>, whilst Jy^^ contains only the complex scalar x- At macroscopic time t = 0, we turn 
on the interactions and allow the system ^ = U to re-thermalize. For our numerical 
analysis, we take for definiteness the thermodynamic temperature to be T = 10 GeV, the 
mass of the heavy $-scalar M = 1 GeV, the mass of the complex x-scalar m = 0.01 GeV, 
and their trilinear coupling g = Q.l GeV. 
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(a) ini^,U(9> <?', i/; (b) 

Figure 5: Non-local one-loop $ and x self-energies. 

The free propagators of the fields $ and x a-t time t = are the equilibrium propagators 
in (5.10) and (C.25), containing the Bose-Einstein distributions at temperature T. We take 
the chemical potential of the complex scalar to be vanishingly small in comparison to its 
mass, i.e. // <C m, such that /;^(|p|,0) = /^(|p|,0) = /b(-Ex(p))- Without solving the 
system of time evolution equations, the form of the the statistical distribution functions of 
the <I> and x scalars is unknown for t ^ Q. We assume that the heat bath of x's is sufficiently 
large so as to remain unperturbed by the addition of the real scalar <I>. Specifically, we 
may consider the number density of x's to remain unchanged and the free equilibrium x 
propagators in (C.25) to persist for all times. 

By the optical theorem, the width r$ of the scalar <I> is defined in terms of the absorp- 
tive part of the retarded $ self-energy Im r via 



where 



r*(gi,92,i/;ti) = -^imn$,R(gi,g2,t/;ii) > (9-4) 

Imn$,R(gi,g2,t/;ii) = ^(n$,>(gi,g2,i/;?i) - n$,<(gi,g2,i/;ti)) • (9.5) 

At the one-loop level, the self-energies ^^^^^{qi, q2,i f]ii) are given by (9.3a). 

Employing the relative and central momenta, Q = qi — q2 and q = {qi + '?2)/2, and 
using the results of Appendix F, the Laplace transform with respect to the macroscopic 
time t of the one-loop ^ width is 

r'i\q + Q/2,q-Q/2,s) = {27r)'l ^^^e'^^^'f S'^'HQ) 



32vr2M TT E,E2 (<^2 _ ^^^^ - a^E^f + ' 
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where we use the short-hand notation {a} for the summation over ai, 0.2 = ±1- For 
the sake of generality, we distinguish the x"*" and x~ decay products by assigning them 
different masses mi and m2, respectively, so that 



El = Ei{k) = J|k|2 + ml , (9.7a) 



E2 = ^2(q-k) = y|k|2 - 2|k||q|cos6' + |q|2 + m^ . (9.7b) 

In our numerical analysis, however, we take mi = m2 = m. 

Performing the inverse Wigner transform with respect to Q of (9.6) in the equal-time 
limit Xq = tf and subsequently the inverse Laplace transform with respect to s, we obtain 
the time-dependent $ width 

{a} 

1 + /B(ai^i) + /B(a2^2)) . (9.8) 



In the limit t — )• 00, the sine function on the RHS of (9.8) yields the standard energy- 
conserving delta function: 

lim -sincr(5o - ai-Ei - 02-^2)*] = ^iqo - aiEi - a2-E'2) , (9.9) 
t— s>oo7r 

thereby recovering the known equilibrium result. 

Figure 6 contains a series of contour plots of the differential one-loop $ width d^F^"* 
evaluated over the dominant region of the \k\-9 phase space, where the four-momentum q^^ 
of the $ scalar is on-shell, i.e. q^ = M^. For late times, the integrand is highly oscillatory 
and we expect the dominant peak of the sine function to approach the region of phase 
space permitted in the limit f — t- 00 [cf. (9.9)], for which 



M2|q|cos(9 + y(|q|2 + M2) [M^ - 4m2(M2 + |q|2sin26')] 
" 2(M2 + |q|2sin2e) ' ^^'^^^ 

For early times, the admissible phase space is greatly expanded. For later times, the 
frequency of oscillation increases and the width of the central peak of the sine function 
narrows, lying along a curve in the phase space. To proceed further, we need to develop a 
method for dealing with the kinematics in the absence of exact energy conservation. 

9.2 Generalized Two-Body Decay Kinematics 

At zero temperature and density, it is convenient to analyze the two-body decay kinematics 
by performing a Lorentz boost to the rest frame of the decaying particle. However, at finite 
temperature, the dependence on the thermodynamic temperature of the heat bath breaks 
the Lorentz covariance of the integral. As such, we cannot eliminate dependence upon the 
the three- momentum of the decaying particle by any such Lorentz boost; the dependence 
will reappear in the 'boosted' temperature: 

T' = 'jT , (9.11) 
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Figure 6: Contour plots of |k| versus 6 for d^r|^^ for discrete values of Mt, with = 
and |q| = 10 GeV. The region of phase-space permitted in the t — )• oo limit is shown by 
the blue dashed line, corresponding to the delta function in (9.9). 



where 

7 = (^1 + ^ j (9.12) 

is the usual Lorentz boost factor for the heavy scalar field $. As a result, we are compelled 
to analyze the kinematics of the two-body decay in the rest frame of the heat bath, which 
we define to be the frame in which the EEV of the three- momentum operator (P) is 
minimized. For an isotropic heat bath, this is the frame in which (P) = 0, that is the 
comoving frame. In this section, we look more closely at the kinematics in the absence of 
energy conservation. 

For this purpose, let us introduce the variable 



u 



{qo - aiEi - a2E2)t , (9.13) 
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which may be interpreted in terms of energy borrowed from or lent to the heat bath. We 
shall hereafter refer to the variable u as the evanescent action of the process, since it 
has the correct dimensions and quantifies the extended kinematically allowed phase-space 
configurations. We also define the evanescent energy 



Quit) = qo 



u 



(9.14) 



which satisfies 



lim Quit) = qo . 



(9.15) 



With this substitution, we obtain the kinematic constraint 



quit) - aiEi - a2-E'2 = . 



(9.16) 



Since u can take large positive values, quit) is not necessarily restricted to positive values 
for early times, even when qo = y^|qp + M"^ > is on-shell. Processes with u/t ^ 
are referred to as evanescent. 

In order to make the coordinate transformation |k| — t- li in (9.8), we must solve (9.16) 
for the magnitude of the three-momentum |k(t)| , which becomes implicitly time-dependent. 
Specifically, we find 



|kW(i) 



1 



Hilit) 



|q|^ COS' 



{Qlit) 



I |2 , 2 

|q| + rrii 



rrio 



|q| cos 6 



+ bquit)J Xiqlit) - |q|2,mf,m|) - 'irnl\c(\'^ slv? 



(9.17) 



with 6 = ±1 and 

Kx,y,z) = ix-y-z)'^ - Ayz . 
After a little algebra, we obtain the energies 



E?it) = Ql 



1 



2(g2(t) _ |q|2cos2^) 



llii) - |qP + ml - ml)quit) 



Ei'\t) = 02 



+ 6 |q| cos^Y A(g2(*) — |qpi"if,m|) — imflql"^ sin^O 
1 



(9.18) 



2(g2(t) _ |q|2cos2e)L 



iqlit) - |q|2cos2e - + m^)^„(t) 



b \q\ cos 9 ^1 Xi^q'^it) — |q|2, m|,m|) — 4r?T,2|q|2 sin^^ 



(9.19) 



(9.20) 



where the overall factors of ai and a2 are necessary to satisfy the initial constraint (9.16). 
It is clear that these results collapse to the kinematics of equilibrium field theory in the 
limit u/t — > 0. 
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ai = 02 = +1 


Ol = —02 


Oi = 02 = —1 


u+{t) 


— oo 


a;o(g,oi,02)t 



+00 
CJo(g,Oi,02)t 



Table 3: Limits of integration of the evanescent action u for each of the four processes 
contributing to the non-equihbrium thermal $ width for t > 0, where the angular frequency 
LOo{q, ai, 02) is defined in (9.23). 

Keeping mi and m2 distinct, the $ width for t > is given by 

r«r V- r^'l, ■ sing 

L [q,t) = ^ ^ „ , ^ > dU dn smcfnl ^ 

^^^^^Jo Lit) [qlit) - H^cos^ef 

\{ql{t) - |qp, m?,m|) - 4m?|qp sin^g^ ^''^ 
i^ui^) ~ IqI^ + "^1 ~ m2)|q|cosg 



+ bQuit)y X{ql{t) — |qp, mfjmg) — 4mf|qPsin^6' 

X 0102(1 + /B(oi^f)(t)) + /B(o2£;f (i))) . (9.21) 

The reality of the loop momentum |k(t)| in (9.17) requires that the discriminant 

\{qlit) - |qp,m?,m^) - 4m?|q|^ sin26' > 0, (9.22) 

which is now a time- dependent kinematic constraint. Furthermore, we require |k*^^^(t)| > 0, 
Ei'\t) > nil and E^\t) > 1712- For t > 0, the limits of integration u±{t) are given in 
Table 3, where we have defined 

/ |q|2 

a;o(g, 01,02) = qo - (oimi + 02771,2) 1 + -2 , (9.23) 

\ (mi + m2) / 

which is the angular frequency of the sine-integral-like oscillations of the integral. For 
t = 0, u±{t) = and the domain of integration over u collapses to zero. Given the 
analytic behaviour of the integrand in the limit t — >• 0, the integral vanishes, as we expect. 

For on-shell decay modes with q^ = M^, the angular frequency iOo{q, oi, 02) in (9.23) 
becomes 



wo(|q|, 01,02) = wo(g,ai, "2)1^2 = ^2 = \/|qP + Af2 - y |q|2 + (mi + 1112) . (9.24) 

Thus, the evolution of the phase space for on-threshold decays with = (mi -|- m2)^ is 
critically damped. We note that in the large momentum limit |q| » M, 

a;o(|q|,ai,a2)||q|»j^,^ ^ ^1^^ -, (9.25) 
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Figure 7: The dependence of a;o(|q|,ai,a2)/^ on |q|/M for processes related to 1 — )• 2 
decay (solid black), Landau damping (blue dashed) and 3 — )• total annihilation (red 
dotted), where = and mi = m2 = m = 0.01 GeV. The latter two of these 
processes become highly oscillatory for large |q|/M. 




Figure 8: Processes contributing to the time-dependent $ width: (a) the familiar 1 — )• 2 
decay; (b) 2 — )• 1 Landau damping; and (c) 3 — t- total annihilation. 



such that the evolution of the phase space for high-momentum modes is similarly damped. 
The momentum dependence of a;o(|q|, ai, 02) is shown in Figure 7. 

The summation over ai and Q2 yields four distinct contributions to the decay width. 
For ai = a2 = +1, we obtain the contribution from the familiar 1—7-2 decay process 
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Figure 9: Contour plots of u versus 9 for the 1 — )• 2 decay contribution to d^F^j, for 
discrete values of Mt, with q'^ = and |q| = 10 GeV. The solid excluded regions to 
the right of the red dotted line lie exterior to the limits of integration over u. 



presented in Figure 8a. For ai = —02 = ±1, we obtain the two 2 — )• 1 Landau-damping 
contributions displayed in Figure 8b. For ai = 02 = —1, we obtain the 3 — )• total 
annihilation process shown in Figure 8c. In the latter, the decay 'products' appear in the 
initial state along with the decaying particle. For late times, the Landau damping and 
total annihilation processes are kinematically disallowed, as we would expect. They are 
only permitted in the evanescent regime at early times. 

Figures 9-10 contain contour plots of u versus 9 for the four contributions to the on- 
shell differential width d^F^^. The equilibrium kinematics are obtained in the large-time 
limit t — )• 00 as follows. As can be seen from the (u, 9) contour plots in Figure 9, for the 
1 — )• 2 decay process, the limits of integration grow to encompass the full range of the 
sine function. At the same time, the u dependence of the phase-space pre-factors vanishes, 
since Quit) — ?• qo as t — )• 00. For the 3 — )• total annihilation process, the domain of 
integration vanishes in the large-time limit t — )• cxd. Given that the integrand is finite 
in the same limit, the contribution therefore vanishes as expected, which is confirmed by 
our {u, 9) contour plots in Figure 10. For the two Landau-damping contributions, the 
large-time behaviour becomes more subtle. As t — )• 00, the domain of integration covers 
approximately all positive u. However, the kinematically allowed phase space cannot be 
attained with any values of u and 9, so that the two Landau-damping contributions also 
vanish in the large-time limit. This behaviour is reflected in our {u,9) contour plots of 
Figure 11. Thus, only the usual 1 — >• 2 energy-conserving decay remains for late times. 

In order to reduce the statistical error in our Monte Carlo integration over the (u, 9) 
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Figure 10: Contour plots of u versus 9 for the 3—7-0 total annihilation contribution to 
d^r^^ for discrete values of Mt, where = and |q| = 10 GeV. The regions to the 
left of the red dotted line are exterior to the domain of integration over u. 



phase space, we use a Gaussian sampling bias to ensure that the majority of sampling 
points fall over the dominant region of the sine function of u in (9.21). We define the 
weight function 

du f iu — no)^ 

d^^ " 



W[U) 



(9.26) 



where 



Erf 



r{u) 



1 U — UQ 

s/2 o-u 



■ "0 



Erf 



1 M+ - «0 

\p2 o-u 



Erf 



1 - mq 



e [0, 1] . 



(9.27) 



After performing the change of the variable n — )• r in (9.21) for uq = 0, the limits of 
integration become time-independent and are identical for the decay, total-annihilation and 
Landau-damping contributions. The dependence upon the limits it+(i) and u^{t) appears 
instead within the transformed integrand of (9.21) in the new variables (r, 6). 

The width of the dominant region of the u phase-space is taken to be the distance 
between the central maximum and the maximum at which the amplitude of the sine func- 
tion has fallen to 0.1% with respect to the central maximum. Hence, we require for the 
distant maximum to be at the location where sincu„ ~ 0.001. The extrema m„ of the sine 
function satisfy the transcendental equation 



tan?i„ , 



(9.28) 
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Figure 11: Contour plots of u versus 6 for the two Landau-damping contributions 
(ai = —a2 = ±1) to d^r|^^'' for discrete values of Mt, assuming = and 
|q| = 10 GeV. The regions to the right of the red dotted line are exterior to the do- 
main of integration over u. The contour plots differ between the two contributions due to 
the asymmetry of the integrands in the ^ three-momentum. 
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Figure 12: The four separate contributions to the ratio in (9.30) versus Mt, for 
on-shell decays with |q| = 1 GeV, 10 GeV and 100 GeV. The two Landau-damping 
contributions are identical up to numerical errors. 



whose solutions for n > 1 may be expressed as [76] 



I 1 

Un = ^nir exp i d^ - arg 



1 + l^in^^ ^ ± ^vrie) + nVe' 

2^ 1 + ^ 2 ^ ' ^ 



(9.29) 



The required extremum is then given by n = 318, the 159*^ maximum of the sine function, 
with U31S = 1000. The Gaussian weight function ^{u) in (9.26) is then taken to have a 
variance of cr^ = (u3i8/2)'^, such that 95% of the sampling points fall within this dominant 
region. 
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Figure 13: The ratio F^^^ in (9.30) versus Mt, for on-shell decays with |q| = 1 GeV (sohd 
black), 10 GeV (blue dotted) and 100 GeV (red dashed). 



Our interest is in the deviation of the one-loop $ width r|^^'* from the known equilibrium 
result. It is therefore convenient to define the ratio f ^^(|q|, t) of the time-dependent width 
to its late-time equilibrium value for = M 



2. 



fS^(|q|,i) 



rS^(|q|,i) 



ri'^'"'(|q|,t) + 2ri^)^^°^^"( 



|q|,t)+rp-^°(|q| 



oo 



ri')(|q|,t^oo) 



(9.30) 

In Figure 12, we plot separately the four contributions from Figure 8 to this ratio as a 
function of Mt for a series of discrete momenta. The evanescent Landau-damping processes 
yield a prompt contribution, which can be as high as 10-20% at early times. The evanescent 
total annihilation process contributes similarly at the level of about 5%. 

In Figure 13, we plot the total ratio r^^(|q|,t) as a function of Mt. The total $ 
width F^^(|q|,i) is vanishing for Mt = 0, as we would expect. This is followed promptly 
by a sharp rise, which is particularly pronounced in the infra-red modes, resulting from 
the sudden switching on of the interactions. This so-called shock regime is followed by 
the superposition transient oscillations of angular frequency cjq with short time scales and 
non-Markovian oscillations of longer time scales. The latter of these oscillations exhibit 
time-dependent frequencies, the origin of which will be discussed in the next section. 
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Figure 14: Non-Markovian oscillations of the 1—7-2 decay contribution to the ratio F^^^ 
in (9.30) against Mt, by performing a moving average over bins of 150Mt. 

9.3 Non-Markovian Oscillations 

As we have seen in Figures 12 and 13, the time-dependent ^ width contains a superposition 
of damped oscillatory contributions. The longer-lived of these oscillations exhibit time- 
dependent frequencies. These non-Markovian oscillations are illustrated more clearly in 
Figure 14 for the |q| = 1 GeV mode in which a moving time average is carried out to 
eliminate the higher-frequency Markovian oscillations. In this section, we describe the 
origin of the time-dependent oscillations and show that they are not a numerical artefact, 
but are instead an intrinsic feature inherent to the dynamics of truly out-of-equilibrium 
systems. To this end, we consider the high-temperature limit T ^ M of the time-dependent 
$ width. 

In this high-temperature limit T ^ M, the Bose-Einstein distribution may be ap- 
proximated as follows: 

fB{E) « I , (9.31) 

Returning to the <I> width F^^ in (9.8) and substituting for (9.31), we may perform the 
angular integration by making the following change of variables: 

2|k| |q| 
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Then, the high-temperature hmit T S> M of the time-dependent width of the heavy 
scalar $ becomes 

El 



■.(1),T> Af 



9 

32ir^M 



Ef°° Ikl 
Qiae / d|k| — 

{q} 



(9.33) 



1 + 



+ 



T 



aiEi qq — aiEi 



Si[(g'o - aiEi - a2E2)t\ 



qo - aiEi 



sin[(go — Qii?i)t] Ci(a2-E'2t) — cos[(go — aii?i)t] 81(02-^2*) 



where Si(x) and Ci(a;) are respectively the sine integral and cosine integral functions. We 
have introduced the short-hand notations: {a} for the summation over ai, 02, ae = ±1; 
and 



El 
E2 



k|2 - 2ae|k||q| + |q|2 + ra\ 



(9.34) 
(9.35) 



In terms of the evanescent action u given in (9.13), the high-temperature limit of the 
one- loop ^ width reads: 

.,2 fu+{t) 



r 



(i),r>M 



9 



327r2Mt 



E 



aia2 



{a},fe = ±l 



«-(*) 



Quit) 



+ 



mT — mr, 



mf + 7712 



mr — nin 



2|q| 



T 



ujf\q,u,t) Quit) - a;^"^(g,n,t) 



+ 



Si{u) 



T 



sm 



ujP{q,u,t)t) Ci{uj^''>{q,u,t)t) 



(b), 



ujf\q,u,t) 
cos(a;|^^(g, u, t)t) Si(c(;2^^(g, u, t)t) 



(9.36) 



where ijO^^\q,u,t) and UJ2"' {q,u,t) are time-dependent non-Markovian frequencies defined 



as 



io'l\q,u,t) = qo 



{Ql{t) - |qP + ml - mfjquit) + bae\q\X^/^{ql{t) - \q\^ ,ml,mf) 



2(5^(t) - |q|2) 



(9.37a) 



^2\Q,u,t) 



illit) 



|q|2 - mf + m'^)qu{t) - baeHX^^"^ {ql{t) - |q|2, mf,m|) 



Hilii) - lql 



(9.37b) 
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Figure 15: {q,u,t)/M versus Mt for the first three extrema of sinc(n): ui = 4.49 
(black), U2 = 7.73 (bkie) and us = 10.90 (red), obtained by (9.29), for = m'^, 
|q| = 10 GeV and ao = +1. Sohd lines correspond to b = 1 and dashed lines, b = —1 
in (9.37a). The dotted lines mark the upper limit of the kinematically disallowed region. 



The time dependence of the frequency uj\ {q, u, t) is shown in Figure 15 for ag = +1, with 
the on-shell constraint = M^. The b = 1 contribution persists for late times, at which 
point the frequency of the modulations have decayed to zero. The 6 = — 1 contribution is 
disallowed for late times, such that the amplitude of this constant-frequency modulation is 
damped to zero. For ag = —1, the b = +1 and b = —1 contributions are interchanged, 
such that the frequency of the 6 = — 1 contribution reduces to zero for late times. Thus, 
we obtain the expected kinematics and equilibrium behaviour in the late time limit. We 
have not plotted uj^\ as its behaviour is indistinguishable from lo^\ for rrii = m2 = m. 

Looking again at (9.36), we observe that these non-Markovian oscillations occur only 
for T 7^ 0. We conclude therefore that this behaviour signifies a genuine non-equilibrium 
statistical effect. 

9.4 Perturbative Time Evolution Equations 

In this section, we carry out a loopwise expansion of the time evolution equations derived 
in Section 8 to leading order in the coupling g, evaluating their one- loop structure for 
our simple scalar model. For perturbatively small couplings, such a loopwise expansion is 
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expected to accurately capture the early-time dynamics of non-equilibrium systems. This is 
a regime, in which the applicability of truncated gradient expansions becomes questionable, 
according to our discussion in Section 4.2. Given the closed analytic form of the free 
CTP propagators, both the amplitudes of contributing processes and the resulting phase- 
space integrals describing the kinematics can be analytically determined. The systematic 
treatment of these kinematic effects is essential to understand the consequences on the 
early-time dynamics. In particular, we illustrate the significance of contributions from the 
energy-non-conserving evanescent regime to the dynamics of the system. 

Let us first return to the Boltzmann-like equation in (8.17), in which we artificially 
imposed energy conservation. The collision terms for the real scalar <I> are 

n$,>(p,t)A$,<(p,t) - n$,<(p,t)A$,>(p,t) . (9.38) 

Since we have assumed that the two subsystems =5^$ and 5^-^ are separately in thermody- 
namic equilibrium at the same temperature T at the initial time t = 0, the <I> propagators 
and self-energies will initially satisfy the KMS relation, i.e. 

A$,>(p,0) = e^P«A$,<(p,0) , n$,>(p,0) = e^f°n$,<(p,0) . (9.39) 

As a consequence of the KMS relation, the collision terms in (9.38) are identical to zero. 
This is true also in the collision terms of the gradient-expanded Kadanoff-Baym kinetic 
equation in (E.5b). With no external sources present to perturb the combined system 
from this non-interacting equilibrium, we are faced with the problem that the statistical 
distribution functions will not evolve from their initial forms. In reality, we anticipate 
that the system should evolve to some new interacting thermodynamic equilibrium for late 
times. We conclude therefore that the energy-non-conserving evanescent regime described 
by our approach is entirely necessary to account correctly for the evolution of this system. 

We now turn our attention to the master time evolution equation (8.14). Truncating to 
leading order in the coupling g(, the RHS of this time evolution equation contains free prop- 
agators and one-loop self-energies. For instance, the one-loop collision terms in (8.14) are 

^«(p+f,p-f,t;0) ^ [i4'^(p,,g,t;0)iA0(g,p2,i;0) 

- inW(pi,(7,t) (iAO (g,p2,t;0) - 2iA^(g,p2,i;0))] . (9.40) 

We insert the one-loop integrals (9.3a) and (9.3b), containing the homogeneous $ and x 
distribution functions, /$(|p|, t) and "^(Ipl, t), into the full time evolution equation (8.14) 
for the $ and x fields. Note that the / distribution functions are real and the g distribution 
functions vanish in this spatially homogeneous limit. We tacitly assume that a system ini- 
tially prepared in a spatially homogeneous state remains spatially homogeneous throughout 
its evolution. This assumption is reasonable if the system has infinite spatial extent. 

We may perform the loop and convolution integrals using the techniques and results of 
Appendix F. After carrying out the po and Pq integrals in (8.14), we obtain the following 
one-loop time evolution equation for the spatially homogeneous statistical distribution 
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function /$(|p|, t) of the real scalar field: 

X ^ sinc[(a£'$(p) — aiii^j(.(k) — a2-Ej(^(p - k))t/2] 
X 1^ + 2Si[{aE^{p) + ai£;^(k) + a2Sx(P " k))i/2] } 
X {(^(-a) + /$(|p|,t))[0(ai)(l + /x(|k|,t)) + 0(-ai)/^(|k|,t) 
X e{a2)[l + 4^(|p-k|,t)) + 0(-a2)/x(|p-k|,t) 
0{a) + /<i,(|p|,t))[0(ai)/^(|k|,t) + 0(-ai)(l + /^(|k|,i) 

"^M/^(|p-k|,i) + 0(-a2)(l + /x(|p-k|,t))]} . 



where Si(x) is the sine integral function and {a} is the short-hand notation for the sum- 
mation over a, ai, 02 = ±1- With the summation over {a}, the statistical factors in 
the braces of the last four lines of (9.41) contain the difference of contributions from the 
four processes shown in Figure 8 and their inverse processes. For early times, all of these 
evanescent processes and inverse processes contribute, as can be seen from Figure 12. The 
presence of Si(x) ensures that the correct late-time limit and kinematics are obtained on 
restoration of energy conservation. This factor is missing in the existing descriptions, see 
e.g. [45, 77]. The dispersive force term and off-shell collision term vanish in the spatially 
homogeneous case, thanks to the symmetry of the self-energy under P — t- —P. 
In the large-time limit t — )• 00, we have 



2tt 



sinc[(a-E'(j)(p) — aiE^{k) — 02 E^{p — 1^)^1/2^ 

X jvr + 2Si[(a£;$(p) + aiE^{k) + a2E^{p - k))t/2] | 

2710 {a)6{E^{p) - aiE^{k) - a2E^{p-k)) 



t — > 00 



(9.42) 



The kinematic constraints then force a = ai 

Sk 



a2 



dthi\p\,t) 



9 



1 



-|-1 and we obtain 
1 1 



2 J (27r)3 2E^{p) 2E^{k) 2E^{p-k) 
X 27r5{E^ip) - E^{k) - E^{p-k)) 



X 



M\p\,t)[l + /x(|k|,t))(l + /^(|p-k|,t 



1 + M\plt) )fxi\Kt)f^{\p-k\,t) 



(9.43) 



Equation (9.43) corresponds to the semi-classical Boltzmann equation [cf. (8.17)]. 

By analogy, the one-loop time evolution equation for the spatially homogeneous sta- 
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tistical distribution function fxilpl^t) of the complex scalar field x is given by 
^2 d^k 1 1 1 



dtfxi\p\,t) 



{a} 
t 



(27r)3 2E$(p) 2Sx(k) 2^;,(p-k) 

. . 2^ sinc[(a£';^(p) - ai£^$(k) - a2-E'x(p - k))t/2] 
X + 2Si[(a-E^(p) + ai^$(k) + a2-E;^(p - k))t/2] } 
X {[^(a)/^(|p|,t) + e(-a)(l + /^(|p|,t))](0(ai) + /$(|k|,t) 



(9.44) 



e{a2)[ 1 + /x(|p-k|,t) ) + e(-a2)/^(|p-k|,t) 



^(a)( 1 + /^(IpM) ) + 0(-a)/5^(|pM) ( e{-a{) + /$(|p|,t) 



^(«2)/x(|p-k|,t) + e{-a2)[l + /^(|p-k|,t))]} . 
In the large-time limit t — t- oo, the kinematics restrict a = ai = —a2 = +1, giving 

-2 /"d^k 1 1 1 



r 

' 2 J (27r)3 2£;$(p) 2^^(k) 2£;^(p - k) 
X 2^5(^$(k) - E^iv) - ^x(P - k)) 



1 + /$(|k|,t) /^(|p|,t)/^(|p-k|,t) 



/<i,(|k|,t)(l + /^(|p|,t))(l + /^(|p-k|,t) 



(9.45) 



This result is consistent with (9.43), differing by an overall sign, as one should expect. 
At i = 0, the semi-classical Boltzmann transport equation in (9.43) becomes 



dtf^{\'P\,t) 



t = o 



d^k 



1 



1 



1 



~ 2 J (27r)3 2E^{p) 2Ex{k) 2£;^(p - k) 
X 27r6{E^{p) - Exik) - i^x(p - k)) 



fB{E^{p))[l + /B(i5^x(k)) + /^(i?x(p-k) 
fB{Ex{k))fE{E{p-k)) 



(9.46) 



where we remind the reader that energy conservation is artificially imposed. By virtue 
of this energy conservation, the first product of statistical factors in (9.46) satisfies the 
identity 

fB{E^{p))(l + fB{Ex{k)) + fE{Ex{p-k))) = fB{Ex{k))fE{E{p-k)) . (9.47) 



As a result, the RHS of (9.43) is exactly zero at t = 0. Consequently, the energy-non- 
conserving evanescent regime plays a fundamental role in the description of the evolution 
of the system J^. 

We now wish to show that the master time evolution equations (9.41) and (9.44) 
of our perturbative approach describe a non-trivial evolution of the system J/". For this 
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Figure 16: Numerical estimates of dtf^{\p\,t), as a function of Mt, assuming that the 

(c) 

X statistical distribution functions /^^ (|p|,t) on the RHS of the time evolution equa- 
tion (9.41) maintain their equilibrium Bose-Einstein form for all times. 



purpose, we consider the time evolution equation of the heavy scalar $ in (9.41) and assume 

(c) 

that the x statistical distribution functions fy^ (|p|,t) of the RHS remain in their initial 
equilibrium forms for all times, i.e. /;^*^^(|p|, t) = /b(-Ex(p))- With this assumption, we 
see from Figure 16 that the RHS of (9.41) is non-zero for early times. Thus, the system 
=y is indeed perturbed from its non-interacting equilibrium by the evanescent processes 
described by our perturbative time evolution equations (9.41) and (9.44). 
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Let us finally have a closer look at the early-time behaviour that immediately follows 
after the switching on of the interactions. We expect that this prompt behaviour is domi- 
nated by the ultra-violet contribution to the phase-space integral on the RHS of (9.41) due 
to the Heisenberg uncertainty principle. From Figure 13, we see that this shock regime 
contributes dominantly to the prompt evolution of the infra-red modes of /$(|p|,t). In this 
regime, we take |p| = and m = in (9.41) and introduce the ultra-violet cut-off A in the 
limits of the phase-space integrals of the RHS. Assuming that the x distribution functions 
are tempered, vanishing faster than a power law for large momenta, we may ignore their 
contribution in the ultra-violet limit. Hence, dtf^{\p\ = 0,t) can be approximated by 

2 

dtM\p\=0,t) ~ -^^^3^(1 + 2/$(|p| =0,t))^lim^Si2(At) . (9.48) 

Clearly, dtf^{\p\ = 0, t) vanishes in the limit t <C 1/A — )■ 0, as we expect. Expanding Si (At) 
about At = 0, i.e. for times infinitesimally close to zero, we obtain 

2 

dtM\p\=0,t) ~ + 2/$(|p|=0,0)) lim {Atf. (9.49) 

For small but finite times, ultra-violet contributions are rapidly varying. It is apparent 
from this discussion and Figures 13 and 16 that the effect of the transient behaviour of the 
system upon its subsequent dynamics is significant, since the quantum memories persist on 
long time scales Mt ^ 1. We may therefore conclude that as opposed to other methods 
relying on a truncated gradient expansion, our approach consistently describes this highly- 
oscillatory and rapidly-evolving early-time behaviour of the system. 

9.5 Inclusion of Thermal Masses 

In this section, we describe how local thermal-mass corrections may be incorporated con- 
sistently into our approach. 

The local part of the one-loop x self-energy shown in Figure 4 has the explicit form 

U'°<'\p,p',if;ii) = -^(27r)%(^)(p - p')(2vr;u)2^e*(^'«-fo)?/ 

"^'^ ^ ' (27r)4<5W(A: - k') 



{2-kY (27r)4 \ k'^ - m? + ie 
+ 27r<5(A;2 - m^)\2kQ\^/'^ f^{k,k\t)e''^'''' ' ^''^^~'f\2kX''^'2^5{k''^ - m^) ] , (9.50) 



with d = A — 2e (cf. Appendix F). The first term yields the standard zero-temperature 
UV divergence, which is usually removed by mass renormalization. The second term yields 

n'°<^\p,p\if-li) = -{27:)H\'\p - p')e^(P^>-Po)ifmUir,U) , (9-51) 
where the time-dependent thermal mass mtii{tf;ti), given by 



2 r. A /■ d^k 1 fd^k' 1 



X (/^(k,k',t)e^[^(i')-^W/ + /f (-k,-k',t)e-*[^('')-^W/) , 
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is UV finite. Notice that a unique thermal mass may not be defined in the spatially 
inhomogeneous case due to the explicit dependence on if . 

From the Schwinger-Dyson equation in (4.32), the inverse quasi-particle CTP propa- 
gator A~^^f^{p,p' ,if;ii) of the complex scalar field x takes the form 

^x!abiP^P'^^f'^^^ = (27r)45f)(p-pOe*^^''-^«^*^[(p''-"^t\(i/;i.))^ab + ^eU] , (9.53) 

in which we have assumed mth(i) ^ fn- If the quartic x self-interaction is switched 
on sufficiently long before t = 0, then for t > 0, we may replace the 5t function in the 
inverse quasi-particle propagator (9.53) by an exact energy-conserving delta function. This 
imposition of energy conservation constitutes a quasi-particle approximation, allowing us 
to invert (9.53) exactly, using the arguments of Sections 3.2 and 4.3, to obtain the following 
quasi-particle x propagators: 

A^^^\p,p\if-U) (9.54) 



(p2 - ml^{if,ii) + ie)-i -i27re(-po)5(p' - ml^iifM 



(27r)^<5(^)(p - p') 

-i27r\2po\'/H{p^ - m?h(ty;t,))/x(P,P',t)e*^^"'"°^*^2^|2p^|i/25(p'2 - m,\(t/;t,)) 



1 1 
1 1 



If the subsystem oS^^ is in a state of thermodynamic equilibrium at the initial time 
t = 0, the thermal mass reduces to the known result 

ml{t = 0) = — , (9.55) 

for m = and fj, <^T. In order to describe completely the dynamics of the combined system 
=5^ = U we couple the evolution of the thermal mass mth(t) to the perturbative 
time evolution equations (9.41) and (9.44) for the $ and x statistical distribution functions. 
This is achieved by differentiating (9.52) with respect to t, such that 

dtmthit) = ^ ^— -( dtfJ\k\.t) + dtf^(\k\,t)) . (9.56) 

^' 2mth(t)y (27r)3 2£;x(k) 2V " ' t-'x i' V ^ ^ 

Here, 9i/;^(|k|,t) and 9(/^(|k|,t) depend implicitly upon mth(t) via the coupled first-order 
differential equations (9.41) and (9.44), derived using the quasi-particle x propagators 
in (9.54). On dimensional grounds, we may estimate that the leading contribution to 
dtmfh{t) is of order X^/^gT. The latter estimate provides firm support that our time 
evolution equations may consistently incorporate thermal masses and so describe the main 
non-perturbative dynamics of the system. 
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10 Conclusions 



We have developed a new perturbative approach to non-equihbrium thermal quantum field 
theory. Our perturbative approach is based upon non-homogeneous free propagators and 
time-dependent vertices, which explicitly break space-time translational invariance and 
properly encode the dependence of the system on the time of observation. The forms 
of these propagators are constrained by fundamental field-theoretic requirements, such as 
CPT invariance of the action, Hermiticity properties of correlation functions, causality 
and unitarity. We have shown that our perturbative approach gives rise to time-dependent 
diagrammatic perturbation series, which are free of pinch singularities. The absence of 
these pinch singularities results from the systematic inclusion of finite-time effects and the 
proper consideration of the time of observation. We emphasize that this is achieved without 
invoking ad hoc prescriptions or effective resummations of finite widths. In our formalism, 
we have derived the new master time evolution equations (8.13) and (8.14) for particle 
number densities and statistical distribution functions, which are valid to all orders in 
perturbation theory and to all orders in gradient expansion. Furthermore, the master time 
evolution equations (8.13) and (8.14) respect CPT invariance and are also invariant under 
time translations of the CTP contour. As opposed to other methods, in our perturbative 
approach, we do not need to employ quasi-particle approximation and no assumption is 
necessary about the separation of time scales. 

We have shown how the effect of a finite time interval since the start of evolution of the 
system leads to violation of energy conservation, as dictated by the Heisenberg uncertainty 
principle. Our approach permits the systematic treatment of the pertinent generalized 
kinematics in this evanescent regime. We have found that the available phase space of 
1 — )• 2 decays increases and would-be kinematically disallowed processes can still take 
place at early times, contributing significantly to our time evolution equations. Within 
a simple scalar model, we have illustrated that these kinematically forbidden evanescent 
processes are the 2 — )• 1 processes of Landau damping, as well as 3 — )• processes of 
total annihilation into the thermal bath. The processes of Landau damping and total 
annihilation have been shown to contribute promptly to the particle width, to a level as 
high as 20%. The switching on of the interactions leads to a sudden deposition of energy 
in the system, which manifests itself as a rapid change in both the particle width and the 
collision terms of the time evolution equations. These early-time effects give rise to an 
oscillating pattern, which persists even at later times. We have demonstrated that these 
late-time memory effects exhibit a non-Markovian evolution characterized by oscillations 
with time-dependent frequencies. The latter constitutes a distinctive feature of proper 
non-equilibrium dynamics, which is consistently predicted by our perturbative approach. 

We note that the rapid transient behaviour of the system makes the method of gradient 
expansion unsuitable for early times. We emphasize that in our approach no assumption 
was made to the relative rate of thermalization of either the statistical or spectral behaviour 
of the system. For the considered initial conditions of thermal equilibrium, we have found 
that the spectral evolution resulting from evanescent contributions is critical to the early- 
time statistical dynamics of the system. Consequently, it would have been inappropriate 
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to assume a separation of time scales at early times. A more accurate numerical solution 
to our time evolution equations turns out to be computationally intensive and may be 
presented elsewhere. 

Finite-time effects, which are systematically incorporated in our perturbative approach, 
are also relevant to many-body systems whenever a natural characteristic time scale for 
perturbations arises in such systems. Thus, in addition to possible applications to reheat- 
ing and preheating, of particular interest are first-order phase transitions. For instance, 
such characteristic time scales result from the bubble wall velocity of nucleation in the 
electroweak phase transition, e.g. see [78, 79]. In the vicinity of the bubble wall, evanescent 
contributions may weaken the constraints of decay and inverse decay thresholds and so af- 
fect the washout phenomena and the generation of relic densities. This evanescent regime 
is particularly relevant to prethermalization [80] and isotropization [81] time scales, which 
are known to be shorter than the time scale of thermalization. 

It is straightforward to generalize our perturbative approach to theories that include 
fermions and gauge fields. Thus, it would be very interesting to extend the classical ap- 
proaches [82, 83] to kinetic equations of particle mixing to non- homogeneous backgrounds, 
by including finite-time evanescent effects in line with the non-equilibrium formalism pre- 
sented in this paper. 
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A Propagator Properties and Identities 



In this appendix, we give a detailed summary of the transformation properties and identities 
that relate the various propagators defined in Section 2. In detail, the pertinent two-point 
correlation functions for the complex scalar field x ^'^^ given by 
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The definitions of the charge-conjugate propagators follow from the unitary transformation 



Ulx{x)Uc = x'-ix) = vx\x), (A.2a) 

ul.x\x)Uc = X^\x) = rTxix) , (A.2b) 

where the complex phase rj satisfies |r/p = 1. 

It follows from the definitions that the propagators satisfy the transformations listed 
below under charge- and Hermitian-conjugation: 

A{x,y) = -A*{y,x) = -A^{y,x), (A.3a) 

Ai(x,y) = -Al{y,x) = Af(y,x), (A.3b) 

Av{x,y) = A*T,{y,x) = A^{y,x) , (A.3c) 

AR(a;,y) = Ag*{x,y) = A\{y,x) = Ag(y,x), (A.3d) 

A>(x,y) = -A*>(y,x) = A^(y,x) = -A^*(x,y), (A.3e) 

AF(x,y) = A^{y,x) = -AJ,(y,x) = -Ag*(x,y), (A.3f) 

where the action of charge-conjugation is trivial in the case of a real scalar field. In the 
double momentum representation, these identities take the form: 

A{p,p') = -A*ip',p) = -A^i-p',-p), (A.4a) 

Ai{p,p') = -Al{p',p) = Af{-p',-p), (A.4b) 

Ar{p,p') = A1p{p\p) = A^{-p',-p), (A.4c) 

Ar(p,p') = Ag*{-p,-p') = A\{p',p) = A<l{-p',-p), (A.4d) 

A^{p,p') = -Al{p',p) = A<^{-p',-p) = -Af{-p,-p'), (A.4e) 

Af{p,p') = Ag{-p',-p) = -AUp,p) = -A^*{-p,-p') . (A.4f) 
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Finally, in the Wigner representation, the propagators satisfy the following properties: 



A(g,X) = -A*{q,X) = -/\^{-q,X), (A.5a) 

Ai(g,X) = -A\{q,X) = A?{-q,X), (A.5b) 

AK{q,X) = Ag*{-q,X) = Al{q,X) = A^(-g,X), (A.5c) 

A>(g,X) = -Al{q,X) = A^(-g,X) = -Af{-q,X), (A.5d) 

AF{q,X) = Ag{-q,X) = -AUq,X) = -Ag*(-g,X). (A.5e) 

We also list the following set of useful identities: 

A{x,y) = A>(x,y) - A<(x,y) = Ar(x,?/) - AA{x,y) , (A.6a) 

= e(a;o - yo)(^AF{x,y) - AY){x,y)^ , (A. 6b) 

Ai{x,y) = A>(x,y) + A<(x,y) = AF{x,y) + Aj){x,y) , (A.6c) 

AR(A){x,y) = AF{x,y) - A<(>)(x,y) = -AD{x,y) + A>(<)(x,y), (A.6d) 

^F(D){x,y) = ^((-)2Ap(x,y) + A>{x,y) + A<(x,y)) , (A.6e) 

Ap(x,y) = ^(^An{x,y) + AA{x,y)^ = ^(^AF{x,y) - Ad(x,2/)). (A.6f) 



We note that analogous relations hold for the corresponding self-energies and that these 
identities and relations are true for free and resummed propagators. 

B Correspondence between Imaginary and Real Time Formalisms 

In this appendix, we briefly outline a number of relevant details of the Imaginary Time 
Formalism (ITF) of thermal field theory. These details are discussed in the context of the 
real scalar field theory introduced in Section 2. Subsequently, we identify the correspon- 
dence of the ITF $-scalar propagator and its one-loop non-local self-energy with results 
calculated explicitly in real time, using the CTP formalism of Section 3 in the equilibrium 
limit discussed in Section 5. 

The equilibrium density operator p^q in (5.6) permits a path-integral representation in 
negative imaginary time. The ITF generating functional is 

Z[J] = y"[d$] exp(^ - Sm + J^^dT.Jd^x Jix)^x)^ , (B.l) 

with action 

= j^T^ j d^^ (^\d^^{x)d^,^{x) + \M^^''{x) + y^\x) + ^\^\x)) . (B.2) 

We emphasize the restricted domain of integration over the imaginary time Tx G [0, /?]. 
Four-dimensional Euclidean space-time coordinates are denoted by a horizontal bar, i.e. 
Xfj_ = (tj;, x). In the limit f3 — )■ oo, (B.l) is precisely the Wick rotation to Euclidean space- 
time of the Minkowski-space generating functional via the analytic continuation xq — )• — ir^. 
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The free imaginary-time propagator A*^ may be written as follows: 



(B.3) 



i 

where 



^(^^^'P) = .2,-9, ..9 , fV/.. ■ _^ • (B-7) 



+ P + ^ 

is the so-called Matsubara propagator. The discrete Matsubara frequencies coi = 2-ir£/[3, 
i G Z, arise from the periodicity of the imaginary time direction in order to satisfy the 
KMS relation [cf. (5.15)] 

A°(x-y) = A°(x-y + /3). (B.5) 

The resummed Matsubara propagator is given by the imaginary-time Schwinger-Dyson 
equation 

A-\iuje,p) = A^^-\iuje,p) + U{ioJe,p) , (B.6) 

where ti{u}£, p) is the imaginary-time self-energy. Equation (B.6) may be inverted directly 
to obtain 

1 

w| + p2 + M2 + n{iuje,p) 

The free Matsubara propagator in (B.4) may also be written in the following spectral 
representation: 

AO.. ^ . fdkp A°(fco,p) 

A iiuji,p) = -I — — , (B.8 

J 2-K loje - ko 

where 

iA^{ko,p) = 2TTe{kQ)5{kl - p^ - M^) (B.9) 

is the single-momentum representation of the free Pauli-Jordan propagator, consistent with 
(2.16a). Making the analytic continuation ioji pQ + ie to real frequencies and comparing 
with the spectral representation of the retarded propagator Ar in (2.22), we may establish 
the correspondence 

A^{iuje^Po + ie,p) = -A^ip) . (B.IO) 

This correspondence must also hold for the resummed Matsubara propagator A via (B.6). 
As such, the analytic continuation of the ITF self-energy 11, 

Iliiuje^Po + i€,p) = -U^iip) , (B.ll) 

yields the equilibrium retarded self-energy IIr. 

Thus, in thermodynamic equilibrium, an exact correspondence can be established be- 
tween the ITF and the real-time approach of the CTP formalism, by means of retarded 
propagators and self-energies. The full complement of propagators exhibited in Table 2 and 
the corresponding self-energies may be obtained using the constraints of causality (2.24) 
and unitarity (2.28) in combination with the KMS relation (5.16) and the condition of 
detailed balance (5.25). The relationships between the retarded and CTP propagators and 
self-energies are listed explicitly in (5.19) and (5.26). 
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To illustrate the correspondence in (B.ll), we consider the one- loop bubble diagram of 
the real scalar theory in (2.1). The real and imaginary parts of the retarded self-energy IIr 
may be calculated in the CTP formalism from the time-ordered 11 and positive-frequency 
Wightman n> self-energies, respectively, using the relations in (5.26). 



The real part of the one-loop retarded self-energy is given by 



Ren{J^(p) = Re 



X-igY [ d^k 



(B.12) 



(2vr)^ 

where iAp(/c) is the equilibrium CTP Feynman propagator in (5.10a). Explicitly, we have 



/2E2 ( 
(p_fc)2-M2 "^(^Q " "1^1) (5 + /b(A;c 

dfco . 2 _ (Po - ^0 - a2E-i)y \ + /b(po-^o) 

. 



in^i^(p) = iA°(A;)iAO(p-A;) . (B.15) 



"2 = ±1 fcVA^2 

where E\ = -E(k) and E2 = £'(p — k). The integral subscripts remind us that the 
integration around the on-shell poles are understood in the Cauchy principle value sense. 
Integration over fco yields the result 

Ren^ (P) = -y Ey po - a,E, - a2E2 ' ^^-^'^ 

{a} 

where we have used the short-hand notation {a} to denote summation over ai, 02 = ±1. 
The one-loop positive-frequency Wightman self-energy is given by 

(2vr) 

Using the signum form of the equilibrium positive-frequency Wightman propagator iA^{k) 
in (5.10b), we obtain 

^>\p) = iT^g^^ jfl^ j '^^^ 4Ei^2 ^^^^ ~ "i^i^'^^^o - ko - a2E2) 

{a} 

X ( 1 + /b(A:o) ) ( 1 + /b(po - A;o) ) . (B.16) 
The product of Bose-Einstein distributions in (B.16) satisfies the following relation: 

/b(A:o)/b(po - ko) = fBiPo) ( 1 + /b(A:o) + /b(po - A^o) ) • (B.17) 
Thus, upon integration over /cq, we find 

{q} 

X (1 + hiaiEi) + /B(a2^2)) . (B.18) 
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Using the relation in (5.26c), the imaginary part of the one-loop retarded self-energy may 
then be written down as 

{a} 

X (l + /B(ai^i) + /B(a2^2)) . (B.19) 



In the ITF, the one- loop self-energy is given by 

\2 ±^ r 



n«(fc.,,p) = ^ A(iu;„,k)A(i(c.,-^„),p-k). (B.20) 



n = — oo 



After performing the summation over n (see for instance [68]), we obtain 

K i^v) 2 Z-^ J (2TTf AE1E2 iui - aiEi - a2E2 ^ ' 

Making the analytic continuation iuji pQ + ie and subsequently extracting the real and 
imaginary parts, the result in (B.21) agrees with (B.14) and (B.19) via the correspondence 
in (B.ll). 



C The Complex Scalar Field 

This appendix describes the generalization of the results in Sections 4.3 and 5 to the case of 
a complex scalar field x- particular, we derive the full complement of non-homogeneous 
free propagators and expand upon the ITF correspondence identified in Appendix B. 
Our starting point is the complex-scalar Lagrangian 

C{x) = d^x\x)d^x{x) - mWx)x{x) - \\[x\x)x{x)f . (C.l) 

In analogy to (2.11), the complex scalar field x{x) may be written in the interaction picture 
as 

X{x-ru) = y"dnp (a(p,0;t,)e-*^(P)^«e*P-^ + 6t(p, O; t,)e^^(P)^e-^P-'' ) , (C.2) 

where a"l"(p,0;ti) and a(p, 0;tj) (6"^(p, 0;tj) and 6(p, 0;ti)) are the interaction-picture par- 
ticle (anti-particle) creation and annihilation operators, respectively. Under C-conjugation 
[cf. (A. 2a)] these creation and annihilation operators satisfy the transformations 

ul;a{ip,t-ii)Uc = 7]b{p,i;ii) , ulb\p,i;ii)Uc = r]a\p,i;U) . (C.3) 

Introducing the four-dimensional LIPS measure from (2.3), (C.2) and its Hermitian conju- 
gate may be recast in the form 

xix;U) = J^e-^P-^xip;U), xHx-k) = j e'^^'^ x^-p-.U) , (C.4) 
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where the Fourier amphtudes are given by 

X{p;ii) = 2tt6{p^ - m2)(0(po)a(p,O;t~) + e{-pop {-p,0-U)^ . (C.5) 

For the real scalar field, the quantization scheme was only depending on the restriction 
placed upon the form of the field commutator. In the case of the complex scalar field, we 
have two degrees of freedom to fix. Thus, we begin with the following two commutators of 
interaction-picture fields: 

[x(^), X{y)] = 0, [x{x), xHy)] = il^\x,y-m^) , (C.6) 

where the Pauli- Jordan propagator has precisely the form in (2.15). In analogy to the real 
scalar field, we may derive from (C.6) the equal-time commutation relations 

iA0(x,y;m2)|^„^^,^- = [x(t,x), x\ly)] = 0, (C.7a) 

a,„iAO(x,y;m2)|^o = ,o = i = [^^(i,x), xHly)] = - y) , (C.7b) 

9^„iAO(x,y;m2)|^o^^o^- = [x(t,x), vr(t,y)] = i(5(3)(x - y) , (C.7c) 

d^,dy,i^\x,y-m'')\^,^y,^-^ = [7rt(t,x), ^(t,y)] = 0, (C.7d) 

where 7r(x) = dx^x^x) is the conjugate momentum operator. The particle and anti-particle 
creation and annihilation operators necessarily satisfy the algebra 

[a(p,t), at(p',t')] = [6(p,t), 6t(p',?)] = (2^)32i^(p)j(3)(p _ p/)g-.E(p)(?- ?) , (c.8) 

with all other commutators vanishing. We stress again that the phase factor e"*^^^^^* ~ on 
the RHS of (C.8) has appeared due to the difference in microscopic times of the interaction- 
picture creation and annihilation operators. 

In addition to the C-conserving propagators listed in Appendix A, we may also define 
C-violating propagators. As an example, the C-violating Hadamard propagator A^0(x,y) 
would read: 

iAi^(x,y) = {{x{x), x{y)]) , (C.9) 

which satisfies 

Ai^(x,y) = Ai^(2/,x) = VAS(a:,y). (C.IO) 

This Hadamard correlation function may, in general, be non-zero for early times, thus per- 
mitting extra C-violating evanescent processes in addition to those described in Section 9. 

In analogy to (4.53), we write the following set of EEVs of two-point products of 
particle and anti-particle creation and annihilation operators: 

{a\p\tf;k)a{p,if-ti))^ = 2^(p, p')/(p, p', t) , (C.lla) 

{h{p',if-ii)a{p,tf-ii))^ = 2^(p,p')<7(p,p',t) , (C.llb) 

(a(p',t/;ti)a(p,t/;ti))^ = 2(^(p, p')/i(p, p', t) , (C.llc) 

{h\p',if-ii)a{v,if,ii))t = 2<f(p,p')d(p,P,t) , (C.lld) 
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where the remaining EEVs are obtained by Hermitian and charge conjugation. The four 
statistical distribution functions /, g, h, d satisfy the following identities: 

/(p,p',t) = r(p',p,t), (C.12a) 

g(p,p',t) = 5^(p',p,t) , (C.12b) 

/i(p,p',t) = /i(p',p,t) , (C.12c) 

d(p,p',t) = r/2d^*(p',p,t). (C.12d) 

The non-homogeneous free propagators of the complex scalar field x may be written 
as listed in Table 2 with the substitution of the following ensemble function: 

~f{p,p\t) = e{po)e{p'^)fip,p',t) + 9(-po)9i-p'o)f*{-p,-p',t) 

+ 0{po)9{-p'o)g{p,-p',t) + e{-po)e{p'o)g''*{-p,p',t) , (C.13) 

satisfying the relations: f{p,p',t) = f^{—p',—p,t) = —p', t) in accordance with 

(A. 4). The free C-violating Hadamard propagator A^0{p,p' ,if;ii), defined in (C.9), may 
be written down as 

A,^{p,p',if,k) = -i27rS{p^ - m^)\2po\'/'2d{p,p',t)e'^P'^-P'o^'f\2pX/'27r5{p'' - m^) , 

(C.14) 

where we have defined the C-violating ensemble function d{p,p',t). Evaluating the EEVs 
directly, we find 

d{p,p',t) = e{po)9{p'o)d{p,p',t) + e{-po)e{-p'o)v^d''*{-p,-p',t) 

+ e{po)e{-p'o)h{p,-p',t) + e{-po)e{p'o)ri^h''*{-p,p',t) . (c.is) 

satisfying the relations: d{p,p',t) = d^{—p\—p,t) = rf'd^*{—p,—p',t). 

The inclusion of the C-violating distribution functions requires the addition of the C- 
violating source lab to the expansion of the kernel of the density operator in the CTP gen- 
erating functional, according to our discussion in Section 4.1. Omitting the t-dependence 
of the sources for notational convenience, the CTP generating functional for the complex 
scalar field x takes the form 

ma,kab,lab\ = -ih\n j j[d{x''\x'')] 

X expjl 

+ ^^xd^x' (^x''\x)Kb(x,x')x\x') + x\x)k^,,{x,x')x'\x') 

+ x%x)li,{x,x')x\x') + x"^(x)/a6(x,x')x'^(x')) + ••• } • (C.16) 

The bi-local sources kab and lab necessarily satisfy the identities 

kab{x,x) = kl^{x',x), lab{x,x) = lba{x' , x) = rj^l^^ {x, x') , (C.17) 



Six^^x''] + / d^x (il(x)x"(x) + xi{x)f{x) 
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to ensure that the exponent of the generating functional is CPT- invariant. The subsequent 
derivation of the effective action then follows analogously to Section 4. 

The Lagrangian in (C.l) is invariant under the global U{1) transformation 



entailing the conserved Noether current 
with corresponding conserved charge 



(C.18) 
(C.19) 



: Qixo) : 



d^x : jo{x) : 



(p,0)a(p,0) - 6t(p,o)&(p,0)) , (C.20) 



where : : denotes normal ordering. The existence of this conserved charge necessitates the 
introduction of a chemical potential and, as such, the equilibrium density operator is of 
the grand-canonical form 



-l3{H-f,Q) 



(C.21) 



In the presence of this chemical potential, the KMS relation in (5.15) generalizes to 

A>(:cO-yO,x-y) = e-^''A<{x^-y^ + i(3,^-y) (C.22) 

or, in the momentum representation, 

Ay{p) = e^(Po-^)A<(p) . (C.23) 

Proceeding as in Section 3, we find that the final constraint on f{p), generalising (5.17), is 

f{p) = e(Po)/B(po) + Oi-po)f^{-po) , (C.24) 

where f^\po) = (e'^'^""^"''^^] — 1)~^ is the particle (anti-particle) Bose-Einstein distribution 
function. In equilibrium, translational invariance is restored and the elements of the free 
CTP propagator may be written in the single momentum representations 

-1 



iAUp) 



H P 



m + le 



+ 2^ ( e(po)/B(po) + e{-po)fE{-po) Wp 



m 



iAlip) 
iAO(p) 



27r 
2tt 



m 



m 



(C.25a) 
(C.25b) 
(C.25c) 



0{po)[l + /b(po)) + 0{-po)fE{-Po)\S{p^ 

'o{po)fBiPo) + 0{-po)[l + f^{-Po))]s{p^ 

In order to define the ITF generating functional for the grand canonical partition 
function in (C.21), one may consider the Hamiltonian form of the path integral directly 
(see for instance [68]). We may write 



[dinlTT)] [d(x^x)] 



exp 









f d^x 







?^(^(t),^{t)) 



- i (^Tr{x)dr^xix) + tt'' {x)dr,xHx)^ - if^ (^T^Hx)xHx) - TT{x)xix)^ 

- jHx)x{x)-j{x)x{x)]] . (C.26) 
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Expanding the fields and conjugate momenta in terms of two real degrees of freedom Xi,2{x) 
and vri^2(2;) as 



Xix) 



V2 



(C.27) 



we may analytically calculate the Gaussian integrals over 711,2(2:), yielding 

-J^^dr, J d^X ( {dr^ + fl)xHx){dr^ - f,)x{x) 

+ Vx"^(x) • Vx(x) + m^xHx)xix) - j\x)xix) - j{x)xHx)^ 



m = [dix^,x)] 



(C.28) 



In order to derive the form of the ITF x-propagator, we insert into (C.28) the Fourier 
transform 

+00 



X{x) 



1 ^ r i\uj,T^+p.y 



'x(«w^,p) 



(C.29) 



where oji are the discrete Matsubara frequencies described in Appendix B. The effect of 
the chemical potential is to shift the poles of the Matsubara propagator, which becomes 



A°(ia;£ - ^,p) 



(C.30) 



In generalization of the correspondence in (B.IO), the equilibrium retarded propagator Ar 
is obtained by the analytic continuation iu^ — fj, ^ pQ + ie via 



A'^iiui - Po + ie,p) = -A^{p) 



(C.31) 



Correspondingly, the equilibrium retarded self-energy IIr is given by 



U{iu}i - fj,^ po + ie) 



(C.32) 



thereby generalising (B.ll). 



D Non-Homogeneous Density Operator 

In (D.l), we highlight the full form of the series expansion of the general Gaussian-like 
density operator described in Section 4.3. Symmetric and asymmetric multi-particle states 
are built up by summing over all possible convolutions of the W amplitudes. To facilitate 
our presentation task, the time dependence and phase-space integrals have been omitted. 
In fact, all momenta are integrated with the LIPS measure given in (2.3). 
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p = (l + ^W^io(pi : 0)T^oi(0 : Pi) + ^H^2o(pi, P2 : 0)W^o2(0 : P2, Pi) + •••){ |0)(0| 
+ ( - W^io(ki : 0) + ^W^ii(ki : qi)W^io(qi : 0) 

+ ^W^2o(ki,qi : 0)W^oi(0 : Qi) + • • • ) |ki)(0| 
+ ( - Woi(0 : k'l) + ^Woi(0 : qi)Wii(qi : k[) 

+ ^Wo2{0 : k;,qi)W^io(qi : 0) + • • • ) |0)(k'i| 
+ (^{27rf2E{ki)d(^\ki - k'l) - W^ii(ki : k[) + Wwi^i : 0)VFoi(0 : 



+ -Wii(ki :qi)VFii(qi :k;] 



+ ^W^2o(ki,qi : 0)W^o2(0 : qi,k'i) + • • • ) |ki)(k'i| 
+ l{- W^2o(ki,k2 : 0) + VFio(ki : 0)W^io(k2 : 0) + ^W^ii(ki : qi)W^2o(qi, k2 : 0) 

+ ^M/ii(k2 :qi)W^2o(qi,ki :0) + • • • ) |ki, k2)(0| 
+ l{- ^02(0 : ki,k'2) + Woi{0 : k[)Woi{0 : k'2) + ^W^o2(0 : k'^, qi)W^ii(qi : k^) 

+ ^Wo2{0 : k^,qi)t^ii(qi : k'^) + • • • ) |0)(k;,k'2| 
+ ^(Wii(ki : k;)Wio(k2 : 0) + W^ii(k2 : k;)W^io(ki : 0) 

+ W2o(ki,k2 : 0)VFoi(0 : k'l) + • • • ) |ki,k2)(k;| 
+ ^(Mi(ki : k'i)W^oi(0 : k^) + W^ii(ki : k'2)W^oi(0 : k'J 

+ Wo2(0 : k;,k'2)W^io(ki : 0) + • • • ) |ki)(k;,k'2| 
+ ^(T^2o(ki,k2:0)W^io(k3:0) + • • • ) |ki, ks, k3)(0| 
^(Wo2(0:k;,k'2)T^oi(0:4) + • • • ) |0)(k;, k'2, 4| 

(27r)32^(ki)5(3)(ki - k'l) - T^ii(ki : k'l) + H^io(ki : 0)Woi(0 : k'^) 

+ ^H^iilki : qi)H^ii(qi : k'^) + ^W^2o(ki,qi : 0)W^o2(0 : qi,k'i) + • • • ) 
: ((27r)=^2£;(k2)(5(3)(k2 - k^) - W^ii(k2 : k^) + W^io(k2 : 0)Woi(0 : k'2) 

+ ^l^ii(k2 : q2)l^ii(q2 : k^) + ^M/2o(k2,q2 : 0)W^o2(0 : q2,k'2) + • • • ) 
^( - VF2o(ki,k2 : 0) + VFio(ki : 0)T^io(k2 : 0) + ^T^ii(ki : qi)W2o(qi,k2 : 0) 



+ 



+ 



+ 4 



+ 2^11(^2 : qi)W2o(qi,ki : 0) + 



( - Wo2{0 : k'i,k^) + ^01(0 : K)Woi{0 : k^) + \wo2{0 : k^, qi)T^ii(qi : k'2) 
+ ^W^o2(0:k'2,qi)W^ii(qi :kl) + •••)] |ki, k2)(k;, k'21 + •••}. (D.l) 
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E Kadanoff— Baym Equations 

In this appendix, we outline the derivation of the KadanofF-Baym equations [52, 53]. The 
resulting time evolution equations are included here for comparison with the new results 
described in Sections 8 and 9. The Kadanoff-Baym equations are obtained by the partial 
inversion and subsequent gradient expansion of the Schwinger-Dyson equation derived in 
Section 4.1; see also our discussion in Section 4.2. We omit the if- and tj-dependence of 
all propagators and self-energies for notational convenience. 

Inverse Fourier transforming (8.5), we obtain the coordinate-space representation of 
the partially inverted Schwinger-Dyson equation 

- (n^. + M2)A^(x,y) +^d^^ (n^(x,z)Ap(z,y) + U^ix, z)A^{z,y)) 

= I [ d^z (n<(x,z)A>(z,y) - n>(x,z)A<(z,y)) , (E.l) 

where we have also included the positive-frequency contribution for completeness. Intro- 
ducing the relative and central coordinates and via (4.38), we may show that the 
d'Alembertian operator may be rewritten as 

□2 = + + hA\ . (E.2) 

Performing a gradient expansion of the Wigner transform of (E.l), as described in Sec- 
tion 4.2, we obtain 



^ exp[-i (q.X + + 20+ ^ 

({n^(g + Q/2,X)}{Ap(g-Q/2,X)} + {np(g + Q/2, X)}{A^(g - Q/2, X)} 

= ^ (Q) exp [ - z ( Q . X + + 20+ ^ 

({n<(<7 + Q/2,X)}{A>(g-Q/2,X)} - {n>(g + Q/2, X)}{A<(g - Q/2, X)} 

(E.3) 

The diamond operators 0+{»}{»} are defined in (4.46) and (4.47) of Section 4.2. Subse- 
quently separating the Hermitian and anti-Hermitian parts of (E.3), we find the so-called 
constraint and kinetic equations 

e -M^ - ^□2,)A^(g,X) + 1^ (2vr)%(^)(Q) cos(q.X + + 20^_^ ) 

({n^(g+f,X)}{Ap(g-f,X)} + {np(g+f,X)}{A^(g-f,X)}) 

= \\wf ^^""^'^^'^^^^ sin[Q.X + 0,-^ + 20+^) 

({n>(g+f,X)}{A<(g-2,X)} + {U^{q + X)}{A^{q - X)}) , (E.4a) 
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({n^(g+f,X)}{Ap(g-f,X)} + {np(g+f,X)}{A^(g-f,X)}) 

({n>(g+f,X)}{A<(g-2 X)} - {n<(g + f , X)}{A>(<? - f , X)} ) , (E.4b) 
respectively. 

In the late-time limit t — )• cxd, the Q dependence is removed and the microscopic 
violation of energy conservation resulting from the uncertainty principle is neglected. As 
was established in Section 9, this microscopic violation of energy conservation is significant 
to the early-time dynamics of the system. Subsequently, keeping terms to zeroth order 
in the gradient expansion, the above expressions reduce to the more familiar form of the 
Kadanoff-Baym equations [1, 42, 43, 52, 53]: 

(g2 _ M^)^^{q,X) = - (n^(g,X)Ap(g,X) + np(g,X)A^(<7,X)) , (E.5a) 
q-dx^^{q,X) = ^(n>(g,X)A<(g,X) - n<(g, X)A>(q, X)) . (E.5b) 

The kinetic equation (E.5b) is to be compared with the semi-classical Boltzmann transport 
equation (6.18) and the energy-conserving limit of our time evolution equations in (8.17). 



F Non-Homogeneous Loop Integrals 

In this final appendix, we outline the techniques we have been using to perform the loop in- 
tegrals that occur with the modified time-dependent Feynman rules of this new perturbative 
approach. In particular, we develop a method for dealing with the energy-non-conserving 
vertices that utilizes the Laplace transform. Lastly, we summarize the time- independent 
equilibrium and time- dependent spatially homogeneous limits of these integrals. 
We begin by defining a generalization of the zero-temperature Bq function [84] 



B^ = ^{q,m,,m2) = (2^^"'^/^ 



1 



1 



k'^ — ml + ie {k — qY — -|- ie 



(F.l) 



To this end, we define the non-homogeneous CTP Bq function, which may be written in 
the 2x2 matrix form as 



5, 



ab 



B(] Bn 



-B* 



(F.2) 
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Figure 17: The non- homogeneous Bq function. 

consistent with the propagators and self-energies of the CTP formalism in Sections 3 and 4. 
The elements of Bq'^ are given by the following integral: 

X r?''^'^AO,(A:i, k[,mi,if, U)A^f{k2, k'^, ms, t/; ^i)r/^^^ 

(F.3) 

where .5^^ is defined in (3.49) and (3.50), A^'''^{p,p' ,t) is the free CTP propagator derived 
in Section 4.3 and the upper-case Roman C denotes charge-conjugate free CTP propagator. 
We reiterate that lower-case Roman CTP indices are raised and lowered by contraction with 
the SO(l, 1) 'metric' rjab = diag(l, —1) and that rjabc = +1 is for a = 6 = c = 1, 
Vabc = — lfora = 6 = c = 2 and r]abc = otherwise. 

In the definition (F.3), we have assumed that the statistical distribution functions are 
cut-off in the ultra-violet. Thus, any ultra-violet divergences anticipated from the super- 
ficial degree of divergence of the integral will be those that result from the homogeneous 
zero-temperature contribution. As a result, the dimensional regularization of the integral 
has been restricted to the ki dependence only. 

In order to deal with the product of time-dependent energy-non-conserving vertices, 
we make the following replacement: 

''^'^''^'^ = ^ ' ^ (X - yy + 4.^ (X + yr + 4.^ ' ^^-'^ 

where the RHS is the inverse Laplace transform and s G C is a complex variable. The 
Bromwich contour is chosen so that o" G M is larger than the real part of the right-most 
pole in the integrand, in order to ensure convergence. We then introduce the representation 
Bq^ {qi, q2,mi,m2, s) of the non-homogeneous CTP Bq function through 

r-a+ioo 

B''\qi,q2,mi,m2,if;ii) = / — e'^ B^''{qi,q2,mi,m2,s) . (F.5) 
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Note that BQ^{qi, q2,'mi, m2, s) is not the exact Laplace transform of -Bo^(q'i, (?2, "ii, m2, i), 
since we have not transformed the t-dependence of the distribution functions. We suppress 
the dependence of i?o''(9i; 92, "ii, "Z2, s) on if and ii for notational convenience. 
After making the replacement (F.4) in (F.3), we obtain the integral 

abr„ „ _ _ „N _ Q^o„,.^4-d i(<7?-oO)t, f fd'^ki d^k[ d^fcs 



i7r2 (27r)4 (27r)4 {2tt)^ 



X r/'^^'^AO,(A;i,A:;,mi,ty;tOA°)^(fe2,4,m2,t/;ti)^^^\ (F.6) 

in which the analytic structure of the product of formerly t-dependent vertex functions is 
now manifest. 

For conciseness, throughout this appendix we use the short-hand notation: 

{a} 

for summation over only the = ±1 and a'^ = ±1 that appear explicitly in an expression. 
On-shell energies are denoted by 



Ei = E,{k,) = ^kf + , E[ = Em = ^Jk'/ + m\ , (F.8) 
and on-shell four-momenta by 

k = kt = {a^E^M), % = k',^^ = {a',Em,k!,). (F.9) 
For the energy factor in (4.54), we use the notation 



Si = Si{k^M) = \l Ei{ki)Er{K) . (F.IO) 

Products of unit-step functions are denoted by 

,z) = e{x)e{y)---e{z) . (F.ll) 

Time-Ordered Functions 

The (a, 6) = (1,1) element of (F.3) coincides with the time-ordered function Bq, which we 
may separate into four contributions (suppressing all arguments): 

Bq = /« + + + . (F.12) 

These four contributions are the zero temperature limit (/*•*-*), the purely thermal term 
(/(»»)) and the cross terms (l(***'^) and /(***'')). 
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(i) The zero-temperature part may be extracted from the product of terms that remain 
in the limit of vanishing statistical distribution functions. This part is 



q2,mi,m2,s) 



d'^h d% d^fca d^A;^ 



in^ (27r)4 (27r)4 {2tt)^ 
X (2^)35(3)(qi - ki - k2)(27r)3<5(3)(q2 - k'^ - k'^) 



X 4s 



— m| + le 



{2TTf5^^\kl - k[ 



m\ + ie 



(27r)^(^(^)(A:2 - k'2) . (F.13) 



The k[^ and k'2^ integrations are performed by means of the delta functions, giving 



d'^-ifci d^k; d^ks d3k'2 . 



iir^ (27r)3 (27r)3 (27r)3 
X (2^)35(3)(qi - ki - k2)(27r)35(3)(q^ _ j^/ _ j^/) 



dfc^ dA:^ 



4s 



- 92^) + 4s^ 
1 



-1 



+ s^ 



(fc?)2 - + ie (A;0)2 



+ ie 



2 '"1 ~ '^2 

(27r)3<5(3)(ki - k'i)(27r)3<5(3)(k2 - k'2; 



(F.14) 



By virtue of the residue theorem, we may perform the k^ and /c2 integrations by closing 
contours in the lower-halves of the ki and ^2 complex planes. After collecting together the 
resulting terms, we find 



q2, mi, 1712,3 



,s 


) = 






EIE2 


1 




vr 





-2(27r)V'-V('??-''^')*/ J] 

{a} 



d'^-iki d^k'i d3k2 d^k'a 



(2^)"'-i (27r)3 (27r)3 (27r)3 



+ 4s' 



-1 



Ql 



ai -El + E2 



IS 



X (27r)35(3)(k, _ k;)(27r)35(3)(k2 _ k'2) , 
or, in terms of the Lorentz-invariant phase-space measures in (2.3), 

I'^'>{qi,q2,mi,m2,s) = -8(27r) V-''e*^^° " E / / dn^^UU^, dU^, dU^, 
X (2^)3<5(3)(qi - ki - k2) (2^)3^(3) (q2 - k'^ - k'2) 



(F.15) 



vr 



«0 - nO 



+ 4s" 



ai ii^i + i?2 



IS 



(27r)32E;<5(^)(ki - k'i)(27r)32E^(^(3)(]^2 _ j^/ 



2; > 



(F.16) 



where the superscript d—\ indicates that the ki integration is to be taken over a (d — 1)- 
dimensional phase-space. 
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The Laplace transform J-{s) of a complex-valued function F{t) has the following prop^ 
erties: 

1 



Ct[ReF]{s) 
Ct[ImF]{s) 



2 L 
1 

2il 



F{s) + [^F{s 
T{s) - (f{s*) 



(F.17a) 
(F.17b) 



We may then identify the dispersive and absorptive parts of the t-dependent Bq function 
with the parts symmetric and anti-symmetric in s, respectively. With this separation, 
(F.16) becomes 



{a} 



■ I dn^-i dHk; dHk, dHk^ 



1 r 



TT 



{ 



qi - 92 ) + 4s 



-1 



ai -El -I- E2 



2 2l~^ 



ai i?i -I- E2 



+ is 



X {27rf2E[5^^\ki - k[){2TTf2E'26^^\k2 - k'2) , (F.18) 

where the delineation between dispersive and absorptive parts is identified within the curly 
brackets of the fourth line. 

Isolating the dispersive part of this result and performing all but one of the remaining 
phase-space integrals, we find 



I)ispBo{qi,q2,mi,m2,s) = (27r) 



-5(3) (c 



)ei(g?-q§)V2 



X ^-(27r/i)^-M d'^-iki 



^1 + ^2 1 r/g?+g§ . ^2 
^ 2 — las 



El + E2 



E1E2 vrL^ 

(F.19) 

where Disp stands for the dispersive part. We compare this result to the form of the 
zero-temperature Bq function (F.l) after the ko integration has been performed: 



Bq °{q,mi,m2) = 
Hence, we may write 



4-d l^d-U 



El + E2I \ 



Disp/(')(g'i,g2,™i,m2,s) = (27r)' 



1 



E1E2 IT 
1 



^'0 



El + E2 



. (F.20) 



qi - q2je 



(</! - Q2) + 

\^ ias,'^^^,mi,m2} , 



(F.21) 



where, for d = 4 — 2e dimensions, 



7E + In 



+ X^^'^ {{q^ — ias)'^ — ,1711, m^) cosh 



mima (^o _ ^q^s) - q 

I f ml + 1712 — {q^ — iasY + q^ 
2mim2 



(2 2\^ ^i 

— rrii) m — ^ 



, (F.22) 
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containing the zero-temperature UV divergence [84]. In (F.22), 7e is the Euler constant 
and is the 't Hooft mass scale. 

(ii) The purely thermal part is more straightforward to analyze: 



d'^ki d% d^k2 



in^ (27r)4 (27r)4 (2^)4 
X (27r)35(3)(qi - ki - k2)(27r)35(3)(q2 _ j^/ _ ^) 

X As\(q^, - ql - k\ + - kl + + 4^2" 

'ql + ql-k\-k[' -kl-k',')' +As'' 

X 2^<5(fci - m2)|2A;2T/'/2^(A:2,A:^,t)e^('=°-^^>'^l2fc^Y/'2vr5(A;^^ 



mi 



(F.23) 



Performing the four zeroth-component momentum integrations by virtue of the on-shell 
delta functions 6{kf — mf) and 5{k'-^ — mf) in (F.23), we find 

/(")(gi,(?2,mi,m2,s) = 8(27r)V-''e*^''°~^"^*^E / J dUf-' dU^, dU^^ dU^>^ 

{a} 

X (27r)35(3)(qi - ki - k2)(27r)3<5(3)(q2 - k'^ - k'2) 

X 2 ' 

Qi - q2 - OLiEi + a^i?! - 02-^2 + a2-E^2 ) + 



2is 



X 2^i/i(fci, fc;,t)e*("i^i-°'i-^i)*-f 2,^2/2^(^2, ^2, i)e^^°''^'""^^ • (F.24) 

Note that the purely thermal part in (F.24) contains only absorptive contributions. 

(iii) The cross-terms yield both dispersive and absorptive contributions. The first of the 

cross terms yields 



d'^ki d% d^k2 d^k'2 



iiT^ (27r)4 (2^)4 (27r)4 
X (27r)35(3)(qi - ki - k2)(27r)35(3)(q^ _ ^ _ 



X 4s 



q',-q'2-kUk[-k'2+k'2') +4s 

2 



qO + gO _ ^0 _ 



/O 



+ 4s" 



1 



A;| — + ie 



{2Tr)^6^^\ki - k[) 



X 27r5(fci - m2)|2A:2T/¥2^(fc2,fc^,t)e^(^§-^^>'/|2fc^Y/'2vr<5(fc^^ - ml). 

(F.25) 

After evaluating the k[^, k2 and ^2'^ integrals by virtue of the delta functions, we perform 
the ki integral by closing a contour in the lower half of the k^ complex plane. Equa- 
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tion (F.25) may then be written in the more compact form 



{a} 



X (2^)35(3)(qi - ki - k2) (2^)35(3) (q2 - k[ - k',) 



1 



X 



vr 

X ai 



Qi - q2 - "2-^2 + a'sE^sj + 4s 

9 9 '-^l I ^1 



IS 



X 



(F.26) 



Separating again into dispersive and absorptive parts as in (F.18), may be written as 



/('"")(gi,g2,mi,m2,s) = 8(2vr)V-"e^^'^^-'^°^*^E / / dH^^^^ dHk; dHk, dHk: 



{a} 



X (27r)3<5(3)(qi - ki - k2)(2^)35(3)(q^ _ _ 



qi - ql - a2^2 + "2^2 ) + 



-1 



aii?i 



02^2 + aLEL ^ 2 



+ is 



(F.27) 



Similarly, for the second cross term 

j(iiib) ;^2) , we obtain 



I^"'''Hqi,q2,mi,m2,s) = 8(27r) V"''e'('^° " 

{a} 

X (27r)35(3)(q^ - ki - k2) (2^)35(3) _ - k'2 



dn^^^Mnk^dnk, dn^^ 



1 

X — 
vr L 



qi - q2 - aiEi + a[E[ j + 

{-a2(^-^^i^^^-a2i^2) +^s} 
(27r)32^^5(3)(k2 _ k'2)2(fi/i(^i,fc;,t)e^(°i-^i-"'i-^i)*/ . 



(F.28) 



Collecting the four contributions from (F.18), (F.24), (F.27) and (F.28) together, we 
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find the full time-ordered function 



Bo{qi,q2,mi,m2,s) = 8(27r)V^ 



^3„4-d,.(,o _ ,0),-^, ^ / . . . / dn^;i dn^ dHk, dn^. 



1 



vr 



Q'? - ^2 ~ ai-Ei + o'l-Ei - a2E2 + a'2-E'2 ) + 4s 



g? + g^' _ aiEi+a[E[ _ 02 £^2 + ^2^2 " ^ 
2 2 2 ' ' ^ 



2 2 2 



(F.29) 



where {k} = {ki, k[, k2, k2} is the set of on-shell four momenta. The statistical structure 
is contained within the distributions F^{{k},t) and F^{{k},t) defined as follows: 



F^{{k},t) = {2TTf5^^^\ki - k[){2TTf6^^\k2 - k'2) (F.30a) 
X (^9{ai, a'l, 02, 02) ~ ^(— ai, — a'^, — 02, — 02) ) 
+ (2^)V3)(ki - k[)(^eiai,a[) - ei-ai,-a[))f^ik2,k'2,t) 
+ h{k,k'i,t)i27rf5(-'Hk2 - k'2) (0(02, 4) - e(-a2,-a'2)) , 
F\{k},t) = {27rf5'^^\ki - k[){27rf6^^\k2 - k'2) (F.30b) 
X (^9{ai,a[, 02,0(2) + "i> — a'l) — "2, — "2) ) 
+ (2^)V3)(ki - k[)(^e{ai,a\) + e{-ai,-a\))f^{k2,k'2,t) 
+ Mki,k[,t){27Tf6'^''Hk2 - k'2) (0(02, 4) + ^(-a2,-a'2)' 
+ 2fi{ki,K,t)f^{k2,k'2,t) . 



For completeness, the anti-time-ordered function, obtained from the RHS (F.29) by 
taking s — >• —s and multiplying by an overall minus sign, is given by 



B*,iq^,q2,m,,m2,s) = 8(27r) V'e*^''" " ^ 

{a} 

X (2^)35(3)(qi - ki - k2)(27r)3<5(3)(q2 - k'^ - k'2) 
X 2S'ie'^"^^^ ~ "i^0*/2(f2e*^"2'^^ " o^2^2)*/ 



dn^^Mn^dHk, dHk^ 



1 



vr 



Qi - - aiEi + a[E[ - a2E2 + aa^^a ) + 



{( 



2 2 2 

gg + g§ _ QiEi +a;g( _ Q2E2 - 

2 2 



-1 



^i^)F^^(W,t) +zsFi(W,t)}, (F.31) 
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Finally, from (4.27a), the 'Hadamard' function Bq is given by 



Bl{qi,q2,mi,m2,s) = 16(2^) V~''e'^'^^ " "^^^^^ J] 

{a} 

X (2^)35(3)(q^ _ _ k2)(2vr)35(3)(q2 - k'^ - k'^ 



dnf' dn^ dHk^ dHk^ 



X — 

TT L 



Qi - - + - «2^2 + a2£^2 ) + 4s 



2 2 



isF\{k},t), (F.32) 



which depends only on F^{{k},t), as one should expect. 



Absolutely-Ordered Functions 

We turn our attention now to the (a, 6) = (2,1) and (a, 6) = (1,2) elements of (F.2). 
These are the positive- and negative-frequency absolutely-ordered functions, respectively: 



B^''^\qi,q2,mi,m2,s) = -8(27r/i)* "^e 



d'^ki d'^k[ d^k2 d% 



in^ (2^)4 (27r)4 (27r)4 
X (2^)35(3)(q^ _ _ k2)(2^)3<5(3)(q2 - - k'2) 



X 4s 



q'l-ql-kl + k',' -kl + k',') +4.2 



X 27T6{kl - ml)\2k^^\^/^(^e{{-)k°^,{-)k[°){2TTf5^^\ki - k[) 
+ ^ik^,k[,t))e^('''^-''r)if\2k[Y'2n6{k[' - m?) 

X 2^5(A;i - ml)\2k',\'/^[e{{-)kl{-)k',''){27rf6^'\k2 - k'2) 

+ f^ik2,k'„t))e^(''°^-''^"^'f\2k^2\'^'27r5{k',' - ml) , (F.33) 



where 6{x,y) is defined in (F.ll). 

After carrying out the zeroth-component momentum integrals, we obtain 



^0^^^^(91,92,^1,^2,5) = 16(27r)V~'*e 



{a} 



■■■ /dn^^Mn^dnk.dnk^ 



X {2TTf5^'\^i - ki - k2)(2^)35(3)(q2 - k'l - k'2) 



1 r 



vr 



- - Qi^i + a'i£^; - a2E2 + a'^E'^ ) + 4s 



2 2 



+ S^ 



isF>«)({^},t) , (F.34) 
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where 

F>«\{k},t) = (0((-)ai,(-)a;)(2^)¥3)(ki - k[) + fiih,k[,t)) 

X (0((-)a2,(-)a'2)(27r)3<5(3)(k2 - k'^) + f^{k,k'2,t)) ■ (F.35) 



We may confirm the following relations between the -F's: 

F\{k},t) = F>i{k},t) - F<i{k},t) , (F.36a) 
F\{k},t) = F>i{k},t) + F<i{k},t) , (F.36b) 

which are in line with our notation for the propagators and the self-energies. 

Causal Functions 

We are now in a position to obtain the causal counterparts. Given the relations in (A.6d), 
the retarded and advanced functions B^^^^ are given by 

Bf^^\q,,q2,m,,m2,s) = -8(27r) V-'e*^"^" " E /••• J dU^-;' dUy^ dU^, dUy^ 

{a} 

X (2^)35(3) (q^ - ki - k2)(27r)35(3)(q^ _ j^/^ _ ^) 



1 

X — 
vr 



Qi - q2 - aiEi + aiE[ - 02-^2 + "2-^2 ) + 4s 



„2 



-1 



X 



l+£ _ a,E,+a[Ei _ a,E, + a',E', +(_),,^ , (F.37) 



where = F^{{k},t). 

The Thermodynamic Equilibrium Limit 

In the thermodynamic equilibrium limit, we expect to be able to recover the results from 
the discussions of Section 5, using the correspondence identified in (5.8). It follows that 
the various distribution functions satisfy the following factorization: 



F{{k},t) ^ (27r)35(3)(ki - k[)(^9{ai,a\) + e{-ai,-a[)j (F.38) 
X (2^)V=^)(k2 - k'2)f 0(02,02) + 9{-a2,-a2))aia2F^^{aiEi,a2E2) , 



where 



F>(ai^i,a2^2) = 1 + /B(ai^i) 1 + 7^(02^2) , (F.39a) 



cq 



F<{aiEi,a2E2) = /B(ai^i)/^(a2^2) • (F.39b) 

Using the above expressions, we may perform all but one of the three-dimensional phase- 
space integrations and both of the 'primed' summations in the various Bq functions. We 
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then obtain the following set of 'equilibrium' results: 



1 



1 



Bo{-B^){qi,q2,mum2,s) = (2^)*- , . n^z 

^ {li -92)+ 
1 aiQ!2 r/ „o 4_„o 



i(9?-9§)i75(3)(, 



qi - q2j 



- aiEi - 02^2 + s 



-1 



(+) 



TT E1E2 

2 ai^i - 02-^2 ) -f^eq(«i-^i> "2-^2) + isFl^{aiEi,a2E2) 



(F.40a) 



-S5''^'^(g'i,g'2,?n'i,'m-2,s) = (27r)*- 



4I L e^('/?-92°)t'/5(3)(q^ _ q^) 



X 2isF>^<'^{aiEi,a2E2) 
B^^^\qi,q2,mi,m2,s) = -(27r)'^- 



-1, 1 aia2 



TT E1E2 - 



TT (nO - q^Y + 4s2 



- aiEi - 02^2 ) + s" 



eiii1-i'2)ifs(^){ci, - q2) 



(F.40b) 



{a} 

X Ff^^^\aiEi,a2E2) , (F.40c) 

with Fc'^(aiSi,a2£;2) = Fg^(aiSi, 02^2)- 

At this point, we caution the reader that (F.40a)-(F.40c) are not the exact equilibirum 
-Bo functions. For late times, we may use the final value theorem 

lim F{t) = lim sCAfMs) (F.41) 



i -> 00 ~ ' s — 5> 

to obtain the true time-invariant equilibrium functions: 



i?o(-i?o)(9i,92,mi,m2) = {2TT)H^''\q^ - g2)(2vr^)'-'^ ^ /d'^'ik 



1, aia2 

1 

E1E2 



+ 



(01^1,02^2) 



2tt q^ — aiEi — 02-^2 



+ 7,^{Qi - - a2E2)F^ JaiEi,a2E2) ] , 



(F.42a) 



aia2 
E1E2 



B>'<'\qi,q2,mi,m2) = i(27r)V^)(gi - (?2)(27r//)4-^ J] /d'^-iki 

X (5(g? - ai^i - a2-E2)-Fcq'^'^(ai^i,a2F2) , (F.42b) 
B^^^\quq2,m,,m2) = -(27r)V^)fe - g2)(2M'-'E /d'^'^i "^"^ 



{a} 



E1E2 



Ff {aiEi,a2E2) . 



2tt q\ - aiEi - a2E2 + (-) ie ^ 



(F.42c) 



The above expressions are consistent with known results calculated in the ITF or equilib- 
rium CTP formalism (see Appendix B). 
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The Homogeneous Limit 

Lastly, we summarize the time-dependent homogeneous limit of the Bq functions. In the 
spatially homogeneous case, the -F's satisfy the following factorization 

F{{k},t) {2TTf6^^\'ki - k\)(^e{ai,a[) + 9{-ai,-a\)^ (F.43) 
X (27r)35(3)(k2 - k'2)(0(a2,a2) + ^(-02, -aa) ) i^hom(ai, "2, ki, k2) . 
The set of homogeneous distributions Fhom can be obtained from 

F^}^\ai,a2MM,t) = (^((-)«i) + e{ai)f{\ki\,t) + 0(-ai)/^(|ki|, t) ) (F.44) 

'eii-)a2) + e{a2)fi\k2\,t) + 0(-a2)/(|k2|,t)) , 



using the relations in (F.36). Substituting these distributions into the non-homogeneous Bq 
functions, we perform all but one of the three-dimensional phase-space integrals and both 
'primed' summations. Subsequently, with the aid of an inverse Laplace transformation 
with respect to s, the following set of t-dependent homogeneous Bq functions are obtained: 



Bo{-B^){qi,q2,mi,m2,if;U) = (27r)V3)( 



qi - q2 



(F.45a) 



X ^ / d'^-iki 

{a} 



E1E2 



9? + 92° 



+ 



aiEi — Q2-E2 



27r (gO - aiEi - 02-^2) (^2 ~ ctiEi - a2-E'2) 



X 



{^t{qi - ql) - 6t[ql + ql - 2aiEi - 202^2) ) Fifom(«l, «2, ki, qi - ki, t) 
^St{qi - aiEi - a2E2)5t{ql - aiEi - a2-E2)Fhom(ai, "2, ki, qi - ki, t) 



B^^<^\qi,q2,mi,m2,tr,ii) = i(27r)V3)(qi - q2)(2^/i)4-V('?i " "2)*/ 

X ^ y '^'^'^^^^[^jMi - - (^2E2)5t[q^ - aiEi - 02^2 



{a} 

^ -^h^m'^(«i'«2,ki,qi - ki,t) 



(F.45b) 



Bf^^\qi,q2,mi,m2,if,ii) = {27t)^6'^''\ 



X ^ /d'^-iki 



E1E2 



qi - q2)(2^^)^-V('??-''2)*^ 

1 - aiEi - a2E2 

27r (qO - aiEi - 02-^2) (^'2 - "1-^1 - a2-£^2) 



X 



Stiq^i - Q2) - HQi + Q2 - 2ai^i - 202^2) 



+ (-) i^HQi - ai^i - a2^2)'^i(g2 - aiEi - 02^2) 

^ -^h!5m^(«i'«2,ki,qi - ki,t) . 
In the limit t — )• 00, using (3.51) and making the replacement 

^hom(ai,a2,ki,qi - ki,t) aia2-?^eq(ai, "2, ki, qi - ki) , 
we recover the time-independent equilibrium results in (F.42). 



(F.45c) 
(F.46) 
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